EDGEWOOD 
RESEARCH, 
DEVELOPMENT  & 
ENGINEERING 
CENTER 


ERDEC-CR-019 


AEROSOL  SAMPLING  MODELS  SURVEY 


Dennis  Metz 
Paul  Harvey 


GENERAL  MANAGEMENT  ASSOCIATES 

Abingdon,  MD  21009 


March  1993 


Approved  for  public  release;  distribution  is  unlimited. 


6  i ..  00  ^ 


93 


93-13574 


f 


Disclaimer 

The  findings  in  this  report  are  not  to  be  construed  as  an 
official  Department  of  the  Army  position  unless  so  designated  by 
other  authorizing  documents. 


REPORT  DOCUMENTATION  PAGE 


form  ApprOiPii 
OMB  No  O.WA  0!BS 


I’.jn'i'-  r«*o  'f  rm-j  0^  mfo^n'ation  tO  -^.Pf  •»•■}<*  ’  D*”  r-  ru^jir.-}  thr  'f'.if.N  ■ti'j  ^''Sf  f  u'  T*''’" 

jnhpf'ng  ir'd  m->inr  MfMn<)  Ih#»  need»vd.  ,ind  CCfT’ptet'f’9  ni  '.•nd  ^cmfr>«vr,i^  rog.ug.ri-^  t;  •' 

f  tion  ‘>!  irifo^i^jnon.  -nc  (tiding  iug'3**<tton^  (Of  f^dufing  this  bufOpn  to  A^^^h.ngfQn  >*»*  -ups.  '  ■•‘^^s'rr^uon  UfM*'  - 

f>,.s  l»'g(  ,Mv,  "iuitP  WC'1  AiJuKjfon.VA  ^10?  .^nd  to  t^i*  0((*'**  ■'(  -ind  Pud  1P»  ntpprA.rfi  H^du' t  d ’(toy  .‘.A 


3  REPORT  TYPE  AND  DATES  COVERED 

Final,  90  Jun  -  90  Dec 


5  FUNDING  NUMBERS 

C-DAAA15-87-D-0021 

TA-99 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

ERDEC-CR-019 


10.  SPONSORING /MONITORING 
AGENCY  REPORT  NUMBER 


1.  AGENCY  USE  ONLY  (leave  blank)  2.  REPORT  DATE 

1993  March 


A  TITLE  AND  SUBTITLE 

Aerosol  Sampling  Models  Survey 


6.  AUTHOR(S) 

Metz,  Dennis,  and  Harvey,  Paul 


7.  PERFORMING  ORGANIZATION  NAME<S)  AND  AODRESS(ES) 

General  Management  Associates 
1308  Continental  Drive,  Suite  K 
Abingdon,  MD  21009 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AOORESS(ES) 

DIR,  ERDEC,*  ATTN:  SCBRD-RTB, 

APG,  MD  21010-5423 


11.  SUPPLEMENTARY  NOTES 

COTR:  Robert  Doherty,  SCBRD-RTB,  (410)  671-2326 

*When  this  work  was  performed,  ERDEC  was  known  as  the  U.S.  Army 
Chemical  Research,  Development  and  (Continued  on  page  2) 


12a.  DISTRIBUTION /AVAILABILITY  STATEMENT  12b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  is 
unlimited. 


13.  ABSTRACT  (Maximum  200  words) 

The  Obscuration  Sciences  Branch  of  Physics  Division,  Research 
Directorate  at  the  U.S.  Army  Chemical  Research,  Development  and 
Engineering  Center,  has  undertaken  a  research  project.  This  effort 
is  aimed  at  identifying  a  comprehensive  algorithm,  computer  program, 
or  theory  that  might  already  exist  and  could  be  used  to  design 
aerosol  sampling  and  transport  systems.  Initially,  a  literature 
survey  was  performed  by  the  Obscuration  Sciences  Branch  which  began 
with  several  hundred  potential  articles  and  has  been  distilled  to 
four  candidate  reports. 

This  report  evaluates  the  four  candidate  aerosol  sampling  models  for 
self  consistency  and  compares  the  numerical  calculations  resulting 
from  four  models.  Finally  a  comparison  of  these  model-derived 
results  is  made  to  the  experimentally-derived  results. 


14.  SUBJECT  TERMS 

Aerosol  science 
Aerosol  particles 
Monodispersed  aerosols 


Polydispersed  aerosols 
Aerosol  sampling  train 
(Continued  on  page  2) 


15.  NUMBER  OF  PAGES 

149 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION  118.  SECURITY  CLASSIFICATION  119.  SECURITY  CLASSIFICATION  I  20.  LIMITATION  OF  ABSTRACT 

OF  REPORT  I  OF  THIS  PAGE  I  OF  ABSTRACT  I 


UNCLASSIFIED 


UNCLASSIFIED 


UNCLASSIFIED 


NSM  7540  01  780-5500 


Standard  form  298  (Rev  2  89) 

pTf*scTtb*’d  bv  ansi  S’d 
;'J8  ro.? 


11.  SUPPLEMENTARY  NOTES  (Continued) 

Engineering  Center,  and  the  Contracting  Officer's  Representative 
was  assigned  to  the  Research  Directorate. 

14.  SUBJECT  TERMS  (Continued) 

Aerosol  sampling  models 
Experimental  verification 
Aerosol  penetrating  efficiency 
Aerosol  sampling  and  transport  systems 


2 


PREFACE 


The  work  described  in  this  report  was  authorized  under 
Contract  No.  DAAA15-87-D-0021,  Task  99.  This  work  was  started  in 
June  1990  and  completed  in  December  1990. 

The  use  of  trade  names  or  manufacturers'  names  in  this 
report  does  not  constitute  an  official  endorsement  of  any 
commercial  products.  This  report  may  not  be  cited  for  purposes 
of  advertisement. 

Reproduction  of  this  document  in  whole  or  in  part  is 
prohibited  except  with  permission  of  the  Director,  U.S.  Army 
Edgewood  Research,  Development  and  Engineering  Center  (ERDEC),* 
ATTN:  SCBRD-RT,  Aberdeen  Proving  Ground,  MD  21010'-5423. 

However,  the  Defense  Technical  Information  Center  and  the 
National  Technical  Information  Service  are  authorized  to 
reproduce  the  document  for  U.S.  Government  purposes. 

This  report  has  been  approved  for  release  to  the 

public. 


Acknowledgments 

The  authors  express  their  appreciation  to  Robert 
Doherty,  ERDEC,  for  his  support  and  guidance  during  the  effort. 
In  addition,  appreciation  is  given  to  Dr.  Andrew  McFarland  and 
his  staff  from  the  Department  of  Mechanical  Engineering  at  Texas 
A&M  University  (College  Station,  TX)  in  performing  the  numerical 
model  calculations  and  conducting  the  wind  tunnel  experiments; 
Dr.  Vladimir  Kogan,  Battelle  (Columbus,  OH)  in  performing 
numerical  model  calculations;  and  Dr.  Narayanan  Rajendran,  IITRI 
(Chicago,  IL)  in  performing  numerical  model  calculations. 


*When  this  work  was  performed,  ERDEC  was  known  as  the  U.S.  Army 
Chemical  Research,  Development  and  Engineering  Center,  and  the 
Contracting  Officer's  Representative  was  assigned  to  the 
Research  Directorate. 


3 


CONTENTS 


Page 

SECTION  1  -  INTRODUCTION 

1.1  Background .  11 

1.2  Objective .  11 

1.3  Study  Approach .  12 

SECTION  2  -  AEROSOL  SAMPLING  MODEL  DESCRIPTION 

2.1  Introduction  .  . .  13 

2.2  Numerical  Model  Developed  by  Texas  A&M 

University .  13 

2.3  Numerical  Model  Developed  by  IITRI  .  18 

2.4  Numerical  Model  Developed  by  Battelle  .  19 

2.5  Numerical  Model  Developed  by  University  of 

Duisburg .  20 

2.5.1  Calculation  of  the  Sampling  Efficiency  .  ,  21 

2.5.2  Calculation  of  Transport  Efficiency  Through 

a  Straight  Tube .  21 

2.5.3  Calculation  of  Transport  Efficiency  Through 

a  90®  Bend  24 

SECTION  3  -  AEROSOL  SAMPLING  TRAIN  CONFIGURATION  FOR 
NUMERICAL  CALCULATIONS  AND  EXPERIMENTAL 
VERIFICATION 

3.1  Physical  Description  .  27 

3.2  Region  of  Aerosol  Sampling  Train  Addressed  by 

Four  Models .  27 

3 . 3  Parameters  of  Model  Calculations  and  Experimental 

Verification  .  27 

SECTION  4  -  RESULTS  OF  MODEL  NUMERICAL  CALCULATIONS 

4.1  Results  from  Texas  A&M  Model .  33 

4.1.1  Results  of  Monodispersed  Simulation  Cases  .  33 

4.1.2  Results  of  Polydispersed  Simulation  Cases  .  33 

4.2  Results  from  IITRI  Model .  37 

4.3  Results  from  University  of  Duisburg  Model  ...  40 

4.3.1  Background .  40 

4.3.2  Results  of  Monodispersed  Simulation  Cases  .  40 

4.4  Results  from  Battelle  Model  .  44 

4.4.1  Results  of  Monodispersed  Simulation  Cases  .  44 

4.4.2  Results  of  Polydispersed  Simulation  Cases  .  44 

4.5  Comparison  of  Numerical  Calculations  .  48 

4.5.1  Numerical  model  calculations  for  Inlet 

Region  of  Aerosol  Sampling  Train  .  48 


5 


4.5.2 


Numerical  Efficiency  Calculations  for 
Horizontal  Region  of  Aerosol  Sampling 
Train .  52 

4.5.3  Numerical  Efficiency  Calculations  for  90® 

Elbow  Region  of  Aerosol  Sampling  Train  .  .  58 

4.5.4  Numerical  Efficiency  Calculation  for 

Vertical  Region  of  Aerosol  Sampling 

Train .  64 

4.5.5  Numerical  Efficiency  Calculations  for  45° 

Elbow  Region  of  Aerosol  Sampling  Train  .  .  73 

4.5.6  Numerical  Efficiency  Calculations  for 

Inclined  Region  of  Aerosol  Sampling  Train  .  79 

4.5.7  Cumulative  Efficiency  Calculations  for 

Overall  Sampling  Train  .  85 

SECTION  5  -■  RESULTS  OF  EXPERIMENTAL  VERIFICATION 

5.1  Experimental  Apparatus  and  Methodology  for 

Aerosol  Penetration  Measurements  .  91 

5.2  Experimental  Results  .  93 

SECTION  6  -  COMPARISON  OF  EXPERIMENTAL  AND  NUMERICAL 
RESULTS 

6.1  Limitations  in  Comparison  of  Experimental  and 

Numerical  Results  .  97 

6 . 2  Discussion  of  Experimental  vs  Numerical 

Results .  97 

G . 2 . 1  Comparison  of  Results  for  0°  Inlet 

Orientation  and  70  1/min  Flow  Rate 
Conditions .  97 

6.2.2  Comparison  of  Results  for  0°  Inlet 

Orientation  and  130  1/min  Flow  Rate 
Conditions . loi 

6.2.3  Comparison  of  Results  for  90°  Inlet 

Orientation  and  70  1/min  Flow  Rate 
Conditions . 101 

6.2.4  Comparison  of  Results  for  90°  Inlet 

Orientation  and  130  1/min  Flow  Rate 
Conditions . 102 

SECTION  7  -  CONCLUSIONS  AND  RECOMMENDATIONS 

7.1  Conclusions . 103 

7.1.1  Capabilities  of  Numerical  Models/ 

Algorithms . 103 


6 


7.1.2  Numerical  Calculations  Generated  by  Four 

Models . 103 

7.1.3  Comparison  of  Experimental  Measurements 

and  Numerical  Calculations . 104 

7.2  Recommendations . 104 

REFERENCES .  107 

APPENDIX  A . . 109 

APPENDIX  B .  117 

APPENDIX  C . . . 127 

APPENDIX  D .  131 

APPENDIX  E .  137 

APPENDIX  F .  14  7 


LIST  OF  TABLES  AND  FIGURES 

Table 


1  Four  Candidates  Models'  Characterization/ 

Simulation  of  Aerosol  Sampling  Train  .  29 

2  Matrix  of  Simulation  Cases  to  be  Generated 

for  Aerosol  Sampling  Train  Configuration  ....  31 

3  Cumulative  Sampling  Efficiency  Results  by 

Texas  A  &  M  Model .  34 

4  Cumulative  Sampling  Efficiency  Results 

Generated  by  IITRI  Model  .  38 

5  Cumulative  Sampling  Efficiency  Results 

Generated  by  University  of  Duisburg  Model.  ...  41 

6  Cumulative  Sampling  Efficiency  Kesults 

Generated  by  Battelle  Model  .  45 

7  Numerical  Model  Calculations  (Monodispersed 

Particle  Sizes)  for  Inlet  Region  of  Aerosol 
Sampling  Train  .  49 

8  Model  Efficiency  Calculations  (Polydispersed 

Particle  Sizes)  for  Inlet  Region  of  Aerosol 
Sampling  Train  .  53 

9  Model  Efficiency  Calculations  (Monodispersed 
Particle  Sizes)  for  Horizontal  Region  of 

Aerosol  Sampling  Train  .  56 

10  Model  Efficiency  Calculations  (Polydispersed 
Particle  Sizes)  for  Horizontal  Region  of 

Aerosol  Sampling  Train  .  59 

11  Numerical  Model  Calculations  (Monodispersed 
Particle  sizes)  for  90®  Elbow  Region  of 

Aerosol  Sampling  Train  .  62 


7 


12  Model  Efficiency  Calculations  (Polydispersed 
Particle  Sizes)  for  90°  Elbow  Region  of 

Aerosol  Sampling  Train  .  65 

13  Numerical  Model  Calculations  (Monodispersed 
Particle  Sizes)  for  Vertical  Region  of 

Aerosol  Sampling  Train  .  68 

14  Model  Efficiency  Calculations  (Polydispersed 
Particle  Sizes)  for  Vertical  Region  of 

Aerosol  Sampling  Train  .  70 

15  Numerical  Model  Calculations  (Monodispersed 
Particle  Sizes)  for  45°  Elbow  Region  of 

Aerosol  Sampling  Train  .  74 

16  Model  Efficiency  Calculations  (Polydispersed 
Particle  Sizes)  for  45°  Elbow  Region  of 

Aerosol  Sampling  Train  .  76 

17  Numerical  Model  Calculations  (Monodispersed 
Particle  Sizes)  for  Inclined  Region  of 

Aerosol  Sampling  Train  .  80 

18  Model  Efficiency  Calculations  (Polydispersed 

Particle  Sizes)  for  Inclined  Region  of  Aerosol 
Sampling  Train  .  82 

19  cumulative  Sampling  Efficiency  Results  Generated 

by  the  Three  Aerosol  Models .  86 

20  Cumulative  Penetration  Efficiency  Results 

Derived  Experimentally  by  Texas  A  &  M .  94 

21  Cumulative  Penetration  Efficiency  Results: 

Model  Predictions  vs.  Experimental  Data .  98 

Figure 

1  Schematic  of  Aerosol  Sampling  Train  .  28 

2  Regions  of  Sampling  Train  Characterized/ 

Simulated  by  Each  Candidate  Model  .  30 

3  Penetration  Efficiency  -  Texas  A  &  M  Model 

Monodispersed  Results  .  35 

4  Penetration  Efficiency  -  Texas  A  &  M  Model 

Polydispersed  Results  .  36 

5  Penetration  Efficiency  -  ITTRI  Model 

Monodispersed  Results  .  39 

6  Penetration  Efficiency  -  Duisburg  Model 

Monodispersed  Results  .  42 

7  Penetration  Efficiency  -  Duisburg  Model 

Polydispersed  Results  .  43 

8  Penetration  Efficiency  -  Battelle  Model 

Monodispersed  Results  .  46 

9  Penetration  Efficiency  -  Battelle  Model 

Polydispersed  Results  .  47 


8 


10  Inlet  Pene^'.ation  Efficiency  -  Monodispersed 

Result?  -  0°  Inlet  50 

11  Inlet  penetration  Efficiency  -  Monodispersed 

Results  -  90°  Inlet  51 

12  Inlet  Penetration  Efficiency  -  Polydispersed 

Results  0°  Inlet  . 54 

3  Inlet  Penetration  Efficiency  -  Polydispersed 

Results  90°  Inlet  55 

14  Penetration  Efficiency  -  Monodispersed  Results  - 

Horizontal  Sec . 57 

15  Penetration  Efficiency  -  Polydispersed  - 

Horizontal  Sec.  -  0°.  60 

16  Penetration  Efficiency  -  Polydispersed  - 

Horizontal  Sec.  -  90°  61 

17  Penetration  Efficiency  -  Monodispersed  Results 

90°  Elbow  63 

18  Penetration  Efficiency  -  Polydispersed  -  90° 

Elbow  -  0°  Inlet  66 

19  Penetration  Efficiency  -  Polydispersed  -  90° 

Elbow  -  90°  Inlet  . 67 

20  Penetration  Efficiency  -  Monodispersed  Results 

Vertical  Sec  . 69 

21  Penetration  Efficiency  -  Polydispersed  - 

Vertical  -  0°  Inlet  71 

22  Penetration  Efficiency  -  Polydispersed  - 

Vertical  -  90°  Inlet  72 

23  Penetration  Efficiency  -  Monodispersed  Results 

45°  Elbow  . 75 

24  Penetration  Efficiency  -  Polydispersed  -  45° 

Elbow  -  0°  Inlet  77 

25  Penetration  Efficiency  -  Polydispersed  -  45° 

Elbow  -  90°  Inlet  78 

26  Penetration  Efficiency  -  Monodispersed  Results  - 

Incline .  81 

27  Penetration  Efficiency  -  Polydispersed  - 

Incline  -  0°  Inlet .  83 

28  Penetration  Efficiency  -  Polydispersed  -  Incline  - 

90°  Inlet  . 84 

29  Penetration  Efficiency  -  Monodispersed  - 

Cumulative  -  0°  Inlet  87 

30  Penetration  Efficiency  -  Polydispersed  - 

Cumulative  -  0°  Inlet  89 

31  Experimental  Results  -  Cumulative  Penetration 

Efficiency . 95 

32  Experimental  Results  -  Cumulative  Penetration 

Efficiency . 96 


33  Penetration  Efficiency  -  70  L/min,  0°  Inlet  .  .  99 

34  Penetration  Efficiency  -  130  L/rain,  0°  Inlet  .  .  100 

9 


AEROSOL  SAMPLING  MODELS  SURVEY 


SECTION  1  -  INTRODUCTION 


1 . 1  BACKGROUND 

Aerosol  sampling  and  transport  systems  are  extensively 
used  for  various  applications  that  include  measuring  airborne 
particle  size  distributions  and  concentrations.  Characterization 
of  obscuration  smokes  is  an  example  of  such  a  use  (Ref.  1), 

In  aerosol  measurements,  as  in  many  other  measurements, 
it  is  often  not  feasible  to  perform  an  in  situ  measurement. 

There  is  a  spatial  distance  between  the  site  of  measurement  and 
the  instrument  in  which  the  actual  analysis  is  performed.  In 
order  to  bridge  this  gap,  a  sample  must  be  taken  since  in  almost 
all  cases  it  is  impossible  to  analyze  the  total  amount  of  aerosol 
(Ref.  2).  The  aerosol  sample  is  then  transported  to  the 
instrument. 

Generally,  aerosol  sampling  devices  include  an  inlet 
sampling  probe  and  a  detector.  Results  of  these  measurements,  as 
observed  by  the  detector,  do  not  usually  reflect  povential  losses 
of  particulate  matter  within  the  probe.  Therefore,  in  many 
cases,  accuracy  of  these  measurements  is  not  known  due  to  the 
removal  of  particles  from  air  streams  in  the  inlet  trains  (Ref. 

1)  . 


The  Obscuration  Sciences  Branch  of  Physics  Division, 
Research  Directorate  at  the  U.S.  Army  Chemical  Research, 
Development  and  Engineering  Center,  has  undertaken  a  research 
project.  This  effort  is  aimed  at  identifying  a  comprehensive 
algorithm,  computer  program,  or  theory  that  might  already  exist 
and  could  be  used  to  design  aerosol  sampling  and  transport 
systems.  Initially,  a  literature  survey  was  performed  by  the 
Obscuration  Sciences  Branch  which  began  with  several  hundred 
potential  articles  and  has  been  distilled  to  four  candidate 
reports . 


1 . 2  OBJECTIVE 

The  objective  of  this  effort  is  two-fold:  the  first  is 
to  evaluate  the  four  candidate  aerosol  sampling  models  for  self 
consistency.  The  second  is  to  compare  the  numerical  calculations 
resulting  from  the  four  models  and  compare  these  model-derived 
results  to  the  experimentally-derived  results. 
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1.3 


STUDY  APPROACH 


The  technical  approach  used  to  accomplish  the  task 
objective  hinged  on  the  implementation  of  four  key  ,teps 
(phases) . 


In  the  first  phase,  information  on  each  of  the  models 
was  reviewed  in  order  to  gain  an  understanding  of  the  theory, 
capabilities,  limitations,  and  input  parameters  to  the  model. 
This  first  phase  involves  a  write-up  describing  each  of  the 
models. 


The  second  phase  required  development  of  an  aerosol 
sampling  train  configuration  upon  which  the  numerical 
calculations  would  be  performed  and  the  experimental  data  would 
be  generated.  This  phase  required  a  write-up  detailing  the 
sampling  train  configuration. 

The  third  phase  required  that  numerical  calculations  be 
made  for  each  of  the  four  models  based  on  the  design  parameters 
of  a  baseline  simulated  set  of  conditions. 

The  fourth  phase  required  a  comparison  of  the  numerical 
calculatior'i  generated  by  the  four  models.  In  addition,  a 
comparison  will  be  made  of  these  model-generated  results  to  the 
experimentally-derived  results.  Findings,  conclusions  and 
recommendations  will  be  compiled  resulting  from  these 
comparisons. 


SECTION  2  -  AEROSOL  SAMPLING  MODEL  DESCRIPTION 

2 . 1  INTRODUCTION 

As  part  of  the  literature  search  conducted  by  the 
Obscuration  Sciences  Branch  to  identify  algorithms,  computer 
programs,  or  theory  which  already  exist  and  could  be  used  to 
design  aerosol  sampling  and  transport  systems,  four  candidate 
models/computer  programs  and/or  theory/algorithms/equations  were 
identified.  The  four  candidates  were  the  following: 

•  A  numerical  model  and  computer  program  developed 
by  Texas  A&M  University,  College  Station,  TX,  to 
predict  particle  deposition  in  aerosol  sampling 
lines  due  to  turbulent  diffusion  and  gravitational 
settling. 

•  An  analytical  model  and  computer  program  (SAMPF) 
developed  by  Battelle,  Columbus,  OH,  to  predict 
aerosol  deposition  in  straight  and  bent  circular 
pipes  under  turbulent  flow  conditions. 

•  A  theoretical  model  and  computer  program  developed 
by  IITRI,  Chicago,  IL,  to  estimate  sampling  errors 
as  a  function  of  various  inlet  geometries  and 
parallel  plates  of  various  openings. 

•  A  report  authored  by  Drs.  Fissan  and  Schwientek, 
University  of  Duisburg,  GE,  which  gives  the 
conditions  for  representative  sampling  in  the  case 
of  a  ducted  aerosol  flow,  as  well  as  equations 
that  allow  a  first  estimation  of  errors. 

A  description  of  the  four  candidate  models  is  presented 
in  the  subsequent  sections  (i.e..  Sections  2.2  to  2.5). 

2.2  NUMERICAL  MODEL  DEVELOPED  BY  TEXAS  A&M  UNIVERSITY 

A  computer-based  model  to  predict  particle  deposition 
in  aerosol  sampling  lines  due  to  turbulent  diffusion  and 
gravitational  settling  was  developed  by  Mr.  N.K.  Anand  and  Dr.  A. 
R.  McFarland  of  Texas  A&M  University  at  College  Station,  TX  (Ref. 
3) .  Models  developed  by  others  for  determining  penetration 
through  tube  bends  are  included. 

In  modeling  the  individual  components  of  the  aerosol 
sampling  train,  consideration  was  first  given  to  the  inlet 
aspiration  ratio  A,  where: 
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Here:  =  aerosol  concentration  in  the  free  stream  and  = 

aerosol  concentration  at  the  entrance  plane  of  the  inlet.  The 
model  of  Vincent  et  al.  (Ref.  4)  was  used  to  compute  the 
aspiration  ratio,  which  gives: 


A  =  1  + 


1 


1 

1  +  1.05-Stk- (cosG  +  4  •  (R-sine)°-^) 


(R-cos©  -  1) 


(2) 


Where:  R  =  W/V,  W  =  airstream  velocity,  V  =  velocity  at  the 

entrance  plane  of  the  inlet,  0  =  angle  between  the  wind  velocity 
vector  and  the  tube  axis  vector  (which  faces  into  the  inlet) ,  and 
Stk  =  Stokes  number.  The  stokes  number  is  given  by: 

C  •  d  •  D  ^  •  W 

Stk  =  -  (3) 

9-m-D^ 

The  additional  parameters  in  the  Stokes  number  are:  C  = 
Cunningham's  slip  correction  (Ref.  5) ,  d  =  particle  density,  D  - 
particle  diameter,  fi  =  air  viscosity,  and  =  inlet  diameter 
(which,  in  this  system,  is  the  same  as  the  tube  inside  diameter) . 


Once  the  aerosol  is  inside  the  inlet  there  are  losses 
in  the  developing  boundary  layer  due  to  gravitational  and 
inertial  forces.  Okazaki  and  Willeke  (Ref.  6)  presented  a 
semiempirical  model  which  predicts  these  inlet  losses.  They  note 
that  supportive  experiments  were  mostly  conducted  with  tube 
Reynolds  numbers  less  than  2000.  Their  model  is: 


P 


1  -  exp  j^'-4.7 


(4) 


where: 


K 


Z  •  Stk  1 

Re°-5  J 


(5) 


The  parameter  Z  used  in  Equation  5  is  the  gravitational  settling 
number  and  is  given  by: 


Z  * 


L*V 


g 


W-D, 


(6) 


where:  L  =  tube  length  over  which  sedimentation  in  the 

developing  boundary  layer  is  important  and  Vg  =  particle 
gravitational  settling  terminal  velocity.  In  turn,  V  is 
calculated  from; 
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(7) 


C-d-D^  -g 

; - 

18  '  fi 

where:  g  =  local  gravitational  constant. 


In  using  the  model  of  Okazaki  and  Willeke,  Texas  A&M 
assumed  that  the  phenomena  were  of  consequence  over  the  first  20 
cm  of  the  first  horizontal  straight  section  of  tubing.  For  the 
remainder  of  this  section  and  the  other  straight  sections,  Texas 
A&M  used  the  model  of  Anand  and  McFarland  (Ref.  3)  which  takes 
into  account  both  gravitational  settling  and  turbulent  deposition 
and  which  is  applicable  to  horizontal,  inclined  and  vertical 
tubes.  The  model  uses  the  expression  for  the  penetration  of 
aerosol  particles  through  a  straight  tube: 


P  =  exp 


7r-D^-V^*L 


Q 


(8) 


where:  =  effective  depositional  velocity  of  aerosol  particles 

to  the  tube  walls,  L  =  length  of  tube  and  Q  =  volumetric  flow 
rate.  The  depositional  velocity  is  comprised  of  that  due  to 
turbulent  diffusion,  that  due  to  gravitational  settling,  v  , 

and  that  due  to  thermal  diffusion,  The  velocity  vectors  for 

thermal  and  turbulent  diffusional  deposition  are  both  normal  to 
the  tube  wall  whereas  the  gravitational  deposition  is 
antiparallel  to  the  gravitational  vector.  In  the  present 
application,  the  depositional  losses  due  to  diffusion  are  of  no 
consequence  and  will  be  neglected.  For  the  model  of  Anand  and 
McFarland  (Ref.  3),  the  effective  depositional  velocity  is 
calculated  from: 


V 


e 


1 


27r 


■  VgeSina)  da 


(9) 


where:  V  =  sin0,  0  =  angle  of  inclination  of  the  pipe 
relative  to  the  vertical  direction,  a  =  angle  in  the  tube  cross 
section  between  a  radius  vector  and  x-axis  of  the  horizontal 
plan.  The  integral  is  subject  to  the  constraint: 

“  '^ge  si"  a  )  >  O  (10) 
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If  the  constraint  is  not  satisfied,  =  O.  This  constraint  is 
necessary  since  otherwise  the  physical  situation  would  be 
equivalent  to  mass  being  transported  from  the  environment  through 
the  top  (relative  to  the  gravitational  vector)  side  of  the  tube 
wall . 


The  Texas  A&M  model  originally  utilized  the  results  of 
Liu  and  Agarwal  (Ref.  7)  to  predict  the  turbulent  depositional 
velocity,  Their  semiempirical  expression  is: 

Vtd  =  (f/2)°-^  (11) 

Here,  f  is  the  friction  factor  of  the  airflow  in  the  tube  which, 
for  a  smooth  wall  pipe,  can  be  represented  by  Blasius'  equation: 


f 


0.316 


4  -  Re 


0.25 


(12) 


The  parameter  is  the  dimensionless  particle  deposition 
velocity  and  is  given  by: 

=  6.9  X  10"'*  for  <  15 

(13) 

=  0.16/  for  >  15 

Here,  is  the  dimensionless  stopping  distance  and  is  given  by: 

T,  =  t-V.2  /  V  (14) 

where:  t  =  particle  stopping  distance,  V,  =  friction  velocity  and 
V  =  kinematic  viscosity  of  air.  In  turn,  these  parameters  may  be 
expressed  as: 


r  =  Vg  /  g  (15) 

V,  =  (f/2)°-5-V  (16) 

where:  V  =  spatial  mean  air  velocity  in  the  tube. 

Recently,  Texas  A&M  has  modified  their  model  to  utilize 
the  predictive  capabilities  of  the  Beal  (Ref.  8)  formulation  for 
turbulent  and  thermal  diffusional  losses  in  straight  tubes - 
Beal's  model  is  much  more  complicated  than  that  of  Liu  and 
Agarwal  and  its  basis  will  not  be  discussed  herein;  however,  for 
particles  in  the  inertial  range,  the  two  models  give  similar 
results.  The  computations  which  were  carried  in  the  present 
study  employed  the  model  of  Beal. 
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Particle  deposition  in  tube  bends  can  be  correlated 
with  the  Stokes  number  (Cheng  and  Wang  (Ref.  9),  Pui  et  al., 
(Ref.  10)).  For  this  task,  Texas  A&M  has  used  the  model  of  Pui 
et  al.,  which  for  90°  bends  is: 

p  =  (17) 

For  45°  bends,  the  numerical  factor  in  the  exponent  of  Equation 
(17)  is  0.482. 

To  calculate  overall  aerosol  penetration  through  the 
tubing  system,  P^,.,  for  monodisperse  aerosols  of  size  Dp.,  Texas 
A&M  used  the  expression  (Ref.  11): 


where:  Pj  .  =  penetration  of  the  ith  particle  size  through  the 

jth  component  (inlet  aspiration,  turbulent  deposition,  etc.). 

Data  are  reported  for  both  the  overall  penetration  and  the 
penetration  for  the  individual  components. 

For  the  calculational  cases  involving  aerosol 
penetration  through  the  individual  tubing  components,  the  order 
of  the  calculations  makes  no  difference  if  the  aerosol  is 
monodisperse;  however,  if  the  aerosol  is  polydisperse,  the 
calculations  must  be  carried  out  for  the  actual  sequence  of 
components  in  the  tubing  system  layout. 

The  calculations  for  polydisperse  aerosols  were 
performed  by  first  subdividing  a  log-normal  distribution  into 
thirty  particle  size  increments.  The  aerosol  penetration  through 
the  jth  component  of  the  system  is  based  on: 

Pj  =  C-  /  C..1  (19) 

where:  Cj  =  relative  aerosol  concentration  at  the  exit  plane  of 

the  jth  component  and  Cj,,  =  relative  aerosol  concentration  at  the 
entrance  of  the  jth  component.  It  was  assumed  that  the  aerosol 
concentration  in  the  free  stream  was  unity.  The  parameter  C-  was 
determined  by  adding  the  concentration  values  associated  witn 
each  of  k  particle  size  increments. 


(20) 


The  change  in  aerosol  concentration  of  the  kth  particle  size 
increment  was  obtained  from  the  appropriate  component  model 
(inlet  aspiration,  elbow,  etc.). 
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2.3  NUMERICAL  MODEL  DEVELOPED  BY  IITRI 

A  theoretical  model  and  computer  program  was  developed 
by  IITRI  (Ref.  12)  to  estimate  sampling  errors  as  a  function  of 
various  inlet  geometries  (e.g.,  circular  tube  of  thick  and  thin 
wall  and  parallel  plates) ,  of  various  openings.  Samplers  whose 
face  is  not  perpendicular  to  the  ambient  flow  are  simulated  by  a 
line  sink/source  in  a  uniform  stream. 

The  computer  program  actually  consists  of  two  separate 
programs  (Ref.  12).  Program  "FLOWFI”  solves  the  fluid  flow  in 
and  around  the  sampling  head  with  circular/parallel  plate 
geometry.  Program  "TRAJEC”  computes  the  particle  trajectories  in 
the  specified  flow  region.  The  IITRI  model  accounts  for  inertial 
and  sedimentation  effects  on  particle  motion.  Types  of  flow 
simulated  include  calm  air,  variable  flow  directions,  and 
turbulent  flow. 

General  equations  of  the  IITRI  model  are  the  following: 

C  Cj  C  SC 

Sampling  efficiency  =  n  =  -  =  -  •  -  =  £•  -  •  - 

Co  c,  c,  u  c. 


Where: 

Cq  =  number  of  particles/unit  volume  in  free  stream 
C  =  actual  number  of  particles  sensed  by  instrument 
Cj  =  concentration  at  inlet  to  probe 
SC  =  loss  of  particles  in  probe  (=  C,  -  C) 

E  =  efficiency  of  capture  of  inlet 
Ujj  =  free  stream  velocity 
U  =  suction  velocity  of  the  sample 

To  determine  the  "C”  terms,  equations  of  motion/particle 
trajectory  equations  are  derived  and  solved,  to  include  solutions 
for  samples  at  various  angles  to  oncoming  flow.  An  iterative 
solution  is  used  to  determine  when  a  particle  in  the  free  stream 
actually  passes  through  the  inlet.  The  equations  are  a  function 
of  Stokes  number,  sedimentation  velocities,  fluid  velocities,  and 
particle  velocities,  all  relative  to  the  orientation  thickness  of 
the  inlet  (i.e.,  the  equations  solved  indicate  whether  a  particle 
goes  into  the  inlet,  hits  the  inlet  wall,  or  goes  "around"  the 
inlet)  . 
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Assumptions  made  for  solving  motion/particle  equations 
are  the  following: 

•  Particles  are  uniformally  distributed  and  move 
with  same  velocity  as  free  stream  when  far  from 
inlet. 

•  Particles  are  spherical  and  do  not  change  in  size 
due  to  agglomeration  or  other  factors. 

•  Particles  are  "small”  in  comparison  to  probe  size 
and  move  as  individual  particles  with  no  hydro- 
dynamic  interactions  among  themselves  or  between 
probe  walls. 

2.4  NUMERICAL  MODEL  DEVELOPED  BY  BATTELLE 

Battelle  had  developed  earlier  an  analytical  model 
which  addressed  aerosol  deposition  in  straight  and  bent  circular 
pipes  under  turbulent  flow  conditions  (Ref.  13) .  In  order  to 
perform  the  model  calculations  required  under  this  effort  for  the 
Obscuration  Sciences  Branch,  Battelle  modified  it's  existing 
computer  models.  The  models  were  modified  to  specifically 
address  aerosol  deposition  in  a  sampling  probe  which  resulted  in 
the  SAMPLE  computer  program  (Ref.  1) .  The  calculation  of 
particle  behavior  in  SAMPLE  is  done  using  turbulent  air  flow 
models  for  straight  pipes  and  round  circular  bends. 

Turbulent  deposition  of  particles  in  straight  pipes  is 
treated  with  the  assumption  that  the  deposition  surface  may  not 
be  ideally  smooth,  and  even  the  particles  deposited  may  result  in 
an  increased  roughness.  The  calculation  of  aerosol  deposition 
from  turbulent  flow  to  a  rough  surface  is  performed  in  SAMPLE 
using  correlations  developed  by  Wood  (Ref.  14) .  According  to 
this  approach,  deposition  of  particles  due  to  molecular  diffusion 
enhanced  by  turbulent  eddies  and  due  to  inertial  effect  are  both 
considered. 

In  bends,  there  exists  a  secondary  flow  that  promotes  a 
more  energetic  hydrodynamic  regime  and  higher  pressure  drop  as 
compared  to  a  smooth  straight  pipe  of  the  same  length.  This 
provides  favorable  conditions  for  the  diffusional  and  inertial 
deposition  of  aerosol  particles  from  the  turbulent  flow  in  a  bent 
pipe. 


The  Battelle  model  uses  an  assumption  that  turbulent 
deposition  velocities  of  suspended  particles  in  a  curved  pipe  and 
in  a  straight  rough  pipe  of  the  same  length  are  equal  if  the 
pressure  drops  across  them,  which  are  caused  by  the  wall 
friction,  are  the  same.  In  other  words,  SAMPLE  calculates 
turbulent  deposition  in  a  round  bend  by  calculating  deposition  in 
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a  section  of  a  pressure-equivalent  straight  pipe  of  identical 
length  and  of  s _ch  a  roughness  that  it  provides  a  identical 
resistance  coefficient. 

In  addition  to  turbulent  deposition,  SAMPLE  considers 
gravitational  settling  of  particles  to  the  horizontal  surfaces  of 
the  pipes.  As  a  first  approximation,  the  Stokes'  formulation  of 
settling  velocity  is  used. 

SAMPLE  calculates  the  steady-state  deposition  and 
transport  of  aerosols  through  the  sampling  system.  It  calculates 
fluid  thermo-  and  hydrodynamics  properties  for  an  arbitrary  set 
of  sample  flow  rates,  temperature,  and  air  humidity,  and  uses 
these  properties  for  evaluating  the  transport  parameters  for 
airborne  particulate  matter. 

Since  aspiration,  or  inlet  efficiency  is  currently  not 
considered  by  the  model,  isokinetic  sampling  is  assumed  in  the 
calculations.  Also,  this  version  of  the  program  does  not 
consider  the  potential  for  resuspension  of  deposited  particles. 

2.5  NUMERICAL  MODEL  DEVELOPED  BY  UNIVERSITY  OF  DUISBURG 

The  Eissan  and  Schwientek  paper  (referred  to  in  this 
report  as  the  University  of  Duisburg  model)  (Ref.  2)  gives  the 
conditions  for  representative  sampling  of  aerosols  for  the  case 
of  sampling  of  an  aerosol  flowing  through  a  duct  using  a  thin 
walled  probe  facing  upstream  (0°  inlet  orientation) .  Equations 
are  presented  which  allow  a  first  estimation  of  errors  caused  by 
anisokinetic  sampling.  In  addition,  equations  are  presented 
which  address  the  particle  deposition  effects  in  the  transport  of 
the  aerosol  to  the  site  of  the  analysis.  The  transport 
deposition  effects  considered,  in  laminar  and  turbulent  flow,  are 
Brownian  diffusion,  gravitational  settling,  and  inertial 
deposition.  Although  these  effects  are  interdependent,  they  are 
treated  separately  in  the  paper. 

The  discussion  of  the  paper  will  be  presented  in  three 
parts.  These  three  parts  relate  specifically  to  different 
regions  of  the  aerosol  sampling  train.  The  three  parts  are  the 
following: 


•  Calculation  of  the  sampling  efficiency 

•  Calculation  of  transport  efficiency  through  a 
straight  tube 

•  Calculation  of  transport  efficiency  through  a  90® 
bend 

A  listing  of  the  computer  program  by  General  Management 
Associates  (GMA)  which  incorporates  and  extends  the  equations 
presented  in  the  Eissan  and  Schwientek  paper  is  attached  at 
Appendix  A. 
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2.5.1 


CALCULATION  OF  THE  SAMPLING  EFFICIENCY 


The  sampling  efficiency,  e^,  is  the  ratio  of  the 
sampled  aerosol  concentration  to  the  ambient  aerosol 
concentration.  This  efficiency  is  a  function  of  the  Stokes 
number,  Stk^,  and  the  ratio  W/V,  where  W  is  the  face  velocity  of 
the  aerosol  and  V  is  the  suction  veloctiy  created  by  the  sampling 
probe.  For  sampling,  the  Stokes  number  is  given  by  the  following 
relation. 


C-d-D  2  -  W 

Stky  =  -  (22) 

9-m-D, 

where  W  =  Face  velocity 

Dg  =  Inside  diameter  of  probe 
M  =  Dynamic  viscosity  of  fluid 
d  =  Density  of  aerosol  particle 
Dp  =  Diameter  of  aerosol  particle 
C  =  Cunningham  slip  correction 

The  Cunningham  slip  correction  is  given  by: 


C  -  1  +  -  [6.32  +  2.01-exp(-.1095-P-D„)  (23) 

P*D  *' 

p 

where  P  =  Ambient  pressure  in  cm  Hg 

Dp  =  Diameter  of  aerosol  particle  in  microns 


The  sampling  efficiency  is  given  by: 

2- (W/V)  +0.62 


e  *  1  +  (W/V  -  1) 


(W/V)-Stk'^  +  2 -(W/V)  +0.62 


(24) 


2.5.2  CALCULATION  OF  TRANSPORT  EFFICIENCY  THROUGH  A  STRAIGHT 

TUBE 


In  laminar  flow,  the  aerosol  particle  losses  occur  by 
two  mechanisms.  Brownian  diffusion  and  gravitational 
sedimentation.  The  losses  through  a  section  of  straight  tube  by 
Brownian  diffusion  are  independent  of  the  orientation  of  that 
tube  within  a  gravitational  field.  The  following  are  the 
equations  describing  aerosol  particle  losses  due  to  Brownian 
diffusion  under  laminar  flow  conditions. 
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n  =  0. 819  •  +  0  -  097  •  +  0.032-e'^^^ 

for  X  >  0.02  (25) 

n  =  1  -  2.56-X^^^  +  1.2 -X  +  0.177-x'^''^ 

for  X  <  0.02  (26) 


where  X  =  ir-D-L/Q 

D  =  Brownian  diffusion  coefficient  of  aerosol  particle 
L  =  Length  of  the  tube 
Q  =  Volumetric  flow  rate  within  tube 

The  Brownian  diffusion  coefficient  of  the  aerosol 
particle  is  calculated  from  the  following  equation. 

k-T-C 

D  =  -  (27) 

3  ■  TT  *  Ai  •  Dp 

where  k  =  Boltzman's  constant 

T  =  Absolute  temperature  within  probe 
C  =  Cunningham  slip  correction 
M  =  Dynamic  viscosity  of  fluid 
Dp  =  Diameter  of  aerosol  particle 

The  aerosol  particle  losses  in  a  horizontal  tube  due  to 
gravitational  sedimentation  is  representated  by  the  following 
equation. 

n  =  1  -  -  2‘0-(l  -  4-  sin’^^''^  -  0’''^-  (1  -  0^''^)’''^ 

n  I 

(28) 

Where  0  =  (3/8)  ■  (L/R)  •  (V^/V) 

L  =  Length  of  the  tube 
R  =  Radius  of  the  tube 

=  Sedimentation  velocity 
V  =  Linear  velocity  in  the  tube 

The  aerosol  particle  sedimentation  velocity  can  be  expressed  as: 

=  0.003 -d -Dpi  (29) 

where  d  =  Aerosol  particle  density  in  gm/cm^ 

Dp  =  Aerosol  particle  diameter  in  microns 

Under  conditions  of  turbulent  flow,  the  aerosol 
particle  losses  are  due  to  Brownian  diffusion  and  inertial 
deposition  effects.  The  model  describing  the  transport 


22 


efficiency  is  the  same  for  both  Brownian  diffusion  and  inertial 
deposition,  differing  only  in  the  calculation  of  u,  the  particle 
deposition  velocity.  This  model  is  given  by  the  following 
relation: 


(30) 


where  R  =  Radius  of  the  tube 
L  =  Length  of  the  tube 
Q  =  Volumetric  flow  rate  within  tube 

For  losses  due  to  Brownian  diffusion,  the  particle 
deposition  velocity  is  given  by  the  following  equation. 

u  =  0.042-V'f’^^-Sc'^''^  (31) 

where  f  =  Fanning  friction  factor 
Sc  =  Schmidt  number 
V  =  Linear  velocity  in  the  tube 

The  Schmidt  number  is  given  by: 

V 

Sc  =  -  (32) 

D 

where  v  =  Kinematic  viscosity  of  the  fluid 

D  =  Brownian  diffusion  coefficient  of  aerosol  particle 

The  Fanning  friction  factor  is  given  by; 


f 


(33) 


where  Re  =  Reynold's  number  for  tube 

The  Reynold's  number  for  flow  through  a  tube  is  given 
by  the  following  relation: 

Dg-V-d 

Re  =  -  (34) 

4 

where  =  Diameter  of  the  tube 

V  =  Linear  velocity  in  the  tube 
d  =  Density  of  the  fluid 
4  =  Dynamic  viscosity  of  the  fluid 
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For  inertial  deposition  of  aerosol  particles,  the 
particle  deposition  velocity  is  given  by: 

u  =  u^-u* 

where  u*  =  (f/2)^''^-V 
=  (6  X  lo"^)  • 

=  Dimensionless  relaxation  time 
The  dimensionless  relaxation  time  is  given  by: 


where  r  =  Particle  relaxation  time 
v  =  Kinematic  viscosity 

The  particle  relaxation  time  may  be  calculated  from  the 
following  equation: 

-d-C 


where  Dp  =  Diameter  of  aerosol  particle 
d  =  Density  of  aerosol  particle 
C  =  Cunningham  slip  correction 
M  =  Dynamic  viscosity  of  fluid 


2.5.3 


CALCULATION  OF  TRANSPORT  EFFICIENCY  THROUGH  A  90®  BEND 


For  conditions  of  laminar  flow,  the  transport 
efficiency  of  an  aerosol  particle  through  a  90°  bend  is  given  by; 

n  =  1  -  Stk  (38) 

where  Stk  =  Stokes  number 


For  turbulent  flow,  the  transport  efficiency  of  an 
aerosol  particle  through  a  90°  bend  is  given  by: 


n  =  10 


■0.963Stk 


where  Stk  =  Stokes  number 

For  the  90°  bend  model,  the  Stokes  number  is  given  as: 


C  •  d  •  Dp2  •  V 
9-m-D, 
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where 


V  =  Linear  velocity  within  tube 
=  Inside  diameter  of  tube 
M  =  Dynamic  viscosity  of  fluid 
d  =  Density  of  aerosol  particle 
Dp  =  Diameter  of  aerosol  particle 
C  =  Cunningham  slip  correction 
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SECTION  3  -  AEROSOL  SAMPLING  TRAIN  CONFIGURATION  FOR 
NUMERICAL  CALCULATIONS  AND  EXPERIMENTAL  VERIFICATION 

3.1  PHYSICAL  DESCRIPTION 

The  geometry  of  the  aerosol  sampling  train  which  was 
configured  for  the  model  calculations  and  the  experimental 
verification  is  presented  in  Figure  1.  As  shown  in  the  figure, 
the  sampling  train  is  subdivided  into  six  consecutive  sections. 
They  are  the  following:  inlet,  three  2-foot  long  straight 
sections  -  horizontal,  vertical,  and  inclined  -  connected  by  a 
90°  elbow  and  a  45°  elbow,  respectively.  Both  elbows  have  3-inch 
radii  as  measured  from  their  centers  of  curvature  to  the  elbow 
centerlines.  The  whole  train  is  made  of  a  nominal  3/ 4- inch 
schedule  40  PVC  smooth  wall  electrical  conduit  which  has  an 
inside  diameter  of  1.05  inches.  The  conduit  is  assumed  to  have  a 
surface  roughness  typical  for  drawn  tubing  (brass,  lead,  glass, 
and  the  like)  (Ref.  15).  Two  inlet  configurations  (inlet 
oriented  parallel  to  the  air  stream  (0°)  and  perpendicular  to  the 
air  stream  (90°))  are  being  considered. 

3.2  REGION  OF  AEROSOL  SAMPLING  TRAIN  ADDRESSED  BY  FOUR 

MODELS 

Each  of  the  four  candidate  "models"  discussed  in 
Section  2  simulate/characterize  different  regions  of  the  six 
section  aerosol  sampling  train.  Each  model's  characterization  of 
the  different  regions  of  the  sampling  train  is  presented  in  Table 
1.  As  Table  1  points  out,  the  Texas  A&M  and  Battelle  models  are 
the  only  two  which  simulate  the  90°  and  45°  bends,  and  the 
horizontal,  vertical  and  inclined  sections.  The  Texas  ASeM  model 
is  the  only  one  which  simulates  the  six  regions  of  the  sampling 
train.  Figure  2  presents  graphically  the  regions  of  the  sampling 
train  which  are  characterized/simulated  in  each  model. 

3.3  PARAMETERS  OF  MODEL  CALCULATIONS  AND  EXPERIMENTAL 

VERIFICATION 

Model  calculations  and  experimental  data  (to  determine 
cumulative  sampling  efficiency  (%  mass))  were  generated  for 
sampling  flow  rates,  70  1/min  and  130  1/min,  and  an  ambient  wind 
speed  of  3  m/sec.  Seven  monodispersed  aerosols  having  aerodynamic 
equivalent  particle  diameter  (AED)  of  3,5,7,10,15,20,  and  25 
micrometers  are  generated  and  sampled.  AED  is  defined  as  the 
diameter  of  a  water  droplet  of  unit  specific  gravity  which  has 
the  same  time  constant,  t,  as  an  arbitrary  particle  of  density  d 
and  equivalent  diameter  Dp.  In  addition,  two  log-normal  aerosol 
distributions  were  generated  having  5  and  10  aerodynamic  mass 
median  diameter  (micrometers)  and  each  with  a  geometric  standard 
deviation  of  1.5  micrometers.  Unit  density  was  assumed  for  all 
airborne  aerosol  particles.  A  matrix  containing  the  36 
simulation  cases  to  be  conducted  is  presented  in  Table  2. 
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Figure  1.  Schematic  of  Aerosoi  Sampiing  Train 


REGIONS  OF  AEROSOL  SAMPLING  TRAIN  SIMULATED/CHARACTERIZED 


Table  1.  Four  Candidates  Models’  CharacterIzatlon/SImulation  of  Aerosol  Sampling  Train 
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Figure  2.  Regions  of  Sampling  Train  Characterized/Simuiated  by  Each  Candidate  Model 


TYPE  OF  PAHTICLE  INLET  PROBE  FLOWRATE  WIND  SPEED 

PARTICLE  DISTRIBUTION  SIZE  (  M  1  ORIgNTATION  ilZmiDl  (m/8W> 


1.  MONOOiSPERSEO 

3 

0° 

70 

3 

2.  MONOOISPERSEO 

5 

0= 

70 

3 

3.  MONOOISPERSEO 

7 

0° 

70 

3 

4.  MONOOISPERSEO 

10 

0° 

70 

3 

5.  MONOOISPERSEO 

15 

0° 

70 

3 

6.  MONOOISPERSEO 

20 

0° 

70 

3 

7.  MONOOISPERSEO 

25 

0° 

70 

3 

8.  MONOOISPERSEO 

3 

0° 

130 

3 

9.  MONOOiSPERSEO 

5 

0° 

130 

3 

10.  MONOOISPERSEO 

7 

0° 

130 

3 

11.  MONOOISPERSEO 

10 

o“ 

130 

3 

12.  MONOOISPERSEO 

15 

0° 

130 

3 

13.  MONOOISPERSEO 

20 

0° 

130 

3 

14.  MONOOISPERSEO 

25 

0° 

130 

3 

15.  MONOOISPERSEO 

3 

90° 

70 

3 

16.  MONOOISPERSEO 

5 

90° 

70 

3 

17.  MONOOISPERSEO 

7 

90° 

70 

3 

18.  MONOOISPERSEO 

10 

90° 

70 

3 

19.  MONOOISPERSEO 

IS 

90° 

70 

3 

20.  MONOOISPERSEO 

20 

90° 

70 

3 

21.  MONOOISPERSEO 

25 

90° 

70 

3 

22.  MONOOISPERSEO 

3 

90° 

130 

3 

23.  MONOOiSPERSEO 

5 

90° 

130 

3 

24.  MONOOISPERSEO 

7 

90° 

130 

3 

10 

90° 

130 

3 

26.  MONOOISPERSEO 

IS 

90° 

130 

3 

27.  MONOOISPERSEO 

20 

90° 

130 

3 

28.  MONOOISPERSEO 

25 

90° 

130 

3 

29.  POLYDISPERSEO 

MMA0«S<  0r«1.5) 

0° 

70 

3 

30.  POLYDISPERSEO 

MMA0«15(C7»1.5) 

0° 

70 

3 

31.  POLYDISPERSEO 

MMA0«5(  a-1-5) 

130 

3 

3Z  POLYDISPERSEO 

MMA0=>t5(a»1-5) 

0° 

130 

3 

MMAO«5<or=1S> 

90° 

70 

3 

34.  POLYDISPERSEO 

MMA0*1S(O'«1S) 

90° 

70 

3 

35.  POLYDISPERSEO 

MMAD«S  (  a -1-5) 

90° 

130 

3 

36.  POLYDISPERSEO 

MMAD-rtS(Cr>1.5| 

90° 

130 

3 

Table  2.  Matrix  of  Simulation  Cases  to  be  Generated  for  Aerosol 

Sampling  Train  Configuration 
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SECTION  4  -  RESULTS  OP  MODEL  NUMERICAL  CALCULATIONS 


4.1  RESULTS  FROM  TEXAS  A&M  MODEL 

As  mentioned  earlier,  numerical  calculations  using  the 
Texas  A&M  model  were  made  for  two  log-normally  distributed 
aerosols  and  seven  monodispersed  aerosols  for  two  flow  rates  and 
two  probe  inlet  orientations-  The  cumulative  sampling  efficiency 
results  (%  mass)  for  these  36  simulation  cases  are  summarized  in 
Table  3  and  graphically  presented  in  Figures  3  and  4 .  The 
complete  set  of  numerical  calculations  produced  by  the  Texas  A&M 
model  are  provided  at  Appendix  B. 

4.1.1  Results  of  Monodispersed  Simulation  Cases 

The  cumulative  sampling  efficiency  results  for  a 
monodispersed  aerosol  distribution  indicate  the  following: 

•  As  the  particle  size  increases,  the  sampling  efficiency 
decreases. 

•  An  increase  in  the  inlet  probe  orientation  (0°  to  90°) 
will  result  in  a  decreased  sampling  efficiency. 

•  As  the  flow  rate  increases  from  70  1/min  to  130  1/roin, 
the  sampling  efficiency  decreases. 

•  The  larger  the  particle  size  the  greater  the  decrease 
in  sampling  efficiency. 

•  Increase  in  the  inlet  probe  orientation  produces  a  more 
dramatic  decrease  in  the  sampling  efficiency  than  does 
an  increase  in  the  flow  rate. 

•  Particle  depositional  losses  in  the  elbows  are 
generally  larger  than  those  in  the  other  components  of 
the  sampling  train. 

4.1.2  Results  of  Polydispersed  Simulation  Cases 

The  cumulative  sampling  efficiency  results  for  a 
polydispersed  aerosol  distribution  indicate  the  following: 

•  As  the  mass  median  aerodynamic  diameter  (MMAD) 
increases,  the  sampling  efficiency  decreases. 

•  As  the  flow  rate  increases,  the  sampling  efficiency 
decreases. 
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TYPE  OF 


PARTICLE 


INLET  PROBE 


TYPE  OF 

PARTICLE  DISTRIBUTION 

PARTICLE 

glZE<,.iU 

INLET  PROBE 
ORIENTATION  t  Q  ) 

FLOW  RATE 
(l/mln> 

WIND  SPEED 
InVec) 

CUMULATIVE 

SAMPLING 

EFFICIENCY 

1.  MONODISPERSED 

3 

0° 

70 

3 

97.4 

2.  MONODISPERSED 

5 

0° 

70 

3 

93.6 

3.  MONODISPERSED 

7 

0° 

70 

3 

88.4 

4.  MONODISPERSED 

10 

0° 

70 

3 

78.6 

5.  MONODISPERSED 

15 

0° 

70 

3 

59.0 

6.  MONODISPERSED 

20 

0° 

70 

3 

39.2 

7.  MONODISPERSED 

25 

0° 

70 

3 

22.8 

8.  MONODISPERSED 

3 

0° 

130 

3 

96J 

9.  MONODISPERSED 

5 

0° 

130 

3 

90.4 

10.  MONODISPERSED 

7 

0° 

130 

3 

82.5 

11.  MONODISPERSED 

10 

0° 

130 

3 

67.8 

12.  MONODISPERSED 

15 

0° 

130 

3 

41.0 

13.  MONODISPERSED 

20 

0° 

130 

3 

16.6 

14.  MONODISPERSED 

25 

0° 

130 

3 

5.6 

15.  MONODISPERSED 

3 

90° 

70 

3 

94.3 

16.  MONODISPERSED 

5 

90° 

70 

3 

85.6 

17.  MONODISPERSED 

7 

90° 

70 

3 

74.7 

18.  MONODISPERSED 

10 

90° 

70 

3 

57.0 

19.  MONODISPERSED 

15 

90^ 

70 

3 

31.5 

20.  MONODISPERSED 

20 

90° 

70 

3 

15.2 

21.  MONODISPERSED 

25 

90° 

70 

3 

6.4 

22.  MONODISPERSED 

3 

90° 

130 

3 

94.3 

23.  MONODISPERSED 

5 

90° 

130 

3 

85.5 

24.  MONODISPERSED 

7 

90° 

130 

3 

74.0 

25.  MONODISPERSED 

10 

90° 

130 

3 

55.1 

26.  MONODISPERSED 

15 

90° 

130 

3 

27.1 

27.  MONODISPERSED 

20 

90° 

130 

3 

8.7 

28.  MONODISPERSED 

25 

90° 

130 

3 

Z4 

29.  POLYOISPERSEO 

MMAO«5(  (T>1.5| 

0° 

70 

3 

91.9 

30.  POLYDISPERSED 

MMAD3l5(0'=1.5) 

0° 

70 

3 

55.5 

31.  POLYOISPERSEO 

MMAD-5  (  0-  -1.5) 

0° 

130 

3 

87.8 

3Z  POLYOISPERSEO 

MMAO»15(a-1.5) 

0° 

130 

3 

40.1 

33.  POLYDISPERSED 

MMAO-S((7-1.S) 

90° 

70 

3 

82.6 

34.  POLYOISPERSEO 

MMAO«15((7«1.S) 

90° 

70 

3 

33.4 

35.  POLYOISPERSEO 

MMAO«5(Or  «1.5) 

90° 

130 

3 

82J2 

36.  POLYOISPERSEO 

MMAO«15(a«1.5) 

90° 

130 

3 

29JI 

Table  3.  Cumulative  Sampling  Efficiency  Results  by  Texas  A  &  M  Model 
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2  6  10  U  18  22  26 

Moss  Median  Aerodynamic  Diameter  (^) 
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•  As  the  inlet  probe  orientation  increases,  the  sampling 
efficiency  decreases. 

•  Increase  in  the  inlet  probe  orientation  produces  a  more 
dramatic  decrease  in  the  sampling  efficie"::y  than  does 
an  increase  in  the  flow  rate. 

•  The  5/Lini  MMAD  polydispersed  particle  distribution 
produces  smaller  sampling  efficiencies  than  does  the 
Sum  monodispersed  particle  distribution.  However,  the 
differences  in  efficiency  (typically  a  few  percent)  are 
considered  insignificant. 

•  The  15/im  MMAD  polydispersed  particle  distribution 
produces  smaller  sampling  efficiencies  than  does  the  15 
fim  monodispersed  particle  distribution  for  only  the  0° 
inlet  probe  orientation  and  70  1/min  flow  rate.  Again, 
these  differences  are  considered  insignificant.  For 
other  combinations  of  inlet  probe  orientation  and  flow 
rate,  the  15  /xm  MMAD  polydispersed  particle 
distribution  produces  greater  sampling  efficiencies 
than  does  the  15  /xm  monodispersed  particle 
distribution.  As  with  the  other  comparisons,  the 
efficiency  differences  are  considered  insignificant. 

4.2  RESULTS  FROM  IITRI  MODEL 

As  discussed  earlier,  the  IITRI  model  calculates 
sampling  efficiencies  based  on  characterization  of  the  inlet 
probe  only.  None  of  the  other  regions  of  the  aerosol  sampling 
train  are  addressed.  The  complete  set  of  numerical  calculations 
produced  by  the  IITRI  model  is  provided  at  Appendix  C. 

Sampling  efficiencies  as  a  function  of  monodispersed 
particle  size  were  calculated  for  28  simulation  cases  and  are 
presented  in  Table  4  and  graphically  presented  in  Figure  5. 

These  results  indicate  the  following: 

•  As  the  particle  size  increases,  the  sampling  efficiency 
decreases  for  the  conditions  of  an  inlet  probe 
orientation  of  0°  and  90°  and  a  flow  rate  of  130  1/min, 
and  also  the  cases  for  an  inlet  probe  orientation  of 
90°  and  a  flow  rate  of  70  1/min.  The  cases  for  a  0° 
inlet  probe  orientation  and  70  1/min  flow  rate  resulted 
in  an  increase  in  the  sampling  efficiency  as  the 
particle  size  increases. 

•  The  following  conditions  in  descending  order  resulted 
in  the  greatest  decrease  in  sampling  efficiency  as  the 
particle  size  increases. 


TYPE  OF 

PARTICLE  DISTRIBUTION 

PARTICLE 

$IZ£,LA-1 

INLET  PROBE 
ORIENTATION  1  a  ) 

FLOW  RATE 
(l/tTtlnl 

WIND  SPEED 
(m/»ecl 

SAMPLING 

EFFICIENCY 

OF  INLET  PROBE 
1%  maM) 

1.  MONODiSPERSED 

3 

0° 

70 

3 

100.7 

2.  MONODISPERSED 

5 

0° 

70 

3 

100.7 

3.  MONODISPERSED 

7 

0° 

70 

3 

100.7 

4.  MONODISPERSED 

10 

0^^ 

70 

3 

100.7 

5.  MONODISPERSED 

15 

0° 

70 

3 

108.4 

6.  MONODISPERSED 

20 

0° 

70 

3 

114.7 

7.  MONODiSPERSED 

25 

0° 

70 

3 

121.2 

8.  MONODISPERSED 

3 

0° 

130 

3 

94.1 

9.  MONODISPERSED 

5 

0° 

130 

3 

94.1 

10.  MONODISPERSED 

7 

0° 

130 

3 

94.1 

11.  MONODISPERSED 

10 

0° 

130 

3 

93.0 

12.  MONODISPERSED 

15 

0° 

130 

3 

91.9 

13.  MONODISPERSED 

20 

0° 

130 

3 

91.9 

14.  MONODiSPERSED 

25 

0° 

13ti 

3 

90.9 

15.  MONODISPERSED 

3 

90° 

70 

3 

89.3 

16.  MONODiSPERSED 

5 

90° 

70 

3 

89.3 

17.  MONODISPERSED 

7 

90° 

70 

3 

88.3 

18.  MONODISPERSED 

10 

90° 

70 

3 

87.2 

19.  MONODISPERSED 

15 

90° 

70 

3 

79.9 

20.  MONODISPERSED 

20 

90° 

70 

3 

75.9 

21.  MONODISPERSED 

25 

90° 

70 

3 

71.3 

22.  MONODISPERSED 

3 

90° 

130 

3 

85.4 

23.  MONODISPERSED 

5 

90° 

130 

3 

85.3 

24.  MONODISPERSED 

7 

90° 

130 

3 

84.7 

25.  MONODiSPERSED 

10 

90° 

130 

3 

84.8 

26.  MONODISPERSED 

15 

90° 

130 

3 

aiJt 

27.  MONODiSPERSED 

20 

90° 

130 

3 

76.1 

28.  MONODISPERSED 

25 

90° 

130 

3 

70.2 

Table  4.  Cumulative  Sampling  Efficiency  Results  Generated  by 

IITRI  Model 
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Figure  5.  Penetration  Efficiency 


IITRI  Model  Monodispersed  Results 


AerodyTKimic  Particle  Diometer  (/jm) 

70  L/min  O’  +130  L/min  0*  P  70  L/min  90*  A  130  L/min  90 
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90°  inlet  probe  orientation  and  130  1/min 
flow  rate 

90°  inlet  probe  orientation  and  70  l/min  flow 
rate 

0°  inlet  probe  orientation  and  130  1/itiin  flow 
rate 

4.3  RESULTS  FROM  UNIVERSITY  OF  DUISBURG  MODEL 

4.3.1  Background 

As  mentioned  earlier,  the  Fissan  and  Schwientek  model 
simulates  only  one  inlet  probe  orientation  (0°)  and  does  not 
account  for  a  45°  bend  or  inclined  section  of  the  sampling  train- 
Furthermore,  it  does  not  predict  cumulative  sampling  efficiency 
resulting  from  a  log-normal  particle  distribution.  GMA,  as  part 
of  its  effort  to  develop  a  computer  program  for  the  German  model, 
extended  the  model's  capability  to  include  a  45°  bend,  an 
inclined  region  of  the  train,  and  log-normal  distribution. 

This  section  of  the  report  is  to  focus  primarily  on  the 
cumulative  sampling  efficiencies  generated  by  the  model. 
Therefore,  the  results  generated  by  the  model  and  presented  in 
this  section  will  not  include  a  45°  bend,  inclined  region,  or 
log-normal  distribution.  However,  Section  6  will  include  a 
comparison  of  the  model  results  to  the  experimental  data  gathered 
by  Texas  A&M.  The  comparison  will  include  the  model  results  for 
the  45°  bend,  inclined  section,  and  log-normal  distribution.  THe 
complete  set  of  numerical  calculations  prcduced  by  the  GMA 
computer  program  is  provided  at  Appendix  b. 

4.3.2  Results  of  Monodispersed  Simulation  Cases 

The  cumulative  sampling  efficiency  results  presented  in 
Table  5  and  Figure  6  for  a  monodispersed  aerosol  distribution 
indicate  the  following: 

•  As  the  particle  size  increases  the  sampling  efficiency 
decreases . 

•  As  the  flow  rate  increases  from  70  l/min  to  130  l/min, 
the  sampling  efficiency  decreases. 

•  The  larger  the  particle  size,  the  greater  the  decrease 
in  cumulative  sampling  efficiency. 

Cumulative  sampling  efficiency  resluts  generated  by  the  GMA 
computer  program  for  a  polydispersed  aerosol  distribution  has 
been  presented  in  Figure  7  for  illustrative  purposes  only. 
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TYPE  OF 

PAnTICUE  DISTRIBUTlOri 

1.  MONOOISPERSEO 

PARTICLE 

SIZE  (  M  ) 

3 

INLET  PROBE 
ORieNTATlQN.t  ) 

0° 

FLOW  RATE 
(l/mln> 

70 

WIND  SPEED 
imiifisi 

3 

CUMULATIVE 

SAMPLING 

EFFICIENCY* 

98.4 

2.  MONOOISPERSEO 

5 

0° 

70 

3 

95.7 

3.  MONOOISPERSEO 

7 

0° 

70 

3 

91.7 

4.  MONOOISPERSEO 

10 

0° 

70 

3 

83.3 

5.  MONOOISPERSEO 

15 

0° 

70 

3 

64.8 

6.  MONOOISPERSEO 

20 

0° 

70 

3 

43.8 

7.  MONOOISPERSEO 

25 

0° 

70 

3 

25.6 

8.  MONOOISPERSEO 

3 

0° 

130 

3 

95.9 

9.  MONOOISPERSEO 

5 

0° 

130 

3 

69.1 

10.  MONOOISPERSEO 

7 

0° 

130 

3 

79.8 

11.  MONOOISPERSEO 

10 

0° 

130 

3 

63.0 

12.  MONOOISPERSEO 

15 

0° 

130 

3 

34.6 

13.  MONOOISPERSEO 

20 

0° 

130 

3 

14.2 

14.  MONOOISPERSEO 

25 

130 

3 

3.9 

15.  POLYOISPERSEO 

MMA0=5(  a  =1.5) 

0° 

70 

3 

94.1 

16.  POLYOISPERSEO 

MMA0=15(CT=1.5) 

0° 

70 

3 

59.9 

17.  POLYOISPERSEO 

MMA0=5(  (7  =1.5) 

0° 

130 

3 

86.0 

18..  POLYOISPERSEO 

MMAO=15(a=1.5) 

0° 

130 

3 

35.8 

*  Cumulativ*  Etflclancy  Reaults  raflact  modlflcationa  mada  to  tha  modal  by  Ganaral  Managameni 
Aaaociataa  to  allow  calculations  for  45°  albow  and  inclinad  aactlon  of  the  aampling  train. 


Table  5.  Cumulative  Sampling  Efficiency  Results  Generated  by 

University  of  Duisburg  Model 
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Percent  Penetration 


Figure  6.  Penetration  Efficiency 


Duisburg  Model  Monodisoersed  Resulte 


Aerodynamic  Particle  Diameter 
O  70  L/min  O’  +  1 30  L/min  ff 
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f'ercont  Penetration 


Figure  7,  Penetration  Efficiency 


Duisburg  Model  Poiydispersed  Results 


2  6  10  14  18  22  26 

Moss  Median  Aerodynamic  Diameter  (^) 

0  70  L/min  +130  L/min  CT 
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4.4  RESULTS  FROM  BATTELLE  MODEL 

For  background  purposes,  the  Battelle  model  does  not 
account  for  an  inlet  probe  orientation.  However,  the  model  does 
account  for  the  three  straight  sections  of  the  train  and  the  90° 
bend  and  45°  bend.  Therefore,  the  model  calculates  cumulative 
sampling  efficiencies  for  18  simulation  cases  which  are 
summarized  in  Table  6  and  graphically  presented  in  Figures  8  and 
9.  The  complete  set  of  numerical  calculations  produced  by  the 
Battelle  model  are  provided  at  Appendix  E. 

4.4.1  Results  of  Monodispersed  Simulation  cases 

The  cumulative  sampling  efficiency  results  for  a 
monodispersed  aerosol  distribution  indicate  the  following; 

•  As  the  particle  size  increases,  the  sampling  efficiency 
decreases . 

•  An  increase  in  the  flow  rate  results  in  a  decrease  in 
the  sampling  efficiency.  However,  the  differences  are 
small  (less  than  2  percent) . 

•  The  larger  the  particle  size,  the  greater  the  decrease 
in  sampling  efficiency. 

•  Particle  depositional  losses  in  the  elbows  are 
generally  larger  than  those  in  the  other  components  of 
the  sampling  train. 

4.4.2  Results  of  Polydispersed  Simulation  cases 

The  cumulative  sampling  efficiency  results  for  a 
polydispersed  aerosol  distribution  indicate  the  following; 

•  As  the  MMAD  particle  size  increases,  the  sampling 
efficiency  decreases. 

•  As  the  flow  rate  increases,  the  sampling  efficiency 
decreases.  However,  the  decrease  is  not  considered 
significant  for  the  range  of  particles  evaluated. 

•  Particle  depositional  losses  in  the  elbows  are 
generally  larger  than  those  in  the  other  compartments 
of  the  sampling  train. 

•  The  5/im  MMAD  polydispersed  particle  distribution 
produces  smaller  sampling  efficiencies  than  does  the 
5^^m  monodispersed  particle  distribution.  However,  the 
differences  in  efficiency  (typically  less  than  one 
percent)  are  considered  insignificant. 
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TYPE  OF 

PARTICLE  DISTRIBUTION 

PARTICLE 

INLET  PROBE 
ORIENTATION  i  O  ) 

FLOW  RATE 
<l/mln) 

WIND  SPEED 
(m/iftcj 

CUMULATIVE 
SAMPLING 
EFFICIENCY* 
(%  mass! 

1.  MONODiSPERSED 

3 

N/A 

7C 

3 

99.0 

2.  MONODISPERSED 

5 

N/A 

70 

3 

97.8 

3.  MONODISPERSED 

7 

N/A 

70 

3 

96.0 

4.  MONODISPERSED 

10 

N/A 

70 

3 

92.4 

5.  MONODISPERSED 

15 

N/A 

70 

3 

83.9 

6.  MONODISPERSED 

20 

N/A 

70 

3 

73.1 

7.  MONODISPERSED 

25 

N/A 

70 

3 

60.7 

8.  MONODISPERSED 

3 

N/A 

130 

3 

98.8 

8.  MONODISPERSED 

5 

N/A 

130 

3 

98.0 

to.  MONODISPERSED 

7 

N/A 

130 

3 

96.7 

t1.  MONODISPERSED 

10 

N/A 

130 

3 

93.9 

12.  MONODISPERSED 

15 

N/A 

130 

3 

83.7 

13.  MONODISPERSED 

20 

N/A 

130 

3 

71.7 

14.  MONODISPERSED 

25 

N/A 

130 

3 

60.2 

15.  POLYDISPERSED 

MMAD=:S(  0=1.5) 

N/A 

70 

3 

97.1 

16.  POLYDISPERSED 

MMAD=15(  0=1.5) 

N/A 

70 

3 

79.8 

17.  POLYDISPERSED 

MMA0»5(  0  =1.5) 

N/A 

130 

3 

97.4 

18..  POLYDISPERSED 

MMAD=15(  0=1.5) 

N/A 

130 

3 

80.7 

*  Efflcltney  from  difforont  probo  oriontation*  is  not  consldarsd  in  the  Battelie  model. 


Table  6.  Cumulative  Sampling  Efficiency  Results  Generated  by 

Battelie  Model 
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Percent  Penetration 


Figure  8.  Penetration  Efficiency 


Bottelle  Model  Monodispersed  Results 


Aerodynamic  Particle  Oiometer  ^jm) 
□  70  L/min  +  130  L/min 
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•  The  I5fxm  MMAD  polydispersed  particle  distribution 
produces  smaller  sampling  efficiencies  than  does  the 
15Min  monodispersed  particle  distribution.  These 
differences  (3  to  4  percent)  are  larger  than  those  seen 
with  the  S/iiti  particle  size  comparisons. 

4.5  COMPARISON  OF  NUMERICAL  CALCULATIONS 

In  order  to  perform  a  comparison  of  the  numerical 
results  among  the  four  models,  it  is  necessary  to  establish  a 
common  basis  upon  which  the  results  can  be  compared.  The  four 
models  individually  do  not  provide  numerical  calculations  for  the 
six  sections  of  the  sampling  train.  Therefore,  the  common  basis 
for  comparison  will  be  the  numerical  results  for  specific  regions 
of  the  sampling  train  characterized/simulated  by  each  model. 

4.5.1  Numerical  Model  Calculations  for  Inlet  Region  of 
Aerosol  Sampling  Train 

4.5. 1.1  Monodispersed  Particle  Distribution  -  inlet  Region 

Table  7  presents  the  three  model's  numerical  results 
for  the  inlet  region  of  the  sampling  train.  Figures  10  and  11 
provide  a  graphical  presentation  of  the  numerical  results.  It  is 
important  to  note  that  only  the  Texas  A&M  model  and  the  IITRI 
model  consider  air  stream  orientations  other  than  parallel  to  the 
inlet.  The  University  of  Duisburg  model  only  considers  a  0° 
inlet  orientation. 

The  inlet  sampling  efficiency  results  for  a 
monodispersed  particle  distribution  indicate  the  following: 

•  Both  the  Texas  A&M  and  IITRI  models  predict  a  decrease 

as  the  inlet  probe  orientation 
to  90®. 

•  All  three  models  predict  a  decrease  in  the  inlet 
efficiency  as  the  flow  rate  increases. 

•  The  Texas  A&M  model  predicts  a  decrease  in  the  inlet 
efficiency  as  the  particle  size  increases  for  all  four 
combinations  of  inlet  orientation  and  flow  rate.  The 
IITRI  and  German  models,  on  the  other  hand,  predict  an 
increase  in  the  efficiency  as  the  particle  size 
increases  for  the  cases  of  0®  inlet  orientation  and  70 
1/min  flow  rate-  However,  for  the  other  combinations 
of  inlet  orientation  and  flow  rate,  these  two  models 
predict  a  decrease  in  the  efficiency  as  the  particle 
size  increases. 


in  the  efficiency 
increases  from  0° 
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PARTICLE 

SIZEf  ^  ) 


INLET  PROBE 
ORIENTATION  (  O  ) 


FLOW  RATE 
fl/mlnl 


SAMPLING  EFFICIENCY  (%  massl  FOR  INLET  REGION 


TEXASAAM 


Him 


UNIVERSITY  OF 

ByiSBUflg 


1. 

3 

0° 

70 

99.4 

100.7 

100.6 

z 

5 

0° 

70 

98.9 

100.7 

101.8 

3. 

7 

0° 

70 

98.4 

100.7 

103.3 

4. 

10 

0° 

70 

97.7 

100.7 

106.2 

5. 

15 

70 

96.4 

108.4 

111.9 

6. 

20 

0° 

70 

94.7 

114.7 

117.4 

7. 

25 

0° 

70 

92.4 

121.2 

122.2 

8. 

3 

0° 

130 

99.3 

94.1 

99.6 

9. 

5 

0° 

130 

98.3 

94.1 

98.9 

10. 

7 

0° 

130 

97.1 

94.1 

98.0 

11. 

10 

0° 

130 

94.9 

93.0 

96.3 

12. 

15 

0° 

130 

90.6 

91.9 

93.2 

13. 

20 

130 

85.7 

91.9 

90.2 

14. 

25 

0° 

130 

80.5 

90.9 

87.7 

15. 

3 

90° 

70 

96.2 

89.3 

N/A 

16. 

5 

90° 

70 

90.5 

89.3 

N/A 

17. 

7 

90° 

70 

83JZ 

88J 

N/A 

18. 

10 

90° 

70 

70.9 

87.2 

N/A 

19. 

15 

90° 

70 

51 

79.9 

N/A 

20. 

20 

90° 

70 

36.6 

75.9 

N/A 

21. 

25 

90° 

70 

26.1 

71.3 

N/A 

22. 

3 

90° 

130 

97.2 

85.4 

N/A 

23. 

5 

90° 

130 

92.9 

85.3 

N/A 

24. 

7 

90° 

130 

87.2 

84.7 

N/A 

25. 

10 

90° 

130 

T7JZ 

84.8 

N/A 

26. 

15 

90° 

130 

59.8 

81.2 

N/A 

27. 

20 

90° 

130 

45.1 

76.1 

N/A 

28. 

25 

90° 

130 

33.8 

70.2 

N/A 

Table  7.  Numerical  Model  Calculations  (Monodlspersed  Particle  Sizes) 
for  Inlet  Region  of  Aerosol  Sampling  Train 
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Percent  Penetration 


Figure  11.  inlet  Penetration  Efficiency 


Monodispersad  Resuite  -  9(7  Inlet 


Aerodynamic  Particle  Diameter  (/jm) 

0  70  L/min  Texos  AisM  +  130  l/min  Texas  A&M  O  70  L/min  IfTRI 

t  130  L/min  IITRI 
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4. 5. 1.2  Polydispersed  Particle  Distribution  -  Inlet  Region 

Table  8  presents  two  model's  numerical  efficiency 
results  for  the  inlet  region  of  the  sampling  train.  Figures  12 
and  13  provide  a  graphical  presentation  of  the  numerical  results. 
It  is  important  to  note  that  only  the  Texas  A&M  model  considers 
air  stream  orientations  other  than  parallel  to  the  inlet.  The 
IITRI  model  does  not  consider  polydispersed  particle 
distribution.  The  University  of  Duisburg  model  as  extended  by 
GMA  considers  polydispersed  particle  distribution  but  only  for  a 
0°  inlet  orientation.  The  inlet  sampling  efficiency  results  for 
polydispersed  particle  distribution  indicate  the  following: 

•  For  the  Texas  A&M  model,  as  the  particle  MMAD 
increases,  the  sampling  efficiency  decreases.  The 
University  of  Duisburg  model,  on  the  other  hand, 
predicts  an  increase  or  a  decrease  in  the  sampling 
efficiency  depending  on  the  flow  rate. 

•  Both  models  predict  a  decrease  in  the  sampling 
efficiency  as  the  flow  rate  increases  when  the  inlet 
orientation  is  parallel  to  the  air  stream.  The  Texas 
A&M  results  indicate  an  increase  in  the  sampling 
efficiency  as  the  flow  rate  increases  when  the  inlet 
orientation  is  normal  to  the  air  stream. 

4.5.2  Nxunerical  Efficiency  Calculations  for  Horizontal  Region 
of  Aerosol  Sampling  Train 

4. 5. 2.1  Monodispersed  Particle  Distribution  -  Horizontal  Region 

Table  9  presents  the  three  model's  numerical  results 
for  the  horizontal  region  of  the  sampling  train.  Figure  14 
provides  a  graphical  presentation  of  the  numerical  results.  It 
is  important  to  note  that  for  these  calculations,  the  horizontal 
region  is  not  affected  by  the  inlet's  orientation  to  the  air 
stream.  Therefore,  only  the  flow  rate  and  particle  size  have 
been  varied. 

The  horizontal  region's  sampling  efficiency  results  for 
a  monodispersed  particle  distribution  indicate  the  following: 

•  All  three  models  predict  a  decrease  in  the  efficiency 
as  the  particle  size  increases. 

•  Both  the  Texas  A&M  and  the  Battelle  models  predict  an 
increase  in  sampling  efficiency  as  the  flow  rate 
increases.  In  addition  the  Texas  A&M  results  indicate 
a  decrease  in  the  efficiency  at  the  20  and  25|um 
particle  sizes  as  the  flow  rate  increases.  The  German 
results,  on  the  other  hand,  indicate  a  decrease  in  the 
efficiency  as  the  flow  rates  increase. 
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MASS  MEDIUM 
AERODYNAMIC 
DIAMETER 
(MMADH  M  ) 

INLET  PROBE 

0ni£NTATlQHC.Q  1 

FLOW  RATE 
|l/mln> 

SAMPLING  EFFICIENCY  FOR  INLET  REGION 

TEXASAAM  illBl  BATTELLE  UNIVERSITY  OF 

P-UiSBURQ 

1.  5 

0° 

70 

98.8 

N/A 

N/A 

102.3 

2.  15 

0° 

70 

95.8 

N/A 

N/A 

112.9 

3.  5 

0° 

130 

98.0 

N/A 

N/A 

98.6 

4.  15 

0° 

130 

89.1 

N/A 

N/A 

9^7 

5.  5 

90° 

70 

88.4 

N/A 

N/A 

N/A 

6.  15 

90° 

70 

51.0 

N/A 

N/A 

N/A 

7.  5 

90° 

130 

91.2 

N/A 

N/A 

N/A 

8.  15 

90° 

130 

58.4 

N/A 

N/A 

N/A 

Table  8.  Model  Efficiency  Calculations  (Polydispersed  Particle  Sizes) 
for  Inlet  Region  of  Aerosol  Sampling  Train 
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PARTICLE 

FLOW  RATE 

SAMPLING  EFFICIENCY  f%  ma»«>  FOR  HORIZONTAL  REGION 

SIZE  1  M  1 

(l/mlnl 

TEXAS  A  A  M 

IITRI 

BATTELLE 

UNIVERSITY  OF 

cuisewRo 

1. 

3 

70 

99.7 

N/A 

99.6 

100.0 

2. 

5 

70 

99.3 

N/A 

98.9 

99.99 

3. 

7 

70 

98.6 

N/A 

97.9 

99.98 

4. 

10 

70 

97i 

N/A 

95.8 

99.9 

5. 

15 

70 

93.8 

N/A 

90.7 

99.6 

6. 

20 

70 

89.1 

N/A 

84.1 

98.8 

7. 

25 

70 

83.2 

N/A 

76.3 

97.2 

8. 

3 

130 

99.9 

N/A 

99.8 

99.99 

9. 

5 

130 

99.6 

N/A 

99.4 

99.96 

10. 

7 

130 

99.2 

N/A 

98.9 

99.9 

11. 

10 

130 

98.3 

N/A 

97.7 

99.4 

15 

130 

95.6 

N/A 

94.9 

17.1 

13. 

20 

130 

88.0 

N/A 

91.1 

91.0 

14. 

25 

130 

80.2 

N/A 

86.4 

79.8 

Table  9.  Model  Efficiency  Calculations  (Monodlspersed  Particle  Sizes) 
for  Horizontai  Region  of  Aerosoi  Sampiing  Train 
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Figure  14.  Penetration  Efficiency 


Monodispensed  Results  —  Horizcntal  Sec. 


Aenodynamic  Particle  Diameter  (^) 

D  70  L/min  Texas  A&M  +  130  L/min  Texas  A&M  0  70  L/min  Bottelle 

A  1 30  L/min  Bottelle  X  70  l/min  Duisburg  7  1 30  L/min  Duisburg 


4. 5. 2. 2  Polydispersed  Particle  Distribution  -  Horizontal  Region 

Table  10  presents  the  three  model's  numerical  results 
for  the  horizontal  region  of  the  sampling  train.  Figures  15  and 
16  provide  a  graphical  presentation  of  the  numerical  results.  It 
is  important  to  note  that  only  the  Texas  A&M  model  considers  air 
stream  orientations  other  than  parallel  to  the  inlet.  The  German 
model  only  considers  a  0°  inlet  orientation.  The  Battelle  model 
does  not  consider  inlet  orientations. 


The  horizontal  sampling  efficiency  results  for  a 
polydispersed  particle  distribution  indicate  the  following: 

•  All  three  models  predict  a  decrease  in  the  sampling 
efficiency  as  the  particle  MMAD  increases. 

•  The  Texas  A&M  and  Battelle  models  predict  an  increase 
in  the  sampling  efficiency  as  the  flow  rate  increases 
for  a  0°  inlet  orientation.  The  German  model  predicts 
a  decrease  in  the  sampling  efficiency  as  the  flow  rate 
increases.  The  Texas  A&M  model  also  predicts  an 
increase  in  the  sampling  efficiency  as  the  flow  rate 
increases  for  a  90°  inlet  orientation. 


•  All  three  models  predict  a  decrease  in  the  sampling 
efficiency  as  the  flow  rate  increases. 

4.5.3  Numerical  Efficiency  Calculations  for  90°  Elbow  Region 
of  Aerosol  Sampling  Train 

4. 5. 3.1  Monodispersed  Particle  Distribution  -  90°  Elbow  Region 

Table  11  presents  the  three  model's  numerical  results 
for  the  90°  elbow  region  of  the  sampling  train.  Figure  17 
provides  a  graphical  presentation  of  the  numerical  results.  It 
is  important  to  note  that  for  these  calculations  the  90°  elbow 
region  is  not  affected  by  the  inlet's  orientation  to  the  air 
stream. 


The  90°  elbow  region's  sampling  efficiency  results 
indicate  the  following: 

•  All  three  models  predict  a  decrease  in  the  efficiency 
as  the  particle  size  increases. 

•  All  three  models  predict  a  decrease  in  the  efficiency 
as  the  flow  rate  increases. 

•  The  Battelle  model  predicts  the  highest  efficiencies. 
The  Texas  A&M  model  and  the  University  of  Duisburg 
model  predict  identical  but  lower  efficiencies. 
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MASS  MEDIUM 

INLET  PROBE 

FLOW  RATE 

SAMPLING  EFFICIENCY  f%  mas8>  FOR  HORIZONTAL  SECTION 

AERODYNAMIC 
DIAMETER 
(MMADH  . II 

ORIENTAIION  1  QJ 

n/min> 

TEXAS  A  t  M 

MHELLE 

UNIVERSITY  OF 
DUiSBURO 

1.  5 

0° 

70 

99.0 

98.5 

99.96 

2.  15 

0° 

70 

91  .a 

88.0 

98.7 

3.  5 

0° 

130 

99.4 

99.2 

99.8 

4.  15 

0® 

130 

92.6 

93.2 

91.9 

5.  5 

90° 

70 

99.1 

98.5 

N/A 

6.  15 

90° 

70 

93.9 

88.0 

N/A 

7.  5 

90° 

130 

99.5 

99.2 

N/A 

8.  15 

90° 

130 

94.4 

93.2 

N/A 

Table  10.  Model  Efficiency  Calculations  (Polydispersed  Particle  Sizes) 
for  Horizontal  Region  of  Aerosol  Sampling  Train 
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PARTICLE 

FLOW  RATE 

SAMPLING  EFFICIENCY  (%  mass)  FOR  90^ 

ELBOW  REGION 

SIZE  (  M  > 

<l/mln) 

lEXASAAM 

BAnELLE 

UNIVERSITY  OF 
DUISBURG 

1. 

3 

70 

99.0 

99.8 

99.0 

2. 

5 

70 

97.3 

99.8 

97.3 

3. 

7 

70 

94.9 

99.7 

94.9 

4. 

10 

70 

89.9 

99.6 

89.9 

5. 

15 

70 

78.9 

99.4 

78.9 

6. 

20 

70 

65.7 

98.8 

65.7 

7. 

25 

70 

51.9 

97.5 

51.9 

8. 

3 

130 

98J2 

99.4 

98.2 

9. 

5 

130 

95.1 

99.3 

95.1 

10. 

7 

130 

90.8 

99.1 

90.7 

11. 

10 

130 

82.1 

98.5 

82.1 

12. 

15 

130 

64.4 

94.2 

64.4 

13. 

20 

130 

45.8 

89.1 

45.8 

14. 

25 

130 

29.6 

84.1 

29.6 

Table  1 1 .  Numerical  Model  Calculations  (Monodispersed  Particle  Sizes) 
for  90°  Elbow  Region  of  Aerosol  Sampling  Train 
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4. 5. 3. 2  Polydispersed  Particle  Distribution  -  90°  Elbow  Region 

Table  12  presents  the  three  model's  numerical  results 
for  the  90°  elbow  region  of  the  sampling  train.  Figures  18  and 
19  provide  a  graphical  presentation  of  the  numerical  results.  It 
is  important  to  note  that  only  the  Texas  A&M  model  considers  air 
stream  orientations  other  than  parallel  to  the  inlet.  The  German 
model  only  considers  a  0°  inlet  orientation.  The  Battelle  model 
does  not  consider  inlet  orientation. 

The  90°  elbow  sampling  efficiency  results  for  a 
polydispersed  particle  distribution  indicate  the  following: 

•  All  three  models  predict  a  decrease  in  the  sampling 
efficiency  as  the  particle  MMAD  increases. 

•  All  three  models  predict  a  decrease  in  the  sampling 
efficiency  as  the  flow  rate  increases. 

•  The  Texas  A&M  model  predicts  an  increase  in  the 
sampling  efficiency  as  the  inlet  orientation  increases 
from  0°  to  90°. 

4.5.4  Mumerlcal  Efficiency  Calculations  for  Vertical  Region 

of  Aerosol  Sampling  Train 

4. 5. 4.1  Monodispersed  Particle  Distribution  -  Vertical  Region 

Table  13  presents  the  three  model's  numerical  results 
for  the  vertical  region  of  the  sampling  train.  Figure  20 
provides  a  graphical  presentation  of  the  numerical  results.  It 
is  important  to  note  that  for  these  calculations,  the  vertical 
region  is  not  affected  by  the  inlet's  orientation  to  the  air 
stream. 

The  vertical  region's  sampling  efficiency  results 
indicate  the  following: 

•  The  Texas  A&M  model  and  the  German  model  predict  a 
decrease  in  the  efficiency  as  the  particle  size 
increases.  The  Battelle  model  predicts  no  change  in 
the  efficiency  as  the  particle  size  increases. 

•  The  Texas  A&M  model  and  the  German  model  predict  a 
decrease  in  the  efficiency  as  the  flow  rate  increases. 
The  Battelle  model  predicts  no  change  in  the  efficiency 
as  the  particle  size  increases. 

4. 5. 4. 2  Polydispersed  Particle  Distribution  -  Vertical  Region 

Table  14  presents  the  three  model's  numerical  results 
for  the  vertical  region  of  the  sampling  train.  Figures  21  and  22 
provide  a  graphical  presentation  of  the  numerical  results.  It 
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MASS  MEDIUM 

INLET  PROBE 

FLOW  RATE 

SAMPLING  EFFICIENCY  <%  massi  FOR  90° 

ELBOW  SECTION 

AERODYNAMIC 
DIAMETER 
IMMADW  fJ.  > 

ORIENTATION  (  O  1 

U/min> 

TEXAS  A  A  M 

BAHEUS 

UNIVERSITY  OF 
DUISBURG 

1.  5 

0° 

70 

96.5 

99.7 

96.4 

2.  15 

0= 

70 

76.6 

98.9 

74.7 

3.  5 

0^= 

130 

93.8 

99.2 

93.5 

4.  15 

0° 

130 

63.8 

93.5 

60.5 

5.  5 

90“ 

70 

96.7 

99.7 

N/A 

6.  15 

90° 

70 

81.1 

98.9 

N/A 

7.  5 

90° 

130 

94.0 

99.2 

N/A 

8.  15 

90° 

130 

68.2 

93.5 

N/A 

Table  12.  Moddl  Efficiency  Calculations  (Polydispersed  Particle  Sizes) 
for  90°  Elbow  Region  of  Aerosol  Sampling  Train 
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PARTICLE 

SIZE  (  U  ^ 

FLOW  RATE 
iizmioi 

SAMPLING  EFFICIENCY  (%  mass)  FOR  VERTICAL  REGION 

TEXAS  AiM  BATTELLE  UNIVERSITY  OF 

DUISBURG 

1. 

3 

70 

100.0 

100.0 

100.0 

2. 

5 

70 

99.99 

100.0 

99.99 

3. 

7 

70 

99.98 

100.0 

99.98 

4. 

10 

70 

99.9 

100.0 

99.9 

5. 

IS 

70 

99.6 

100.0 

99.6 

6. 

20 

70 

98.8 

100.0 

98.8 

7. 

25 

70 

96.9 

100.0 

97.2 

8. 

3 

130 

99.99 

100.0 

99.99 

9. 

5 

130 

99.96 

100.0 

99.96 

10. 

7 

130 

99.9 

100.0 

99.9 

11. 

10 

130 

99.4 

100.0 

99.4 

12. 

15 

130 

96.8 

100.0 

97.1 

13. 

20 

130 

84.4 

100.0 

91.0 

14. 

25 

130 

73.7 

100.0 

79.6 

Table  13.  Numerical  Model  Calculations  (Monodispersed  Particle  Sizes) 
for  Vertical  Region  of  Aerosol  Sampling  Train 
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MASS  MEDIUM 

INLET  PROBE 

FLOW  RATE 

SAMPLING  EFFICIENCY  massi  FOR  VERTICAL  SEC TIO N 

AERODYNAMIC 

DIAMETER 

ORIENTATION  ( O  1 

ll/min) 

TEXASAAM 

BATTELLE 

UNIVERSITY  OF 

} 

DUISBURG 

1.  5 

0° 

70 

99.98 

100.0 

99.96 

2.  15 

0° 

70 

99.3 

100.0 

98.7 

3.  S 

0° 

130 

99.9 

100.0 

99.6 

4.  15 

0° 

130 

95.7 

100.0 

91.9 

5.  5 

90° 

70 

99.98 

100.0 

N/A 

6.  15 

90° 

70 

99.6 

100.0 

N/A 

7.  5 

90° 

130 

99.9 

100.0 

N/A 

8.  15 

90° 

130 

96.7 

100.0 

N/A 

Table  14.  Model  Efficiency  Calculations  (Polydispersed  Particle  Sizes) 
for  Vertical  Region  of  Aerosol  Sampling  Train 


70 


Pairc«nt  P®notrottc:n 


Figure  22.  Penetration  Efficiency 


Polydispersed  -  Vertical  -  90*  Inlet 


Mass  Median  Aeroc^nomic  Diameter  (/;m) 

0  70  L/min  Texas  A&M  +  130  L/min  Texas  A&M  0  70  L/min  Battelle 

ii  130  L/min  Bottelie 
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is  important  to  note  that  only  the  Texas  A&M  model  considers 
inlet  orientations  other  than  parallel  to  the  air  stream.  The 
German  model  only  considers  a  0°  inlet  orientation.  The  Battelle 
model  does  not  consider  inlet  orientation. 

The  vertical  region  sampling  efficiency  results  for  a 
polydispersed  particle  distribution  indicate  the  following: 

•  The  Texas  A&M  and  German  models  predict  a  decrease  in 
the  sampling  efficiency  as  the  particle  MMAD  increases. 
The  Battelle  model  predicts  no  change  in  efficiency. 

•  The  Texas  A&M  and  German  models  predict  a  decrease  in 
the  sampling  efficiency  as  the  flow  rate  increases. 

The  Battelle  model  predicts  no  change  in  efficiency. 

•  The  Texas  A&M  model  predicts  an  increase  in  the 
sampling  efficiency  as  the  inlet  orientation  increases 
from  0°  to  90°. 

4.5.5  Numerical  Efficiency  Calculations  for  45°  Elbow  Region 

of  Aerosol  Sampling  Train 

4. 5. 5.1  Monodispersed  Particle  Distribution  *-  45°  Elbow  Region 

Table  15  presents  the  two  model's  numerical  results  for 
the  45°  elbow  region  of  the  sampling  train.  Figure  23  provides  a 
graphical  presentation  of  the  numerical  results.  It  is 
important  to  note  that  for  these  calculations,  the  45°  elbow 
region  is  not  affected  by  the  inlet's  orientation  to  the  air 
stream.  The  45°  elbow  region's  sampling  efficiency  results 
indicate  the  following: 

•  Both  models  (Texas  A&M  and  Battelle)  predict  a  decrease 
in  the  efficiency  as  the  particle  size  increases.  The 
Texas  A&M  model  predicts  a  greater  decrease  in 
efficiency  than  does  the  Battelle  model  at  the  larger 
particle  sizes. 

•  Both  models  predict  a  decrease  in  the  efficiency  as  the 
flow  rate  increases.  The  Texas  A&M  model  predicts  a 
greater  decrease  in  efficiency  than  does  the  Battelle 
model  at  the  larger  particle  sizes. 

4. 5. 5. 2  Polydispersed  Particle  Distribution  -*  45°  Elbow  Region 

Table  16  presents  the  three  model's  numerical  results 
for  the  45°  elbow  region  of  the  sampling  train.  Figures  24  and 
25  provide  a  graphical  presentation  of  the  numerical  results.  It 


PARTICLE 

FLOW  RATE 

SAMPLING  EFFICIENCY  (%  massi  FOR  45^  ELBOW  REGION 

SIZEUi  \ 

fl/mlnl 

TEXAS. A.  AM 

BATTELLE 

1. 

3 

70 

99.5 

99.9 

2. 

S 

70 

98.7 

99.9 

3. 

7 

70 

97.4 

99.9 

4. 

10 

70 

94.8 

99.8 

5. 

15 

70 

88.8 

99.7 

6. 

20 

70 

81.0 

99.4 

7. 

25 

70 

72.0 

98.8 

8. 

3 

130 

99.1 

99.7 

9. 

5 

130 

97.5 

99.7 

10. 

7 

130 

95.3 

99.6 

11. 

10 

130 

90.6 

99.2 

12. 

15 

130 

80.2 

97.1 

13. 

20 

130 

87.7 

94.4 

14. 

25 

130 

54.4 

91.7 

Table  15.  Numerical  Model  Calculations  (Monodispersed  Particle  Sizes) 
for  45^  Elbow  Region  of  Aerosol  Sampling  Train 
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□  70  L/min  Texos  Ai«M  +  130  L/min  Texas  AisM  0  70  L/min  Battelle 


6  130  L/min  Bottelle 
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MASS  MEDIUM 

INLET  PROBE 

FLOW  RATE 

SAMPLINCf  FFFICIENCY  f%  massi  FOR  45’^ 

ELBOW  SECTION 

AERODYNAMIC 

ORIENTATION  <  O'  ) 

liymln) 

DIAMETER 
(MMAPK  H-  1 

TEXASA&M 

BAHELLE 

UNIVERSITY  OF 
DUISBURG 

1.  s 

0° 

70 

98.3 

99.9 

96.4 

2.  15 

0° 

70 

88.9 

99.5 

74.7 

3.  5 

0° 

130 

96.9 

99.6 

93.5 

4.  15 

0° 

130 

83.6 

96.8 

60.5 

5.  5 

90° 

70 

98.4 

99.9 

N/A 

6.  15 

90° 

70 

90.9 

99.5 

N/A 

7.  5 

90° 

130 

97.0 

99.6 

N/A 

8.  15 

90° 

130 

85.2 

96.8 

N/A 

Table  1 6.  Model  Efficiency  Calculations  (Polydispersed  Particle  Sizes) 
for  45°  Elbow  Region  of  Aerosol  Sampling  Train 


is  important  to  note  that  only  the  Texas  A&M  model  considers  air 
stream  orientations  other  than  parallel  to  the  inlet.  The  German 
model  only  considers  a  0°  inlet  orientation.  The  Battelle  model 
does  not  consider  inlet  orientation. 

The  45°  elbow  sampling  efficiency  results  for  a 
polydispersed  particle  distribution  indicate  the  following; 

•  All  three  models  predict  a  decrease  in  the  sampling 
efficiency  as  the  particle  MMAD  increases. 

•  All  three  models  predict  a  decrease  in  the  sampling 
efficiency  as  the  flow  rate  increases. 

•  The  Texas  A&M  model  predicts  an  increase  in  the 
sampling  efficiency  as  the  inlet  orientation  increases 
from  0°  to  90°. 

4.5.6  Mximerical  Efficiency  Calculations  for  Inclined  Region 

of  Aerosol  Sampling  Train 

4.5. 6.1  Monodispersed  Particle  Distribution  -  Inclined  Region 

Table  17  presents  the  three  model's  numerical  results 
for  the  inclined  region  of  the  sampling  train.  Figure  26 
provides  a  graphical  presentation  of  the  numerical  results.  It 
is  important  to  note  that  for  these  calculations,  the  vertical 
region  is  not  affected  by  the  inlet's  orientation  to  the  air 
stream. 

The  inclined  region's  sampling  efficiency  results 
indicate  the  following: 

•  All  three  models  predict  a  decrease  in  the  efficiency 
as  the  particle  size  increases. 

•  The  Battelle  and  German  models  predict  an  increase  in 
the  efficiency  as  the  flow  rate  increases.  In  addition 
the  German  model  predicts  a  decrease  in  the  sampling 
efficiency  at  the  largest  particle  size  (25/xin)  .  The 
Texas  A&M  model  also  predicts  an  increase  in  the 
efficiency  as  the  flow  rate  increases  except  at  the  two 
largest  particle  sizes  (20  and  25/im)  . 

4. 5. 6. 2  Polydispersed  Particle  Distribution  ~  Inclined  Region 

Table  18  presents  the  three  model's  numerical  results 
for  the  inclined  region  of  the  sampling  train.  Figures  27  and  28 
provide  a  graphical  presentation  of  the  numerical  results.  It  is 
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PARTICLE 

FLOW  RATE 

SAMPLING  EFFICIENCY  (%  mass!  FOR  INCLINED  REGION 

S.lZEiJi.I 

(l/mln) 

TEXAS  A  ft  M 

SAHELLE 

UNIVERSITY  OF 
DUISBURG* 

1. 

3 

70 

99.7 

99.7 

99.7 

2. 

5 

70 

99.3 

99.2 

99.2 

3. 

7 

70 

98.5 

98.5 

98.5 

4. 

10 

70 

97.1 

9f..99 

97.0 

5. 

15 

70 

93.4 

93.3 

93.4 

6. 

20 

70 

88.4 

88.5 

88.5 

7. 

25 

70 

82.0 

82.6 

82.2 

8. 

3 

130 

99.9 

99.8 

99.9 

9. 

5 

130 

99.6 

99.6 

99.6 

10. 

7 

130 

99.1 

99.2 

99.1 

11. 

10 

130 

98.1 

98.4 

98.1 

12. 

15 

130 

94.8 

96.4 

95.0 

13. 

20 

130 

84.1 

93.6 

90.7 

14. 

25 

130 

73.7 

90.2 

79.6 

*  Elflclancy  caleulationa  baaad  on  modal  modificationa  provided  by  General  Management  Aasociatea 


Table  17.  Numerical  Model  Calculations  (Monodispersed  Particle  Sizes) 
for  inclined  Region  of  Aerosol  Sampling  Train 
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Parcont  Ponelratton 


100 

98 

96 

94 

92 

90 

88 

86 

84 

82 

80 

78 

76 

74 

72 

70 

□ 
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MASS  MEDIUM 

INLET  PROBE 

FLOW  RATE 

SAMPLING  EFFICIENCY  (%  mass)  FOR  INCLINED  SECTION 

AERODYNAMIC 
DIAMETER 
(MMADK  P-  ) 

ORIENTATION  IG  1 

il'mini 

TEXASA&M 

BATTELLE 

UNIVERSITY  OF 
DUISBURG 

1. 

5 

0^ 

70 

99.0 

98.96 

98.9 

2. 

15 

0° 

70 

93.7 

92.2 

91.0 

3. 

5 

0° 

130 

99.4 

99.4 

99.4 

4. 

15 

0° 

130 

95.4 

95.6 

90.8 

5. 

5 

90'^ 

70 

99.1 

98.96 

N/A 

6. 

15 

90° 

70 

94.9 

92.2 

N/A 

7. 

5 

90° 

130 

99.5 

99.4 

N/A 

8. 

15 

90° 

130 

96.1 

95.6 

N/A 

Table  18.  Model  Efficiency  Calculations  (Poiydispersed  Particle  Sizes) 
for  Inclined  Region  of  Aerosol  Sampling  Train 
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Mass  Median  Aerodynamic  Diameter  l>ym) 

70  L/min  Texas  A4M  +  130  L/min  Texas  AijM  0  70  L/min  Batteile 

130  L/min  Botlelle  X  70  L/min  Duisburg  7  130  L/min  Duisburg 


Purcent  ■tratlan 


Figure  28.  Penetration  Efficiency 


Potydisperssd  —  Incline  -  9Cr  Inlet 


Mass  Median  Aerodji^ncmic  Diameter  (^) 

0  70  L/min  Texas  A&M  +  130  L/min  Texas  AdcM  0  70  L/min  Battelte 

A  130  L/min  Bottelle 


8 


important  to  note  that  only  the  Texas  A&M  model  considers  inlet 
orientations  other  than  parallel  to  the  air  stream.  The  German 
model  only  considers  a  0°  inlet  orientation.  The  Battelle  model 
does  not  consider  inlet  orientation. 


The  inclined  sampling  efficiency  results  for  a 
polydispersed  particle  distribution  indicate  the  following: 


4.5.7 


All  three  models  predict  a  decrease  in  the  sampling 
efficiency  as  the  particle  MMAD  increases. 

All  three  models  predict  an  increase  in  the  sampling 
efficiency  as  the  flow  rate  increases.  However,  the 
German  model  predicts  a  decrease  in  the  efficiency  as 
the  flow  rate  increases  for  the  ISfim  MMAD. 

The  Texas  A&M  model  predicts  an  increase  in  the 
sampling  efficiency  as  the  inlet  orientation  increases 
from  0°  to  90°. 

cumulative  Efficiency  calculations  for  Overall  Sampling 
Train 


4.5.7. 1  Monodispersed  Particle  Distribution  -  Overall 
Efficiencies 


Table  19  presents  the  three  model's  numerical 
efficiency  results  for  the  overall  sampling  train.  Figure  29 
provides  a  graphical  presentation  of  the  numerical  results  for 
the  monodispersed  conditions.  It  is  important  to  note  that  for 
these  calculations,  the  Battelle  model  does  not  consider  inlet 
orientation,  the  German  model  considers  only  a  0°  inlet 
orientation  and  the  Texas  A&M  model  considers  variable  inlet 
orientations . 


The  cumulative  efficiency  results  for  monodispersed 
aerosol  particles  penetrating  through  the  overall  sampling  train 
indicate  the  following; 

•  All  three  models  predict  a  decrease  in  the  cumulative 
efficiency  as  the  particle  size  increases. 

•  Both  the  Texas  A&M  and  University  of  Duisburg  models 
predict  a  decrease  in  the  cumulative  sampling 
efficiency  as  the  flow  rate  increases.  The  Battelle 
model  overall  predicts  a  decrease  in  the  cumulative 
efficiency  as  the  flow  rate  increases  except  at  the  5 
and  7^m  particle  sizes  where  an  increase  occurs  in  the 
efficiency. 
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TYPE  OF 
^  PARTICLE  ^ 

PARTICLE 

SiZEliLl 

INLET  PROBE 
ORIENTATIPN(ft) 

FLOW  RATE 
<l/mlnl 

AIR 

VELOCITY 

(m/iee) 

CUMULATIVE  SAMPLING  EFFICIENCY  (%  mas5- 

TEXAS  A  &M  BATTELLE  UNIVERSITY  OF 

1.MONODISPERSEO 

3 

0° 

70 

3 

97.4 

99.0 

DUISBURG 

98.4 

2.  MONODISPERSED 

5 

0° 

70 

3 

93.6 

97.8 

95.7 

3.  MONODISPERSED 

7 

0° 

70 

3 

88.4 

96.0 

91.7 

4.  MONODISPERSED 

10 

0° 

70 

3 

78.6 

92.4 

83.3 

5.  MONODISPERSED 

15 

0° 

70 

3 

59.0 

83.9 

64.6 

6.  MONODISPERSED 

20 

0° 

70 

3 

39.2 

73.1 

43.8 

7.  MONODISPERSED 

25 

0° 

70 

3 

22.8 

60.7 

25.6 

8.  MONODISPERSED 

3 

0° 

130 

3 

96.3 

98.8 

95.9 

9.  MONODISPERSED 

S 

0° 

130 

3 

90.4 

98.0 

89.1 

10.  MONODISPERSED 

7 

0° 

130 

3 

82.5 

96.7 

79.8 

11.  MONODISPERSED 

10 

0° 

130 

3 

67.8 

93.9 

63.0 

12.  MONODISPERSED 

15 

0° 

130 

3 

41.0 

83.7 

34.6 

13.  MONODISPERSED 

20 

0° 

130 

3 

16.6 

71.7 

14.2 

14.  MONODISPERSED 

25 

0° 

130 

3 

5.6 

60.2 

3.9 

15.  POLYDISPERSED 

MMAD=S 

(  a  =1.5) 

0° 

70 

3 

91.9 

97.1 

94.1 

16.  POLYDISPERSED 

MMAO=15 

(  cr  =1.5) 

0° 

70 

3 

S5.S 

79.8 

59.9 

17.  POLYDISPERSED 

MMAD=S 
(  a  =1.5) 

0° 

130 

3 

87.8 

97.4 

86.0 

18.  POLYDISPERSED 

MMAD^IS 

(  cr  =1.5) 

0° 

130 

3 

40.1 

80.7 

35.8 

Table  19.  Cumulative  Sampling  Efficiency  Results  Generated  by  the 

Three  Aerosol  Models 
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Percent  Penetrotton 


.Figure  29.  Penetration  Efficiency 

Monodispersed  —  Cumulativ«  -  Qf  'nlet 

100 

90 
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40 

30 

20 
10 

0 

2  6  10  14  18  22  26 

Aerodynomic  Pariicie  Diameter  (/^m) 

0  70  L/min  Texas  AiM  +  130  L/min  Texas  AitM  0  70  L/min  Bottelle 

A  130  L/min  Bottelle  X  70  L/min  Duisburg  7  130  L/min  Duisburg 
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•  Both  the  Texas  A&M  and  German  models  produce  very 
comparable  predicted  results  across  all  seven  particle 
sizes  and  for  the  two  flow  rates.  The  Battelle  model 
as  compared  to  the  other  two  models  produces  very 
comparable  predicted  results  at  the  3  and  Sjim  particle 
sizes  but  produces  significantly  higher  predicted 
efficiencies  at  the  other  five  particle  sizes. 

4. 5. 7. 2  Polydispersed  Particle  Distribution  -  Overall 
Efficiencies 

Table  19  cited  earlier  also  presents  the  three  model's 
numerical  efficiency  results  for  the  overall  sampling  train  based 
on  a  polydispersed  particle  distribution.  Figure  30  provides  a 
graphical  presentation  of  the  numerical  results  for  the 
polydispersed  condition.  It  is  important  to  note  that  for  these 
calculations,  the  Battelle  model  does  not  consider  inlet 
orientation,  the  University  of  Duisburg  model  considers  only  a  0° 
inlet  orientation,  and  the  Texas  A&M  model  considers  variable 
inlet  orientations. 

The  cumulative  efficiency  results  for  polydispersed 
aerosol  particles  penetrating  through  the  overall  sampling  train 
indicate  the  following: 

•  All  three  models  predict  a  decrease  in  the  cumulative 
efficiency  as  the  particle  MMAD  increases. 

•  Both  the  Texas  A&M  and  German  models  predict  a  decrease 
in  the  cumulative  efficiency  as  the  flow  rate 
increases.  The  Battelle  model  predicts  an  increase  in 
the  cumulative  efficiency  as  the  flow  rate  increases. 

•  Both  the  Texas  A&M  and  German  models  produce  very 
comparable  predicted  results  across  the  two  particle 
MMADs  and  two  flow  rates.  The  Battelle  model  as 
compared  to  the  other  two  models  produces  fairly 
comparable  predicted  results  at  the  5/im  MMAD  and  70 
1/min  flow  rate.  However,  the  model  produces  higher 
predicted  efficiencies  at  che  other  three  combinations 
of  flow  rates  and  particle  MMAD. 
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2  6  10  14  18  22  26 


Mass  Median  Aero<^nomic  0i«meter  O^m) 

0  70  L/min  Texos  AicM  +  1 30  L/min  Texas  AicM  0  70  L/min  Battel  le 

&  130  L/min  Bottelle  X  70  L/min  Duisburg  V  130  L/min  Duisburg 
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SECTION  5.  RESULTS  OF  EXPERIMENTAL  VERIFICATION 


5.1  EXPERIMENTAL  APPARATUS  AND  METHODOLOGY  FOR  AEROSOL 

PENETRATION  MEASUREMENTS 

The  aerosol  sampling  setup  as  presented  earlier  in 
Figure  l  was  the  experimental  conf igureation  used  to  measure  the 
penetration  (losses)  of  monodispersed  aerosol  particles  through 
the  sampling  train.  Seven  monodispersed  particle  sizes  were 
generated  in  the  experiment.  They  were  3,  5,  7,  10,  15,  20  and 
25nra  aerodynamic  equivalent  diameter  (AED)  .  Two  inlet  configura¬ 
tions  were  considered.  They  were  0°  and  90°  whereby  the  inlet  of 
the  transport  system  is  parallel  to  the  wind  tunnel  air  stream 
and  perpendicular  to  the  air  stream,  respectively. 

The  aerosol  transport  system  was  constructed  from  a 
single  piece  of  1-inch  Schedule  40  PVC  pipe.  The  inlet  was 
chamfered  at  3  0°  relative  to  the  tube  axis  and  was  machined  to  a 
sharp  edge.  Elbows  were  formed  by  electrically  heating  the 
appropriate  region  of  the  tube,  carefully  bending  it,  and 
clamping  the  tube  into  a  mold  as  it  cooled.  The  mold  prevented 
flattening  of  the  tube  at  the  elbow.  The  elbow  curvature  of  the 
mold  was  patterned  from  electrical  conduit  bends. 

An  aerosol  wind  tunnel  with  a  basic  600  mm  x  600  mm 
cross  section  was  used  for  the  testing.  Monodisperse  aerosol 
droplets  were  created  with  a  vibrating  jet  atomizer  (Ref.  16)  at 
the  tunnel  entrance.  A  Stairmand  disc  is  used  to  create  large 
scale  mixing  at  the  tunnel  entrance  and  to  help  uniformize  the 
aerosol  concentration  p'^ofiles  across  the  center  2/3  of  the  wind 
tunnel.  Downstream  from  the  Stairmand  disc  is  a  fan  which 
further  stirs  the  aerosol  to  help  obtain  uniform  concentration 
profiles.  A  perforated  plate  is  placed  downstream  from  the  fan 
to  reduce  the  scale  of  turbulence  and  to  uniformize  the  velocity 
profile.  At  the  test  section,  the  wind  tunnel  body  is  expanded 
to  a  cross  section  of  1.2  m  x  1.2  m  and  the  air  stream  is  reduced 
in  cross  section  to  360  mm  x  360  mm.  This  arrangement  reduces 
blockage  effects  that  could  arise  from  the  presence  of  artifacts 
in  the  test  section. 

Aerosol  was  generated  from  a  mixture  of  oleic  acid  and 
an  analytical  tracer  (sodium  fluorescein)  dissolved  in  alcohol. 
Immediately  after  generation,  the  alcohol  evaporates,  leaving  a 
residual  droplet  that  consists  of  12.4%  (m/m)  analytical  tracer 
and  oleic  acid.  Particle  size  of  the  aerosol  was  determined 
microscopically  with  the  droplets  collected  on  glass  slides  which 
had  been  treated  with  an  oil  phobic  agent.  Size  observed  under 
the  microscope  was  converted  to  an  equivalent  sperical  diameter 
through  use  of  the  flattening  factor  of  Olan-Figuroa  et  al.  (Ref. 
17) .  This  factor  takes  into  account  the  phenomenon  that  gravity 
flattens  droplets  on  the  slide  and  causes  them  to  appear  larger 


than  they  were  in  the  aerosol  state.  The  equivalent  spherical 
diameter  was  converted  to  aerodynamic  equivalent  diameter  (AED) 
through  use  of  the  expression: 

Ca  =  Cp  Pp  Dp'  (41) 

Where:  =  Cunningham's  slip  correction  for  the  aerodynamic 

equivalent  diameter  size,  =  AED,  =  density  of  water  (1000 
kg/m^)  ,  Cp  ==  Cunningham's  slip  correction  for  the  equivalent 
spherical  diameter.  Dp  =  equivalent  spherical  diameter,  and  Pp  = 
particle  density.  Cunningham's  correction  C,  for  the  two  sizes 
was  calculated  from: 

2.52  L 

=  1  + -  (42) 

D,- 


Where:  L  =  mean  free  path  of  the  air  molecules  (taken  as  0.068 

Min)  .  The  aerosol  density  is  kg/m^  as  calculated  from  the  mass 
fractions  of  oleic  acid  and  sodium  fluorescein. 

In  conducting  a  test,  the  aerosol  would  be  sampled  in 
the  wind  tunnel  with  an  isokinetic  probe  which  consists  of  a 
sharp  edged  inlet,  a  conical  diffuser  and  a  filter  collector. 
Either  in  parallel  with  operation  of  the  isokinetic  probe  or 
subsequent  to  use  of  the  isokinetic  probe,  we  would  draw  an 
aerosol  sample  through  the  model  transport  system.  Aerosol  which 
was  transmitted  through  the  system  was  collected  with  a  filter. 

In  an  analysis  laboratory,  sodium  fluorescein  was 
extracted  from  the  filter  of  the  isokinetic  probe,  from  the 
internal  walls  of  the  isokinetic  probe,  and  from  the  filter  of 
the  aerosol  transport  system.  The  concentration  of  fluorescein 
in  the  extracts  was  determined  with  a  Sequoia-Turner  Model  450 
Fluorometer  (Sequoia-Turner  Corp. ,  Mountain  View,  CA) . 

Penetration  of  aerosol  through  the  transport  system  was 
calculated  from: 


Where:  0-^^  =  concentration  of  aerosol  in  the  wind  tunnel  based 

on  the  fluorescein  collected  by  the  isokinetic  filter  and  the 
fluorescein  deposited  on  the  internal  walls  of  the  isokinetic 
probe,  =  concentration  of  aerosol  calculated  from  the 
fluorescein  collected  by  the  filter  at  the  downstream  end  of  the 
transport  system. 
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5 . 2  EXPERIMENTAL  RESULTS 

Tests  were  conducted  with  a  range  of  monodispersed 
aerosol  particle  sizes  of  approximately  3-25Mm  AED,  with  the 
inlet  of  the  transport  system  parallel  to  the  wind  tunnel  air 
stream  (0°)  and  with  the  inlet  perpendicular  to  the  air  stream 
(90°)  .  Tests  were  conducted  with  two  flow  rates  through  the 
system  (70  and  130  1/min) .  The  cumulative  sampling  efficiency 
results  from  these  tests  are  presented  in  Table  20  and  Figures  31 
and  32.  The  complete  set  of  experimental  data  are  provided  at 
Appendix  F.  The  cumulative  sampling  efficiency  results  for  a 
monodispersed  particle  distribution  indicate  the  following: 

•  As  the  particle  size  increases,  the  sampling 
efficiency  decreases. 

•  As  the  flow  rate  (70  1/min  to  130  1/min) 
increases,  the  sampling  efficiency  decreases. 

•  The  larger  the  particle  size,  the  greater  the 
decrease  in  sampling  efficiency. 

•  As  the  inlet  orientation  increases,  the  sampling 
efficiency  appears  to  decrease  based  on  a 
comparison  of  similar  sized  particles  across  the 
two  inlet  orientations. 

•  Increase  in  the  inlet  orientation  appears  to 
produce  a  more  dramatic  decrease  in  the  sampling 
efficiency  than  does  an  increase  in  the  flow  rate 
based  on  a  comparison  of  similiar  sized  particles 
across  the  two  inlet  orientations. 

At  the  completion  of  the  regular  test  program,  the 
sampling  train  was  to  be  segmentized  in  straight  sections, 
elbows,  and  the  inlet.  The  system  was  then  to  be  reassembled 
using  couplings  in  order  to  run  two  tests,  one  at  a  flow  rate  of 
130  1/rain  and  one  at  a  flow  rate  of  70  1/min,  in  which  the  wall 
deposits  would  be  extracted  for  the  individual  components  and 
determine  the  regional  losses.  The  particle  size  to  be  cnosen 
for  these  tests  were  those  which  give  penetration  values  on  the 
order  of  50%  for  the  overall  system.  This  phase  of  testing  was 
never  performed  due  to  time  and  money  constraints.  Therefore,  no 
experimental  data  was  gathered  to  determine  the  penetration  of 
aerosols  through  each  of  the  tubing  sections  of  the  sampling 
train.  In  addition,  no  experiments  were  conducted  which  utilized 
polydispersed  particle  distributions. 
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SIZE  OF 

MONOOISPERSED 
PARTICLE!  M  > 


INLET 

ORIENTATION  ^ 

f  «  ) 


FLOW  RATE 
fl/mln> 


CUMULATIVE  SAMPLING  PENETRATION 
EFFICIENCY  (%  masil 


3.0 

0° 

70 

105.3 

5.1 

0° 

70 

96.3 

7.0 

0° 

70 

89.9 

9.8 

0° 

70 

83.0 

12.0 

0° 

70 

73.7 

15.2 

0° 

70 

55.8 

19.4 

0° 

70 

32.8 

26.7 

0° 

70 

6.5 

3.0 

0° 

130 

93.1 

5.1 

0° 

130 

90.4 

7.0 

0° 

130 

86.4 

9.8 

0° 

130 

69.6 

12.0 

0° 

130 

57.6 

15.2 

0° 

130 

31.5 

19.4 

0° 

130 

11.1 

26.7 

0° 

130 

1.1 

3.5 

90° 

70 

93.1 

4.9 

90° 

70 

92.5 

7.2 

90° 

70 

76.V 

7.4 

90° 

70 

75.4 

10.1 

90° 

70 

61.9 

11,5 

90° 

70 

44.7 

14,9 

90° 

70 

20.4 

19.6 

90° 

70 

4.3 

25.4 

90° 

70 

1.9 

3.5 

90° 

130 

92.4 

4.9 

90° 

130 

87.6 

72. 

90° 

130 

82.9 

7.4 

90° 

130 

81.8 

10.1 

90° 

130 

62.5 

11.5 

90° 

130 

56.0 

14.9 

90° 

130 

27.1 

19.6 

90° 

130 

4.0 

25.4 

90° 

130 

0.5 

^  Rctativc  to  tho  Wind  Tunnoi  Air  Flow  OIroctlon 


Table  20.  Cumulative  Penetration  Efficiency  Resuits  Derived 
Experimentally  by  Texas  A  &  M 
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Figure  31,  Experimental  Results 
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Figure  32.  Experimental  Results 
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Aerodynamic  Particle  Diameter  (irr\) 
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SECTION  6.  COMPARISON  OF  EXPERIMENTAL  AND  NUMERICAL  RESULTS 


6.1  LIMITATIONS  IN  COMPARISON  OF  EXPERIMENTAL  AND  NUMERICAL 

RESULTS 

As  discussed  earlier  in  Section  5,2,  only  the 
cumulative  penetration  efficiency  of  monodispersed  aerosol 
particles  through  the  sampling  train  was  measured  in  the 
experiments.  Therefore,  no  comparison  could  be  made  to  determine 
how  realistic  the  model's  predicted  penetration  efficiency 
through  each  section  of  the  sampling  train  was  to  experimental 
data.  Furthermore,  only  the  Texas  ASM  and  German  models 
predicted  cumulative  penetration  efficiencies  through  all  six 
regions  of  the  sampling  train.  The  Battelle  model  predicted 
cumulative  penetration  efficiencies  through  five  sections  (inlet 
region  was  n^t  included)  of  the  sampling  train  and  the  IITRI 
model  predicted  cumulative  penetration  efficiencies  through  one 
region  (inlet  only)  of  the  sampling  train.  Consequently, 
comparison  of  the  experimental  results  will  be  made  to  only 
numerical  results  for  the  Texas  A&M,  Battelle,  and  German  models. 
It  is  important  to  note  that  the  experimental  apparatus  was  not 
able  to  achieve,  for  each  experimental  run,  the  specified 
particle  size  to  be  disseminated  at  the  inlet  section  of  the 
train.  However,  the  model  results  are  based  on  a  specified 
particle  size  being  disseminated  at  the  inlet  section. 

Therefore,  the  model  predictions  are  compared  to  experimental 
data  which  may  be  for  similar  and  not  necessarily  identical 
particle  size  disseminated  in  the  wind  tunnel.  Comparison  of 
the  model  results  to  the  experimental  measurements  are  presented 
in  Table  21  and  Figures  33  and  34. 

6.2  DISCUSSION  OF  EXPERIMENTAL  VS  NUMERICAL  RESULTS 

6.2.1  Comparison  of  Results  for  o°  Inlet  Orientation  and  70 

1/min  Flow  Rate  Conditions 

The  Texas  A&M  model  and  the  University  of  Duisburg 
model  simulate  a  0°  inlet  orientation  and  70  1/min  flow  rate. 
However,  the  Battelle  model  docs  not  consider  different  inlet 
orientations  but  instead  assumes  100%  inlet  efficiency.  The 
Battelle  model  does  consider  different  flow  rates. 

Both  the  Texas  A&M  and  German  models  provide 
predictions  of  cumulative  penetration  efficiencies  for  the  seven 
monodispersed  particle  distributions  which  are  reasonably  close 
to  those  measured  in  the  wind  tunnel.  The  Battelle  model 
predicts  cumulative  efficiency  values  which  are  in  close 
agreement  with  the  measured  values  at  the  3.0  and  5.0  ^m  particle 
sizes.  However,  at  the  15,  20  and  25  ij-tr  particle  sizes,  the 
Battelle  model  predictions  differ  significantly  from  the  measured 
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SIZE  OF  INLET  FLOW  RATE  CUMULATIVE  SAMPLING  PENETRATION  EFFICIENCY  (%  maul 

MONOOiSPERSEO  ORIENTATION  ILmln) 

PARTICLE  I M )  (  a  )  EXPERIMENTALLY  MOPEL  PREPICTIONS 

MEASURED  . 

TEXAS-AAid  gEflMAN  °  BATTEU.E  ~ 


3.0  0°  70  105.3  97.4  98.4  99.0 


5.1 

0° 

70 

96.3 

93.6 

95.7 

97.8 

7.0 

0° 

70 

69.9 

88.4 

91.7 

96.0 

9.8 

70 

83.0 

78.6 

83.3 

924 

12.0 

0° 

70 

73.7 

— 

■ — 

— 

15.2 

0° 

70 

55.8 

59.0 

64.6 

83.9 

10.4 

0° 

70 

323 

393 

43.8 

73.1 

25.7 

70 

6.5 

22.8 

25.6 

60.7 

3.0 

0° 

130 

93.1 

96.3 

95.9 

98.8 

5.1 

0° 

130 

90.4 

90.4 

89.1 

98.0 

7.0 

0° 

130 

88.4 

823 

79.8 

96.7 

9.8 

0° 

130 

69.6 

67.8 

63.0 

93.9 

12.0 

0° 

130 

57.6 

— 

— 

— 

15.2 

0° 

130 

31.5 

41.0 

34.6 

83.7 

19.4 

0° 

130 

11.1 

16.8 

143 

7V7 

26.7 

0° 

130 

1.1 

5.6 

3.9 

603 

3.5 

90° 

70 

93.1 

94.3 

— 

— 

4.9 

90° 

70 

923 

85.6 

— 

— 

7Jl 

90° 

70 

76.9 

74.7 

— 

— 

7.4 

90° 

70 

75.4 

— 

— 

— 

10.1 

90° 

70 

61.9 

57,0 

— 

— 

11.5 

90° 

70 

44.7 

— 

— 

— 

14.9 

90° 

70 

20.4 

31.5 

— 

— 

19.6 

90° 

70 

4.3 

153 

— 

— 

25.4 

90° 

70 

1.9 

6.4 

— 

— 

33 

90° 

130 

923 

94.3 

— 

— 

4.9 

90° 

130 

873 

853 

— 

— 

7J2 

90° 

130 

823 

74.0 

— 

— 

7.4 

90° 

130 

813 

— 

— 

— 

10.1 

90° 

130 

623 

55.1 

— 

— 

113 

90° 

130 

583 

— 

— 

— 

14.9 

90° 

130 

27.1 

27.1 

— 

— 

19.6 

90° 

130 

4.0 

8.7 

— 

— 

253 

90° 

130 

03 

23 

— 

— 

*  TlM  BatteUa  MocM  doM  not  conaMar  InM  Ortontatlon.  Th«r«for*,  th*  °  ThOMprodletiorarotloetthooxtonsionmi^byQanorai 
prodictlom  moot  clooaly  eorrotpond  to  a  0°  InM  Orlontalion  Managomont  Aaaoeiataa  to  Induda  a  45°  band  and 

Inciinad  aaction  pradletiona 


Table  21.  Cumulative  Penetration  Efficiency  Resuits: 
Model  Predictions  vs.  Experimentai  Data 
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values,  the  model  overpredicts,  and  as  the  particle  size 
increases  the  discrepancy  (percentage  difference)  between  the 
predicted  and  measured  values  becomes  larger. 


6.2.2  Comparison  of  Results  for  0°  Inlet  Orientation  and  130 

1/min  Flow  Rate  Conditions 

The  Texas  A&M  model  and  the  University  of  Duisburg 
model  simulate  a  0°  inlet  orientation  and  130  1/min  flow  rate. 

The  Eattelle  model  only  simulates  the  130  1/min  flow  rate  and 
assumes  100%  inlet  efficiency. 

The  Texas  A&M  model  provides  predictions  of  cumulative 
penetration  efficiencies  for  the  3,  5,  7  and  10  isra  particle  sizes 
which  are  reasonably  close  to  penetration  efficiencies  measured 
in  the  wind  tunnel.  The  model's  cumulative  efficiency 
predictions  for  the  15,  20  and  25  /im  particle  sizes  differ 
significantly  from  the  measured  values,  are  over  predictions,  and 
as  the  particle  size  increases,  the  discrepancy  (percent 
difference)  between  the  predicted  and  measured  values  becomes 
larger. 


The  University  of  Duisburg  model  provides  predictions 
of  cumulative  penetration  efficiencies  for  the  3,  5,  7,  10  and  15 
(j.m  particle  sizes  which  are  reasonably  close  to  penetration 
efficiencies  measured  in  the  wind  tunnel.  The  model's  cumulative 
efficiency  predictions  for  the  20  and  25  particle  sizes 
differs  significantly  for  the  measured  values,  are  over 
predictions,  and  as  the  particle  size  increases,  the  discrepancy 
(percent  difference)  between  the  predicted  and  measured  values 
become  larger.  However,  it  is  interesting  to  note  that  for  these 
two  pai’ticle  sizes,  the  German  model's  predicted  values  are  in 
closer  agreement  to  the  measured  values  than  are  the  Texas  A&M 
model's  predictions. 

The  Battelle  model's  predictions  for  the  cumulative 
penetration  efficiencies  at  the  3  and  5  (j.m  particle  sizes  are  in 
agreement  to  penetration  efficiencies  measured  in  the  wind 
tunnel.  However,  the  model's  cumulative  efficiency  predictions 
for  the  7,  10,  15,  20  and  25  fim  particle  sizes  which  differ 
significantly  from  the  measured  values  are  overpredictions  and  as 
the  particle  size  increases,  the  discrepancy  (percent  difference) 
between  the  predicted  and  measured  values  becomes  larger. 

6.2.3  Comparison  of  Results  for  90°  Inlet  Orientation  and  70 

1/min  Flow  Rate  Conditions 

The  Texas  A&M  model  is  the  only  model  of  the  three 
which  provides  predictions  of  cumulative  penetration  efficiencies 
based  on  a  90°  inlet  orientation  and  70  1/min  flow  rate.  These 
predicted  results  are  in  close  agreement  to  the  efficiencies 
measured  in  the  wind  tunnel  at  the  3,  5,  7  and  10  urn  particle 
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sizes.  The  model  predictions  which  differ  significantly  from  the 
measured  values  at  the  15,  20  and  25  particle  sizes  are  over 
predictions  and  as  the  particle  size  increases,  the  discrepancy 
(percent  difference)  between  the  predicted  and  measured  values 
becomes  larger. 

6.2.4  Comparison  of  Results  for  90°  Inlet  Orientation  and  130 

1/min  Flow  Rate  Conditions 

The  Texas  A&M  model  is  the  only  model  of  the  three 
which  provides  predictions  of  cumulative  penetration  efficiencies 
based  on  a  90°  inlet  orientation  and  13  0  1/min  flow  rate.  These 
predicted  results  are  in  close  agreement  to  the  efficiencies 
measured  in  the  wind  tunnel  at  the  3,  5,  7,  10  and  15  /im  particle 
sizes.  The  model  predictions  which  differ  significantly  from  the 
measured  values  at  the  20  and  25  ^tm  particle  sizes  are  over¬ 
predictions  and  as  the  particle  size  increases,  the  discrepancy 
(percent  difference)  between  the  predicted  and  measured  values 
taecomes  larger. 
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SECTION  7 


CONCLUSIONS  AND  RECOMMENDATIONS 


7 . 1  CONCLUSIONS 

In  general,  aerosol  sampling  models  are  available  which 
can  be  used  to  design  aerosol  sampling  and  transport  systems. 

Four  candidate  models/algorithms  were  reviewed  which  showed  the 
greatest  potential  for  this  application.  They  were  the 
following: 

•  Texas  A&M  model 

•  Battelle  model 

•  IITRI  model 

•  Equations  compiled/developed  by  the  University  of 
Duisburg 

Each  candidate  model/algorithms  was  evaluated  relative 
to  an  aerosol  sampling  train  configuration  upon  which  the  model 
calculations  would  be  generated  and  compared  and  the  experimental 
data  would  be  compared.  The  sampling  train  consisted  of  six 
regions.  They  were  the  following:  (1)  inlet,  (2)  horizontal 
section,  (3)  vertical  section,  (4)  inclined  section,  (5)  90® 
elbow,  and  (6)  45®  elbow. 

7.1.1  Capabilities  of  Numerical  Models/Algorithms 

Of  the  four  models  evaluated,  the  Texas  A&M  model  was 
the  only  one  which  simulated  the  six  regions  of  the  sampling 
train.  The  Battelle  model  simulated  all  regions  of  the  sampling 
train  except  for  the  inlet  area.  For  this  model,  the  inlet 
efficiency  was  assumed  to  be  100%  for  all  inlet  orientations.  The 
University  of  Duisburg  equations  addressed  four  sections  of  the 
sampling  train.  However,  the  45°  elbow  and  inclined  sections 
were  incorporated  in  GMA's  version  of  the  University  of  Duisburg 
computer  program.  The  IITRI  model  only  simulated  the  inlet 
region.  Each  model  predicts  penetration  efficiency  through  each 
region  simulated  as  well  as  cumulative  penetration  efficiency 
through  all  regions  simulated.  Penetration  efficiency 
predictions  for  each  of  the  models  revealed  self-consistency, 

7.1.2  Numerical  Calculations  Generated  by  Four  Models 

The  Texas  A&M  model,  the  Battelle  model  and  the 
University  of  Duisburg  model  generated  penetration  efficiencies 
for  both  monodispersed  and  polydispersed  (log-normal)  particle 
distributions.  The  IITRI  model  generated  numerical  calculations 
for  only  monodispersed  particles. 
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The  penetration  efficiency  calculations  produced  by 
both  the  Texas  A&M  model  and  the  University  of  Duisbuurg  model 
are  very  comparable  across  the  monodispersed  particle 
distributions  and  the  polydispersed  distributions.  The  Battelle 
model  predicts  penetration  efficiency  values  comparable  to  those 
generated  by  the  Texas  A&M  and  University  of  Duisburg  models, 
however,  only  at  the  smaller  monodispersed  particle  size 
distributions.  For  the  larger  particle  sizes,  the  Battelle  model 
produces  significantly  higher  penetration  efficiency  values  than 
those  predicted  by  the  Texas  A&M  and  the  University  of  Duisburg 
models . 


The  IITRI  model  only  predicts  efficiency  values  at  the 
inlet  regions  of  the  sampling  train.  Results  from  the  IITRI 
model  compared  to  the  inlet  efficiency  results  produced  by  the 
Texas  A&M  model  and  the  University  of  Duisburg  model  indicate 
comparability  with  the  University  of  Duisburg  model  results  but 
noncomparability  with  the  Texas  A&M  model  results. 

7.1.3  Comparison  of  Experimental  Measurements  and  Mximerical 

Calculations 

Only  the  cumulative  penetration  efficiency  of 
monodispersed  aerosol  particles  through  the  sampling  train  was 
measured  in  the  experiments.  Consequently,  no  comparison  could 
be  made  to  determine  how  realistic  the  model's  predicted 
penetration  efficiency  through  each  section  of  the  sampling  train 
was  to  experimental  data.  Instead,  only  a  comparison  could  be 
made  between  predicted  and  measured  cumulative  penetration 
efficiency. 

Overall,  the  Texas  A&M  model  and  the  University  of 
Duisburg  model  predictions  produce  the  closest  agreement  to  the 
cumulative  penetration  efficiencies  measured  in  the  wind  tunnel. 
The  Battelle  model  predictions  on  the  other  hand,  produce  the 
greatest  difference  to  the  measured  values.  At  the  larger 
particle  sizes,  the  Battelle  model  significantly  overpredicts  the 
penetration  efficiency. 

7 . 2  RECOMHEMDAT lOMS 

The  primary  focus  of  this  study  was  not  only  to 
compare  numerical  calculations  among  the  four  candidate  models 
but  also  to  compare  these  model  results  to  the  experimental  data. 
Although  both  aspects  were  addressed,  a  comparison  of  model 
calculations  to  experimental  data  on  the  penetration  efficiency 
for  each  individual  section  of  the  aerosol  sampling  train  was 
never  gathered  due  to  time  and  money  constraints.  Without  this 
experimental  data,  a  determination  of  the  model's  realism  in 
predicting  penetration  efficiency  for  individual  sections  of  the 
sampling  train  could  not  be  made.  Therefore,  it  is  recommended 
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that  the  following  be  accomplished  in  order  to  identify  a 
comprehensive  model  that  already  exists  and  could  be  used  to 
design  aerosol  sampling  and  transport  systems: 

A.  Texas  A&M  conduct  wind  tunnel  experiments  to 
measure  monodispersed  particle  penetration  efficiencies  through 
individual  components  of  the  sampling  train.  Experimental 
conditions  identical  to  the  ones  used  in  this  study  should  be 
used  as  well  as  a  limited  set  of  additional  conditions. 

B.  The  Texas  A&M  model  as  the  preferred  model,  should 
be  used  to  compare  its  predicted  results  to  this  new  set  of 
experimental  data. 

C.  Further  investigation  should  be  given  to  extending 
the  University  of  Duisburg  equations  to  incorporate  theoreti¬ 
cally/empirically  based  equations  which  will  predict  penetration 
efficiencies  through  a  45°  elbow  section  and  an  inclined  section. 
Calculations  of  monodispersed  penetration  efficiencies  produced 
by  these  "new"  equations  for  single  sections  of  a  sampling  train 
should  be  compared  to  the  new  experimental  data. 

D.  Inlet  efficiency  predictions  for  the  IITRI  mode] 
should  be  compared  against  the  new  experimental  data. 

E.  Further  investigation  should  be  made  into 
understanding  why  the  Battelle  model,  when  compared  to 
experimental  data,  significantly  overpredicts  penetration 
efficiency  for  monodispersed  particle  sizes  in  the  10  to  25  nm 
region. 
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10  REM  **Th±s  program  calculates  the  "vertice  diameters"  of  a  ** 

20  REM  **log-normal  distribution  of  diameter  sizes,  using  ** 

30  REM  **equal  logarithmic  intervals.  This  program  reads  in  ** 

40  REM  **the  mass  median  diameter,  standard  deviation,  number  ** 

50  REM  **of  intervals,  and  the  name  of  the  output  file.  There** 

60  REM  **can  be  as  many  intervals  as  possible,  but  if  this  ** 

70  REM  **program  is  to  be  used  in  conjunction  with  the  program** 

80  REM  **"average-bas" ,  the  number  of  intervals  must  be  20-  ** 

90  INPUT"Mass  Median  Diameter  — >",  MMD 

100  INPUT "Standard  Deviation  - >",  SIGMA 

110  INPUT"Number  of  Intervals  - >",  N 

120  lNPUT"Diameter  Filename  - >",  LOGFILES 

130  OPEN  "a:"+LOGFILE$  FOR  OUTPUT  AS  #1 

140  PRINT#1,"MMD  •" ;MMD; "fm" ; "  Standard  Deviation  -";S1GMA; "£m" ; "  Size  [*]  fm" 
150  DMIN  -  MMD  *  (SIGMA  ‘  -3.62) 

160  DMAX  -  MMD  *  (SK^IA  *  3.62) 

170  INTLIN  -  (DMAX  -  DMIN)  /  N 

180  INTLOG  -  (LOG(DMAX)  -  LOG(DMIN))  /  (N  *  LOG(lO)) 

190  REM  **DIAMETERS  USING  EQUAL  LOGARITHMIC  INTERVALS** 

200  FOR  J  -  0  TO  N 

210  LD  -  (LOG(DMIN)  /  LOG(IO))  +  (J  *  INTLOG) 

220  D1  «  EXP(LD  *  LOG(IO)) 

230  PRINT#1,D1 
240  NEXT  J 
250  CL0SE#1 
260  END 
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10  REM  **This  program  calculates  the  geometric  average  particle** 
20  REM  **diameter  between  two  adjacent  vertice  diameters.  The  ** 
30  REM  **input  to  this  program  is  the  output  from  the  program  ** 
40  REM  **” diameter. bas " .  This  program  prompts  for  the  names  of** 
50  REM  **input  and  output  files.  ** 

60  DIM  D(21) 

70  DIM  GEOMEAN (20) 

80  INPUT"Input  Filename  - >",IFILE$ 

90  INPUT-Output  Filename  - >",0FILES 

100  OPEN  "a:"+IFILES  FOR  INPUT  AS  #1 

110  INPUT# 1, AS 

120  FOR  I  -  0  TO  20 

130  INPUT#1,  D{I) 

140  NEXT  I 
150  CLOSE  #1 

160  OPEN  "a:"+OFILES  FOR  OUTPUT  AS  #1 
170  PRINT#1,A$ 

180  FOR  1  -  1  TO  20 

190  GEOMEAN(I)  -  S0R(D(I-1)  *  D(I ) ) 

200  PRINT#1,  GEOMEAN(I) 

210  NEXT  I 
220  CLOSE  #1 


10  REM  **This  progratm  calculates  the  sampling  and  transport  efficiency** 

20  REM  **o£  an  aerosol.  The  program  prompts  for  the  input  filename,  ** 

30  REM  **output  filename  and  sampling  rate.  The  input  filename  is  the** 

40  REM  **name  of  the  data  file  that  contains  the  particle  sizes.  The  ** 

50  REM  **output  filename  is  the  name  of  the  file  that  the  efficiencies** 

60  REM  **are  output  to.  ** 

70  CLS 

80  REM  **VALUE  OF  PROCESS  CONSTANTS** 

90  REM  **AMBIENT  TEMPERATURE,  KELVINS** 

100  TEMP  -  298 

110  REM  **AMBIENT  PRESSURE,  CM  HG** 

120  PRESS  -  76 

130  REM  **AER0S0L  DENSITY,  GM/CM3** 

140  RHO  •  1 

150  REM  **PR0BE  DIAMETER,  INCHES** 

160  DS  -  1.05 

170  REM  **AIR  VISCOSITY,  KG/M/S** 

180  NU  -  .0000185 

190  REM  **AIR  DENSITY,  KG/M3** 

200  RHOAIR  •  1.1769 

210  REM  **H0RIZ0NTAL  TUBE  LENGTH,  FT** 

220  HLENGTH  -  2 

230  REM  **VERTICAL  TUBE  LENGTH,  FT** 

240  VLENGTH  «  2 

250  REM  **ACCELERATION  OF  GRAVITY,  ra/s2** 

260  G  -  9.810001 

270  REM  **INPUT  OF  PROCESS  VARIABLES** 

280  'WS  -  3 

290  INPUT  "Input  Filename  ";IFILES 
300  INPUT  "Output  Filename  ";OFILES 
310  INPUT  "Flow  Rate  (1/min)  ";VDOT 
320  OPEN  "a:"+IFILE$  FOR  INPUT  AS  #1 
330  INPUT#1,  AS 

340  OPEN  "a:"+OFILE$  FOR  OUTPUT  AS  #2 
350  PRINT#2,  AS 

360  PRINT#2, "Volumetric  Flow  Rate  «" ;VDOT; "1/min" 

370  PRINT#2,"  FRACTIONAL  PENETRATION  OF  AEROSOL" 

380  PRINT#2,"  SIZE  INLET  HORIZ.  90x  BEND  VERTICAL  45x  BEND  INCLINED  0 
VERALL" 

390  FOR  I  -  1  TO  20 

400  INPUT#1 ,  DP 

410  IF  DP  <-  0  THEN  GOTO  1120 

420  REM  **CALCULATION  OF  CUNNINGHAM  SLIP  CORRECTION** 

430  A  «  -.1095  *  PRESS  *  DP 
440  B  -  6.32  +  2.01  *  EXP(A) 

450  C  •  2  *  B  /  PRESS  /  DP 
460  CC  -  1  +  C 

470  RH4  **CALCULATION  OF  GRAVITATIONAL  SETTLING  VELOCITY** 

480  VS  -  .003  *  RHO  *  DP  *  DP 

490  REM  **CALCULATION  OF  BROWNIAN  DIFFUSION  COEFFICIENT,  MKS** 

500  DIFF  -  1.46529E-18  *  TEMP  *  CC  /  NU  /  DP 
510  REM  **CALCULAT10N  OF  SAMPLING  EFFICIENCY** 

520  STKW  -  4.37445E-09  *  CC  *  RHO  *  DP  *  DP  *  WS  /  NU  /  DS 

530  OMEGA  -  30.40245  *  WS  *  DS  *  DS  /  VDOT 

540  A  -  (OMEGA  /  STKW)  +  (2  *  OMEGA)  +  .62 

550  B  -  (2  *  OMEGA)  +  .62 

560  C  •  OMEGA  -  1 

570  SAMPEFF  -  1  +  (C  *  B  /  A) 

580  REM  **CALCULATION  OF  REYNOLD'S  NUMBER** 

590  NRE  -  8.35459E-04  *  VDOT  *  RHOAIR  /  DS  /  NU 
600  REM  **CALCULATION  OF  STOKES  NUMBER  FOR  BENDS** 

610  STK  -  1.438849E-10  *  CC  *  RHO  *  DP  *  DP  *  VDOT  /  NU  /  DS  /  DS  /  DS 
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620  IF  NRE  *<  2300  THEN  640 
630  IF  NRE  >  2300  THEN  740 

640  REM  **TRANSP0RT  EFFICIENCY  LOSS  BY  BROWNIAN  DIFFUSION,  LAMINAR  FLOW** 
650  XSI  =  57453.446#  *  DIFF  *  HLENGTH  /  VDOT 
650  IF  XSI  *<  .02  THEN  GOTO  1150 
670  IF  XSI  >  .02  THEN  GOTO  1210 

680  REM  **GRAVITATIONAL  SEDIMENTATION  LOSSES,  LAMINAR  FLOW** 

690  LENGTH  -  HLENGTH 
700  GOSUB  1270 

710  REM  **LAM1NAR,  90x  BEND** 

720  GOSUB  1380 
730  GOTO  800 

740  REM  **LOSSES  IN  TURBULENT  FLOW** 

750  REM  **TURBULENT,  HORIZONTAL  TUBE** 

760  LENGTH  -  HLENGTH 
770  GOSUB  1410 

780  REM  **TURBULENT,  90x  BEND** 

790  GOSUB  1600 

800  IF  NRE  «<  2300  THEN  930 

810  IF  NRE  >  2300  THEN  820 

820  INLET  -  SAMPEFF 

830  HORIZ  -  TBREFF  *  TINEFF 

840  BEND90  -  TBENDEFF 

850  VERT  •  TBREFF  *  TINEFF 

860  BEND45  «  TBENDEFF 

8'’0  REM**CALCULATION  OF  INCLINED  PIPE  SCALING  FACTOR** 

880  IF  VDOT  «  130  THEN  GOSUB  1650 
890  IF  VDOT  «  70  THEN  GOSUB  1720 
900  INCLINE  «  TBREFF  *  TINEFF  *  FACTOR 

910  TOTAL  •  INLET  *  HORIZ  *  BEND90  *  VERT  *  BEND45  *  INCLINE 

920  GOTO  1020 

930  INLET  »  SAMPEFF 

940  HORIZ  -  BREFF  *  LGREFF 

950  BEND90  -  LBNDEFF 

960  VERT  •  BREFF 

970  BEND45  -  LBNDEFF 

980  GOSUB  1720 

990  INCLINE  «  BREFF  *  LGREFF  *  FACTOR 

1000  TOTAL  -  INLET  *  HORIZ  *  BEND90  *  VERT  *  BEND45  *  INCLINE 
1010  GOTO  1020 

1020  PRINT#2,  USING"##. ##"; DP; 

1030  PRINT#2,  USING"##.### . ; INLET; 

1040  PRINT#2,  USING"##.### . ; HORIZ; 

1050  PRINT#2,  USING"##.### . ;BEND90; 

1060  PRINT#2,  USING"##.### . ;VERT; 

1070  PRINT#2,  USING"##.### . ;BEND45; 

1080  PRINT#2,  USING"##.### . ; INCLINE; 

1090  PRINT#2,  USING"##.### . ;TOTAL 

1100  GOTO  1110 
1110  NEXT  I 
1120  CLOSE  #1 
1130  CLOSE  #2 
1140  END 

1*50  REM  **CALCULATION  OF  BROWNIAN  LOSSES,  XBI  «<  .02** 

1160  A  -  .177  *  (XSI  "  1.3333333333333#) 

1170  B  -  1.2  *  XSI 

1180  C  -  2.56  *  {XSI  “  .66666667#) 

1190  BREFF  -1-C+B+A 
1200  GOTO  680 

1210  REM  **CALCULATION  OF  BROWNIAN  LOSSES,  XSI  >.02** 

1220  A  -  .032  *  EXP(-57  *  XSI) 

1230  B  -  ,097  *  EXP(-22.3  *  XSI) 
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1240  C  -  .819  *  EXP( -3.657  *  XSI ) 

1250  BREFF  -  A  +  B  +  C 
1260  GOTO  680 

1270  REM  **LAMINAR,  GRAVITATIONAL  SEDIMENTATION  LOSSES** 

1280  PHI  -  4.56036  *  LENGTH  *  VS  *  DS  /  VDOT 
1290  IF  PHI  >  1!  THEN  GOTO  1630 

1300  REM**For  explanation  of  line  1010  see  fig.  6  In  Flssan  paper** 

1310  PHIONE  -  PHI" (1/3) 

1320  PHITWO  -  PHI "(2/3) 

1330  A  -  PHIONE  *  SQR( 1  -  PHITWO) 

1340  B  -  ATN( PHIONE  /  SQR( -PHIONE  *  PHIONE  +1)) 

1350  C  «  2  *  PHI  *  SQR(1  -  PHITWO) 

1360  LGREFF  -  1  -  .6366198  *  (C  +  B  -  A) 

1370  RETURN 

1380  REM  **LAMINAR,  90x  BEND  CALCULATIONS** 

1390  LBNDEFF  -  1  -  STK 
1400  RETURN 

1410  REM  **TURBULENT  FLOW  MODELS** 

1420  REM  **CALCULATION  OF  SCHMIDT  NUMBER** 

1430  SC  -  NU  /  RHOAIR  /  DIFF 

1440  REM  **CALCULATION  OF  FANNING  FRICTION  FACTOR** 

1450  FF  «  .316  /  4  /  (NRE‘.25) 

1460  REM  **CALCULATION  OF  PARTICLE  DEPOSITION  VELOCITY,  BROWNIAN  DIFFUSION** 

1470  UBR  -  .0013814677#  *  VDOT  *  (FF".5)  *  ( SC"-. 6666666667# )  /  DS  /  DS 
1480  REM  **CALCULATION  OF  LOSSES  DUE  TO  BROWNIAN  MOTION** 

1490  ARGTBR  -  1459.318  *  DS  *  LENGTH  *  UBR  /  VDOT 
1500  TBREFF  -  EXP( -ARGTBR) 

1510  REM  **TURBULENT  LOSSES  INERTIAL  EFFECTS** 

1520  TAU  -  5.555556E-11  *  RHO  *  DP  *  DP  *  CC  /  NU 

1530  USTAR  -  2.325822E-02  *  (FF".5)  *  VDOT  /  DS  /  DS 

1540  TAUPLUS  -  TAU  *  USTAR  *  USTAR  *  RHOAIR  /  NU 
1550  UPLUS  «  6.000001E-04  *  TAUPLUS  *  TAUPLUS 
1560  UIN  -  UPLUS  *  USTAR 

1570  ARGTIN  ■  1459.318  *  DS  *  LENGTH  *  UIN  /  VDOT 
1580  TINEFF  «  EXP( -ARGTIN) 

1590  RETURN 

1600  REM  **TURBULENT  FLOW  BEND  EFFICIENCY** 

1610  TBENDEFF  -  10" (-.963  *  STK) 

1620  RETURN 
1630  LGREFF  «  01 
1640  GOTO  1370 

1650  REM**INCLINED  PIPE  FACTOR  FOR  130  l/ra±n** 

1660  IF  DP  <  15  GOTO  1680 
1670  IF  DP  ->  15  GOTO  1700 

1680  FACTOR  -  1.000899  -  ( 1. 234803E-04  *  DP)  -  (1.994214E-04  *  DP  *  DP)  +  (7.352 
604E-06  *  DP  *  DP  *  DP) 

1690  RETURN 

1700  FACTOR  -  .8500409  +  (1.249575E-02  *  DP)  -  ( 2. 598943E-04  *  DP  *  DP) 

1710  RETURN 

"20  REM  ** INCLINED  PIPE  FACTOR  FOR  70  1/mln** 

730  FACTOR  -  .9994884  +  (2.10849E-04  *  DP)  -  (3.401449E-04  *  DP  *  DP)  +  (3.4422- 
3E-06  *  DP  *  DP  *  DP) 

1740  RETURN 
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10  REM  **'This  program  calculates  the  mass-weighted  average  efficiency** 

20  REM  **of  a  log-normal  distribution  of  particles  with  equal  ** 

30  REM  **logarithmic  intervals.  This  program  only  works  for  twenty  ** 

40  REM  **intervals.  The  input  to  this  program  is  the  output  from  the** 

50  REM  **"fissan.bas"  program.  This  program  prompts  for  the  input  ** 

60  REM  **and  output  filenames.  ** 

70  DIM  FRACMASS(20) 

80  DIM  BINEFF(20,8) 

90  DIM  AVGEFF(8) 

100  INPUT"Input  Filename  - >",IFILES 

110  INPUT"Output  Filename  - >",OFILE$ 

120  OPEN  "a:"+IFILE$  FOR  INPUT  AS  #1 
130  INPUT#1,  A$ 

140  1NPUT#1,BS 
150  INPUT#1,C$ 

160  INPUT#1,D$ 

170  REM  **L0AD  ARRAY  WITH  EFFICIENCIES  FOR  EACH  BIN** 

180  FOR  I  -  1  TO  20 
190  FOR  J  -  1  TO  8 
200.  INPUT#1,  BINEFF(I,J) 

210  NEXT  J 
220  NEXT  I 
230  CLOSE  #1 

240  REM  **LOADlNG  MASS  FRACTION  ARRAY** 

250  FOR  I  -  1  TO  20 
250  READ  MASSFRAC 
270  FRACMASS(I)  -  MASSFRAC 
280  NEXT  I 

290  REM  **WEIGHTED  AVERAGE  OF  EFFICIENCIES  TO  OUTPUT  FILE** 

300  FOR  J  -  2  TO  8 
310  AVGEFFO  «  01 
320  FOR  I  «  I  TO  20 

330  AVGEFF(J)  -  AVGEFFO  +  (BINEFF(I,J)  *  FRACMASS(I)) 

340  AVGEFFO  -  AVGEFF(J) 

350  NEXT  I 
360  NEXT  J 

370  REM  **OUTPUT  AVERAGE  EFFICIENCIES  TO  OUTPUT  FILE** 

380  OPEN  "a:"+OFILES  FOR  OUTPUT  AS  #1 
390  PRINT#1,A$ 

400  PRINT#1,BS 
410  PRINT#!,"  " 

420  PRINT#1, "MASS-WEIGHTED  FRACTIONAL  PENETRATION  OF  LOG-NORMAL  DISTRIBUTION  OF 
AEROSOLS" 

430  PRINT#1,"  INLET  HORIZ.  90x  BEND  VERTICAL  45x  BEND  INCLINED  OVERAL 
L" 

440  FOR  I  •  2  TO  8 

450  PRINT#1, USING  "##.### . ;  AVGEFF(I); 

460  NEXT  I 
470  CLOSE  #1 

480  DATA  .0005,  .0013,  .0038,  .0093,  .0201,  .0384,  .0644,  .0979,  .1236,  .1406,  . 
1406,  .1236,  .0979,  .0644,  .0384,  .0201,  .0093,  .0038,  .0013,  .0005 
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COMPLETE  SET  OF  NUMERICAL  CALCULATIONS  GENERATED 
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TablA  B-1  Transmission  of  Monodisperse  Aerosol 
Particles  through  Individual  Components  of  the  Aerosol 
Transport  System.  Case  1:  Flow  rate  »  70  L/min,  Inlet 
Orientation  *  90*.  The  column  entitled  "inlet"  includes 
the  aspiration  efficiency  (Vincent  et  al.;  1986)  and 
losses  in  the  developing  boundary  layer  (Okazaki  and 
Willeke,  1987) .  Identification  of  components  is  given 
in  Figure  1. 


Aero¬ 

dynamic 

Particle 


Diameter 

Inlet 

PI 

P2 

P3 

El 

E2 

3  um 

96.21% 

99.741 

100.00% 

99.72% 

99.02% 

99.51% 

5 

90.49 

99.29 

99.99 

99.25 

97.34 

98.66 

7 

83.16 

98.62 

99.98 

98.54 

94.91 

97.42 

10 

70.86 

97.21 

99.92 

97.05 

89.94 

94.89 

15 

51.52 

93.80 

99.62 

93.43 

78.88 

88.82 

20 

36.60 

89.11 

98.80 

88.41 

65.67 

81.04 

25 

26.08 

83.18 

96.91 

81.99 

51.89 

72.04 
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Table  B-2  Transmission  of  Monodisperse  Aerosol 
Particles  through  Individual  Components  of  the  Aerosol 
Transport  System.  Case  2:  Flow  rate  =  130  L/min,  Inlet 
Orientation  =  90*.  The  column  entitled  "inlet"  includes 
the  aspiration  efficiency  (Vincent  et  al.;  1986)  and 
losses  in  the  developing  boundary  layer  (Okazaki  and 
Willeke,  1987) .  Identification  of  components  is  given 
in  Figure  1.  ^ 


Aero¬ 

dynamic 


r'amcie 

diameter 

Inlet 

PI 

P2 

P3 

El 

E2 

3  um 

97.22 

99,86 

99.99 

99.85 

98.18 

99.09 

5 

92.92 

99.60 

99.96 

99.57 

95.12 

97.53 

7 

87.22 

99.21 

99.85 

99.14 

90.75 

95.26 

10 

77.16 

98.30 

99.39 

98.11 

82.13 

90.63 

15 

59.89 

95.55 

96.79 

94.80 

64.37 

80.23 

20 

45.13 

87.97 

84.41 

84.13 

45.79 

67.67 

25 

33.82 

80.24 

73.65 

73.65 

29.57 

54.38 
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Table  B-3  Transmission  of  a  Polydisperse  Aerosol 
through  Individual  Components  of  the  Transport  System. 
Case  3;  Flow  rate  *»  70  L/min,  Inlet  Orientation  =  90*. 
The  column  entitled  "inlet"  includes  the  aspiration 
efficiency  (Vincent  et  al.;  1986)  and  losses  in  the 
developing  boundary  layer  (Okazaki  and  Willeke,  1987) . 


Mass 

Median 


dynamic  - Componant 


Diameter 

Inlet 

PI 

P2 

P3 

El 

E2 

5  um 

88.36% 

99.11% 

99.98% 

99.10% 

96.74% 

98.39% 

15 

51.00 

93.92^ 

99.56 

94.87 

81.14 

90.93 

Table  B-4  Transmission  of  a  Polydisperse  Aerosol 

through  Individual  Components  of  the  Transport  System. 
Case  4;  Flowrate  «  130  L/min,  Inlet  Orientation  *=  90*. 
The  column  entitled  "inlet"  includes  the  aspiration 
efficiency  (Vincent  et  al.;  1986)  and  losses  in  the 
developing  boundary  layer  (Okazaki  and  Willeke),  1987. 

Mass 

Median 


dynamic  - Component 


Diameter 

Inlet 

PI 

P2 

P3 

El 

E2 

5  um 

91.15% 

99.47% 

99.90% 

99.47% 

93.97% 

97.04% 

15 

58.40 

94.35 

96.70 

96.12 

68.17 

85.22 
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Table  B-5  Transmission  of  Monodisperse  Aerosol 
Particles  through  Individual  Components  of  the  Aerosol 
Transport  System.  Case  5:  Flow  rate  «  70  L/min,  Inlet 
Orientation  =  0*.  The  column  entitled  "inlet"  includes 
the  aspiration  efficiency  (Vincent  et  al.;  1986)  and 
losses  in  the  developing  boundary  layer  (Okazaki  and 
Willeke,  1987} .  Identification  of  components  is  given 
in  Figure  1. 


Aero¬ 

dynamic 

Particle 


Diameter 

Inlet 

PI 

P2 

P3 

El 

E2 

3  um 

99.43% 

99.74% 

100.00% 

99.72% 

99.02% 

99.51% 

5 

98.94 

99.29 

99.99 

99.25 

97.34 

98.66 

7 

98.44 

98.62 

99.98 

98.54 

94.91 

97.42 

10 

97.71 

97.21 

99.92 

97.05 

89.94 

94.84 

15 

96.40 

93.80 

99.62 

93.43 

78.88 

88.82 

20 

94.73 

89.11 

98.80 

88.41 

65.67 

81.04 

25 

92.44 

83.18 

96.91 

81.99 

51.89 

72.04 
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Tabl«  B-6  Transmission  of  Monodisperse  Aerosol 
Particles  through  Individual  Components  of  the  Aerosol 
Transport  System.  Case  6:  Flow  rate  «  130  L/min,  Inlet 
Orientation  =  0*.  The  column  entitled  "inlet"  includes 
the  aspiration  efficiency  (Vincent  et  al.;  1986)  and 
losses  in  the  developing  boundary  layer  (Okazaki  and 
Willeke,  1987).  Identification  of  components  is  given 
in  Figure  i. 


Aero¬ 

dynamic 


Diameter 

Inlet 

PI 

P2 

P3 

El 

E2 

3  un 

99.26% 

99.89% 

99.99% 

99.85% 

98.18% 

99.09% 

5 

98.33 

99.60 

99.96 

99.57 

95.12 

97.53 

7 

97.14 

99.21 

99.85 

99.14 

90.75 

95.26 

10 

94.97 

98.30 

99.39 

98.11 

82.13 

90.63 

15 

90.60 

95.55 

96.79 

94.80 

64.37 

80.23 

20 

85.67 

87.97 

84.41 

84.13 

45.79 

67.67 

25 

80.54 

80.24 

73.65 

73.65 

29.57 

54.38 
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Tabl«  B-7  Penetration  of  Monodlsperse  Aerosols  through 
the  Overall  Tubing  System. 


Case 

Case 

Case 

Case 

1 

2 

5 

6 

90* 

Inlet 

Orientation 

0* 

Aerodynamic 

Particle 

Diameter 

Flow  Rate, 

L/min 

70 

130 

70 

130 

3 

94.28% 

94.29% 

97.43% 

96.27% 

5 

85.64 

85.46 

93.63 

90.43 

7 

74.70 

74.04 

88.43 

82.46 

10 

56.98 

55.05 

78.57 

67.76 

15 

31.51 

27.12 

58.96 

41.02 

20 

15.16 

8.74 

39.24 

16.59 

25 

6.44 

2.37 

22.84 

5.64 
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Tabl*  B-8  Transmission  of  a  Polydisperse  Aerosol 
through  Individual  Components  of  the  Transport  System. 
Case  7:  Flowrate  -  70  L/min,  Inlet  Orientation  »  0*. 
The  column  entitled  -  "inlet"  includes  the  aspiration 
efficiency  (Vincent  et  al.;  1986)  and  losses  in  the 
developing  boundary  layer  (Okazaki  and  Willeke,  1987) . 


Mass 
Median 
Aero¬ 
dynamic  — - 

Diameter  Inlet _ PI 


Component 
P2  P3 


El  E2 


5  urn  98.83%  99.03%  99.98%  99.03%  96.49%  98.26% 

15  95.58  91.82  99.26  93.71  76.47  88.91 


Table  B-9  Transmission  of  a  Polydisperse  Aerosol 
through  Individual  Components  of  the  Transport  System. 
Case  8:  Flowrate  *  130  L/min,  Inlet  Orientation  =  0*. 
The  column  entitled  "inlet"  includes  the  aspiration 
efficiency  (Vincent  et  al.;  1986)  and  losses  in  the 
developing  boundary  layer  (Okazaki  and  Willeke) ,  1987). 


Mass 

Median 


dynamic  - Component 


Diameter 

Inlet 

PI 

P2 

P3 

El 

E2 

5  um 

97.96% 

99.44% 

99.88% 

99.44%  - 

93.67% 

96.90% 

15 

89.10 

92.58 

95.65 

95.42 

63.78 

83.57 
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Table  B-lO  Penetration  of  Polydisperse  Aerosols  through 
the  Overall  Tubing  System.  The  aerosol  in  the  free 
stream  is  assumed  to  be  log-normally  distributed  with  a 
geometric  standard  deviation  of  1.5. 


Case  Case  Case  Case 

_3 _ £ _ 7 _ 8 

Inlet  Orientation 
«  • 

90 _  _ 0 _ 

_ Flow  Rate,  L/min _ 

70 _ no _ 70 _ 130 

5  um  82.58%  82.15%  91.88%  87.82% 

15  33.38  29.75  55.51  40.13 


Mass 

Median 

Aerodynamic 

Particle 

Diameter 
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APPENDIX  C 

COMPLETE  SET  OF  NUMERICAL  CALCULATIONS  GENERATED 

BY  IITRI  MODEL 
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rabl«  C-1  Sampling  Efficiency  as  a  Function  of  Particle  Size 
(Alpha  =  0  degrees /Velocity  Ratio  =  1.3) 


0 

1  DIAMETER,  2 
micron 

STOKES  t 

3  CONC. 
RATIO 

1 

3 

0.0058 

0.940611 

2 

5 

0.0162 

0.940611 

3 

7 

0.0318 

0.940611 

4 

10 

0.0649 

0.930008 

5 

15 

0.1459 

0.919466 

6 

20 

0.2595 

0.919466 

7 

25 

0.4054 

0.908984 

fable  c-2  Sampling  Efficiency  as  a  Function  of  Particle  Size 
(Alpha  =0/  Velocity  Ratio  *  0.7) 


0 

1  DIAMETER  2 
micron 

STOKES  # 

3  CONC 
RATIO 

1 

3 

0.0058 

1.007484 

2 

5 

0.0162 

1.007484 

3 

7 

0.0318 

1.007484 

4 

10 

0.0649 

1.007484 

5 

15 

0.1459 

1.083860 

6 

20 

0.2595 

1.146970 

7 

25 

0.4054 

1.211865 
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Table  c-3  Sampling  Efficiency  as  a  Function  of  Particle  size 
(Alpha  =  90  degrees;  Velocity  Ratio  =  1.3) 


0 

1  DIAMETER  2 
micron 

STOKES  # 

3  CONC 
RATIO 

1 

3 

0.0058 

0.854107 

2 

5 

0.0162 

0.852980 

3 

7 

0.0318 

0.847058 

4 

10 

0.0649 

0.848458 

5 

15 

0.1459 

0.811773 

6 

20 

0.2595 

0.760660 

7 

25 

0.4054 

0.701941 

fable  C-4  Sampling  Efficiency  as  a  Function  of  Particle  Size 
(Alpha  =  90  degrees;  Velocity  Ratio  =  0.7) 


0 

1  DIAMETER  2 
micron 

STOKES  t 

3  CONC 
RATIO 

1 

3 

0.0058 

0.892986 

2 

5 

0.0162 

0.893216 

3 

7 

0.0318 

0.883381 

4 

10 

0.0649 

0.871970 

5 

15 

0.1459 

0.799119 

6 

20 

0.2595 

0.759128 

7 

25 

0.4054 

0.713426 

APPENDIX  D 

COMPLETE  SET  OF  NUMERICAL  CALCULATIONS 
GENERATED  BY  UNIVERSITY  OF  DUISBURG  MODEL 
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Table  D-1.  Penetration  Efficiency  Calculations  For  Monodi spersed  Particle 
Distributions  -  0*  Inlet  Orientation  &  70  L/min  Flow  Rate 


Monodi spersed  Aerosol  Size  [=]  pm 

Volumetric  Flow  Rate  =  70  L/min,  Wind  Speed  =  3  m/sec 


SIZE 

3.00 

5.00 

7.00 

10.00 

15.00 

20.00 

25.00 


INLET 
1 .007E+00 
1 .018E+00 
1 ,033E+00 
1 .062E+00 
1 . 1 19E+00 
1 . 174E+00 
1 .222E+00 


FRACTIONAL  PENETRATION  OF  AEROSOL 


HORIZ. 

1 .OOOE+00 
9.999E-01 
9.99aE-01 
9.993E-01 
9.963E-01 
9.884E-01 
9.721E-01 


90*  BEND 
9.901E-01 
9. 734E-01 
9.490E-01 
8.994E-01 
7.888E-01 
6.567E-01 
5. 189E-01 


VERTICAL 

1 .OOOE+00 

9.999E-01 

9.998E-01 

9.993E-01 

9.963E-01 

9.884E-01 

9.721E-01 


45*  BEND 
9. 90  IE-01 
9.734E-01 
9.490E-01 
8.994E-01 
7.88aE-01 
6.567E-01 
5. 189E-01 


INCLINED 

9.971E-01 

9.924E-01 

9.853E-01 

9.703E-01 

9.343E-01 

8.848E-01 

8.224E-01 


CUM. 

9.841E-01 
9.569E-01 
9. 165E-01 
8.327E-01 
6.455E-01 
4.376E-01 
2.556E-01 


Table  D-2.  Penetration  Efficiency  Calculations  For  Monodi spersed  Particle 
Distributions  -  0*  Inlet  Orientation  &  130  L/min  Flow  Rate 

Monodispersed  Aerosol  Size  [  =  ]  pfn 

Volumetric  Flow  Rate  =  130  L/min,  Wind  Speed  =  3  m/sec 


FRACTIONAL  PENETRATION  OF  AEROSOL 


SIZE 

INLET 

HORIZ. 

90*  BEND 

VERTICAL 

45*  BEND 

INCLINED 

CUM. 

3.00 

9.960E-01 

9.999E-01 

9.818E-01 

9.999E-01 

9.818E-01 

9.989E-01 

9.588E-01 

5.00 

9.894E-01 

9.996E-01 

9.512E-01 

9.996E-01 

9.512E-01 

9.958E-01 

8.908E-01 

7.00 

9.804E-01 

9.985E-01 

9.074E-01 

9.985E-01 

9.074E-01 

9.913E-01 

7 .979E-01 

10.00 

9.635E-01 

9.940E-01 

8.213E-01 

9.940E-01 

8.213E-01 

9.812E-01 

6.301E-01 

15.00 

9.318E-01 

9.706E-01 

6.437E-01 

9.706E-01 

6.437E-01 

9.502E-01 

3.456E-01 

20.00 

9.020E-01 

9. 104E-01 

4.579E-01 

9. 104E-01 

4.579E-01 

9.068E-01 

1  .421E-01 

25.00 

8.769E-01 

7.958E-01 

2.957E-01 

7.958E-01 

2.957E-01 

7.958E-01 

3.864E-02 
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Table  D-3.  Penetration  Efficiency  Calculations  for  Pol ydi spersed  Particle 
Distributions  -  MMD  =  ISum.  Inlet  Orientation  =  0‘, 

Flow  Rate  =  70  L/min 


MMD  =  15  Mm  o  =  1.5  Mm  Size  [-]  Mm 

Volumetric  Flow  Rate  =  70  l/min,  Wind  Speed  =  3  m/sec 

FRACTIONAL  PENETRATION  OF  AEROSOL 

SIZE  INLET  HORI2.  90*  BEND  VERTICAL  45*  BEND  INCLINED  CUM. 

3.72  1.010E+00  1.000E+00  9.850E-01  1.000E+00  9.850E-01  9.957E-01  9.758E-01 

4.31  1.013E+00  1.000E+00  9.801E-01  1.000E+00  9.801E-01  9.943E-01  9.678E-01 

4.99  1.018E+00  9.999E-01  9.735E-01  9.999E-01  9.735E-01  9.924E-01  9.571E-01 

5.78  1.023E+00  9.999E-01  9.648E-01  9.999E-01  9.648E-01  9.899E-01  9.428E-01 

6.69  1.031E+00  9.998E-01  9.533E-01  9.998E-01  9.533E-01  9.865E-01  9.236E-01 

7.75  1.040Er00  9.997E-01  9.380E-01  9.997E-01  9.380E-01  9.820E-01  8.981E-01 

8.97  1.052E+00  9.995E-01  9.180E-01  9.995E-01  9.180E-01  9.760E-01  8.643E-01 

10.39  1.067E+00  9.991E-01  8.919E-01  9.991E-01  8.919E-01  9.680E-01  8.198E-01 

12.04  1.085E+00  9.985E-01  8.580E-01  9.985E-01  8.580E-01  9.573E-01  7.621E-01 

13.94  1.107E+00  9.972E-01  8.146E-01  9.972E-01  8.146E-01  9.430E-01  6.887E-01 

16.14  1.132E+00  9.951E-01  7.599E-01  9.951E-01  7.599E-01  9.241E-01  5.981E-01 

18.69  1.160E+00  9.912E-01  6.923E-01  9.912E-01  6.923E-01  8.990E-01  4.911E-01 

21.65  1.191E+00  9.842E-01  6.111E-01  9.842E-01  6.111E-01  8.656E-01  3.728E-01 

25.07  1.222E+00  9.718E-01  5.169E-01  9.718E-01  5.169E-01  8.214E-01  2.534E-01 

29.04  1.254E+00  9.499E-01  4.130E-01  9.499E~01  4.130E-01  7.629E-01  1.473E-01 

33.63  1.284E+00  9.119E-01  3.057E-01  9.119E-01  3.057E-01  6.865E-01  6.853E-02 

38.94  1.312E-f00  8.473E~01  2.043E-01  8.473E-01  2.043E-01  5.890E-01  2.315E-02 

45.10  1.336E+00  7.425E-01  1.190E-01  7.425E-01  1.190E-01  4.700E-01  4,900E-03 

52.23  1.357E+00  5.857E~01  5.762E-02  5.857E-01  5.762E-02  3.356E-01  5.187E-04 

60.49  1.374E+00  3.824E-01  2.180E-02  3.824E-01  2.180E-02  2.025E-01  1.933E-05 


MMD  =  15  Mm  o  =  1.5  Mm  Size  [=]  Mm 

Volumetric  Flow  Rate  =  70  l/min.  Wind  Speed  =  3  m/sec 

MASS-WEIGHTED  FRACTIONAL  PENETRATION  OF  LOG-NORMAL  DISTRIBUTIOf  OF  AEROSOLS 
INLET  HORIZ.  90*  BEND  VERTICAL  45*  BEND  INCLINED  CUM. 

1.129E+00  9.870E-01  7.466E-01  9.870E-01  7,466E-01  9.104E-01  5.990E-01 
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Table  D-4,  Penetration  Efficiency  Calculations  For  Pol ydi spersed  Particle 
Distributions  -  MMD  =  15  pm,  Inlet  Orientation  =  O', 

Flow  Rate  =  130  L/min 


MMD  =  16  nm  o  =  1.5  pm  Size  C  =  1  Ui" 

Volumetric  Flow  Rate  =  130  l/min.  Wind  Speed  =  3  m/sec 

FRACTIONAL  PENETRATION  OF  AEROSOL 

SIZE  INLET  HORIZ.  90*  BEND  VERTICAL  45*  BEND  INCLINED  CUM. 
3.72  9.940E-01  9.999E-01  9.724E-01  9.999E-01  9.724E-01  9.979E-01  9.376E-01 

4.31  9.920E-01  9.998E-01  9.634E--01  9.998E-01  9.634E-01  9.970E-01  9.175E-01 

4.99  9.895E-01  9.996E-01  9.514E-01  9.996E-01  9.514E-01  9.958E-01  8.912E-01 

5.78  9.862E-01  9.993E-01  9.357E-01  9  993E-01  9.357E-01  9.943E-01  8.572E-01 

6.69  9.819E-01  9.988E-'01  9.150E-01  9.988E-01  9.150E-01  9.921E-01  8.136E-01 

7.75  9.764E-01  9.978E-01  8.880E-01  9.978E-01  8.880E-01  9.892E-01  7.584E-01 

8.97  9.696E-01  9.961E-01  8.531E-01  9.961E-01  8.531E-01  9.852E-01  6.898E-01 

10.39  9.610E-01  9.931E~01  8.085E-01  9.931E-01  8.085E-01  9.795E-01  6.069E-01 

12.04  9.507E-01  9.876E-01  7.524E-01  9.876E-01  7.524E-01  9.712E-01  5.099E-01 

13.94  9.386E-01  9.779E-01  6.833E-01  9.779E-01  6.833E-01  9.587E-01  4.018E-01 

16.14  9.247E-01  9.608E-01  6.006E-01  9.608E-01  6.006E-01  9.45^r-01  2.911E-01 

18.69  9.093E-01  9.308E-01  5.052E-01  9.308E-01  5.052E-01  9.241E-01  1.858E-01 

21.65  8.931E-01  8.792E-01  4.006E-01  8.792E-01  4.006E-01  8.781E-01  9.730E-02 

25.07  8,766E-01  7.937E-01  2.936E-01  7,937E-01  2.936E-01  7.937E-01  3.778E-02 

29.04  8.605E-01  6.604E-01  1.936E-01  6.604E~01  1.936E-01  6.563E-01  9.227E-03 

33.63  8.454E-01  4.747E-01  1 . 107E-01  4.747E-01  1.107E-01  4.634E-01  1.082E-03 

38.94  8.318E-01  2.622E-01  5.23SE-02  2.622E-01  5.235E-02  2.472E-01  3.875E-05 

45.10  8.201E-01  9.027E-02  1.918E~02  9.027E-02  1.918E-02  7.989E-C2  1.964E-07 

52.23  8.102E-01  1 . 328E-02  4.991E-03  1.328E-02  4.991E-03  1.054E-02  3.749E-11 

60.49  8.020E-01  4.235E-04  8.206E-04  4.235E-04  8.206E-04  2.774E-04  2.686E-17 


MMD  =  15  pm  o  =  1.5  pm  Size  [=]  pm 

Volumetric  Flow  Rate  =  130  l/min,  Wind  Speed  =  3  m/sec 

MASS-WEIGHTED  FRACTIONAL  PENETRATION  OF  LOG-NORMAL  DISTRIBUTION  OF  AEROSOLS 
INLET  HORIZ.  90’  BEND  VERTICAL  45*  BEND  INCLINED  CUM. 
9.266E-01  9.192E-01  6.049E-01  9.192E-01  6.049E-01  9.083E-01  3.581E-01 
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Table  D-5.  Penetration  Efficiency  Calculations  For  Pol ydi spersed  Particle 
Distributions  -  MMD  =  5  urn,  Inlet  Orientation  =  O’, 

Flow  Rate  =  70  L/min 


MMD  =  5  pm  o  =  1 , 5  Mm  Size  [=]  Mm 

Volumetric  Flow  Rate  =  70  l/min,  Wind  Speed  =  3  m/sec 


FRACTIONAL  PENETRATION  OF  AEROSOL 


SIZE 

INLET 

HORIZ. 

90’  BEND 

VERTICAL 

45*  BEND 

INCLINED 

CUM. 

1 .24 

1  .OCE+OO 

1 .OOOE+00 

9.982E-01 

1  .OOOE+00 

9.982E-01 

9.992E-01 

9.967E-01 

1 .44 

1  .'  D2E-^00 

1 .OOOE+00 

9.976E-01 

1  .OOOE+00 

9.976E-01 

9.991E-01 

9.958E-01 

1.66 

1 .002E-t-00 

1 .OOOE+00 

9,968E-01 

1  .OOOE-t-00 

9. 968E-01 

9.989E-01 

9.946E-01 

1 .93 

1 .003E+00 

1 .OOOE+00 

9.958E-01 

1  .OOOE+00 

9.958E-01 

9.986E-01 

9.930E-01 

2.23 

1 .004E-^00 

1 .000E-^00 

9.944E-01 

1  .OOOE+00 

9.944E-01 

9.983E-01 

9.909E-01 

2.58 

1 .005E+00 

1 .OOOE+00 

9.926E-01 

1 .OOOE+OO 

9 , 926E-01 

9 . 978E-01 

9.880E-01 

2.99 

1 .007E+00 

1 .OOOE+00 

9.902E-01 

1 .OOOE+OO 

9.902E-01 

9.971E-01 

9.841E-01 

3  46 

1 .009E-f00 

1 .OOOE+00 

9.870E-01 

1 .OOOE+OO 

9.870E-01 

9.962E-01 

9. 790E-01 

4.01 

1  .012E+00 

1 .OOOE+00 

9.827E-01 

1 .OOOE+OO 

9.827E-01 

9.950E-01 

9.720E-01 

4.65 

1  .015E+00 

9.999E-01 

9.770E-01 

9.999E-01 

9. 770E-01 

9.934E-01 

9.627E-01 

5.38 

1  .020E+00 

9.999E-01 

9.694E-01 

9.999E-01 

9.694E-01 

9.912E-01 

9.503E-01 

6.23 

1 .027E+00 

9.999E-01 

9.593E-01 

9.999E-01 

9.593E-01 

9.883E-01 

9. 336E-01 

7,22 

1 .035E+00 

9.998E-01 

9.460E-01 

9.998E-01 

9.460E-01 

9.844E-01 

9. 1  14E-01 

8.36 

1  .046E+00 

9.996E-01 

9.284E-01 

9.996E-01 

9.284E-01 

9.791E-01 

8.819E-01 

9.68 

1 .059E-^00 

9.993E-01 

9.054E-01 

9.993E-01 

9.054E-01 

9.721E-01 

8.429E-0  1 

11.21 

1 .076E+00 

9.988E-01 

8.755E-01 

9.988E-01 

8. 755E-01 

9.628E-01 

7.919E-CI 

12.98 

1 .096E-f00 

9,979E-01 

8,370E-01 

9.979E-01 

8.370E-01 

9.504E-01 

7.264E-t;» 

15.03 

1  .  1 19E+00 

9.963E-01 

7.880E-01 

9.963E-01 

7 . i  iOE-01 

9.340E-01 

6 .442E-01 

17,41 

1  .  146E-»-00 

9.933E-01 

7 .268E-01 

9.933E-01 

7 .268E-01 

9. 121E-01 

5.448E-01 

20.16 

1  .  176E-^00 

9.881E-01 

6.522E-01 

9.881E-01 

6.522E-01 

8,830E-01 

4.311E-01 

MMD  =  5  Mm  o  =  1 . 5  Mm  Size  (=]  Mm 

Volumetric  Flow  Rate  =  70  l/min,  Wind  Speed  =  3  m/sec 

MASS -WEIGHTED  FRACTIONAL  PENETRATION  OF  LOG-NORMAL  DISTRIBUTION  OF  AEROSOLS 
INLET  HORIZ.  90*  BEND  VERTICAL  45*  BEND  INCLINED  CUM. 

1.023E+00  9.996E-01  9.639E-01  9.996E-01  9.639E-01  9.893E-01  9.408E-01 
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Table  D-6.  Penetration  Efficiency  Calculations  For  Pol ydispersed  Particle 
Distributions  -  MMD  =  5  pm.  Inlet  Orientation  =  0*, 

Flow  Rate  =  130  L/min 


MMD  =  5  pm  o  =  1.5  pm  Size  t=]  pm 

Volumetric  Flow  Rate  =  130  l/min,  Wind  Speed  =  3  m/sec 


FRACTIONAL  PENETRATION  OF  AEROSOL 


SIZE 

INLET 

HORIZ. 

90*  BEND 

VERTICAL 

45*  BEND 

INCLINED 

CUM. 

1.24 

9.993E-01 

1  .OOOE+00 

9.966E-01 

1  .OOOE+00 

9.966E-01 

1 .OOOE+OO 

9.928E-01 

1.44 

9.990E-01 

1 .OOOE+00 

9.956E-01 

1  .OOOE+00 

9.956E-01 

1 .OOOE+OO 

9.904E-01 

1.66 

9.987E-01 

1 .OOOE+00 

9.941E-01 

1 .OOOE+00 

9.941E-01 

1 .OOOE+OO 

9.871E-01 

1 .93 

9.983E-01 

1 .OOOE+00 

9.922E-01 

1 .OOOE+OO 

9.922E-01 

9.999E-01 

9.827E-01 

2.23 

9.977E-01 

1 .OOOE+00 

9.897E-01 

1  .OOOE+OO 

9.897E-01 

9.997E-01 

9.769E-01 

2.58 

9.970E-01 

9.999E-01 

9.863E-01 

9.999E-01 

9.863E-01 

9.993E-01 

9.692E-01 

2.99 

9.960E-01 

9.999E-01 

9.819E-01 

9.999E-01 

9.819E-01 

9.989E-01 

9.590E-01 

3.46 

9.947E-01 

9.999E-01 

9.759E-01 

9.999E-01 

9.759E-01 

9.983E-01 

9.456E-01 

4.01 

9.930E-01 

9.998E-01 

9.681E-01 

9.998E-01 

9.681E-01 

9.975E-01 

9.279E-01 

4.65 

9.908E-01 

9.997E-01 

9.576E-01 

9.997E-01 

9.576E-01 

9.965E-01 

9.048E-01 

5.38 

9.879E-01 

9.995E-01 

9.438E-01 

9.995E-01 

9.438E-01 

9.951E-01 

8. 748E-01 

6.23 

9.841E-01 

9.991E-01 

9.257E-01 

9.991E-01 

9.257E-01 

9.932E-01 

8.360E-01 

7.22 

9.793E-01 

9.984E-01 

9.020E-01 

9.984E-01 

9.020E-01 

9.907E-01 

7.867E-01 

8.36 

9.731E-01 

9.971E-01 

8.71  IE-01 

9.971E-01 

8.71  IE-01 

9.873E-01 

7.248E-01 

9.68 

9.654E-01 

9.948E-01 

8.315E-01 

9.948E-01 

8.315E-01 

9.825E-01 

6.489E-01 

1 1.21 

9.560E-01 

9.907E-01 

7.812E-01 

9.907E-01 

7.812E-01 

9.756E-01 

5.586E-01 

12.98 

9.447E-01 

9.833E-01 

7. 185E-01 

9.833E-01 

7. 185E-01 

9.664E-01 

4.553E-01 

15.03 

9.316E-01 

9.703E-01 

6.424E-01 

9.703E-01 

6.424E-01 

9.501E-01 

3.439E-01 

17.41 

9. 169E-01 

9.474E-01 

5.529E-01 

9.474E-01 

5.529E-01 

9.368E-01 

2.357E-01 

20.  16 

9.01  IE-01 

9.076E-01 

4.521E-01 

9.076E-01 

4.521E-01 

9.043E-01 

1 .372E-01 

MMD  =  5  pm  o  =  1 . 5  pm  S i ze  [ = ]  pm 

Volumetric  Flow  Rate  =  130  l/min.  Wind  Speed  =  130  L/min 

MASS-WEIGHTED  FRACTIONAL  PENETRATION  OF  LOG-NORMAL  DISTRIBUTION  OF  AEROSOLS 
INLET  HORIZ.  90*  BEND  VERTICAL  45*  BEND  INCLINED  CUM. 
9.861E-01  9.984E-01  9.349E-01  9.984E-01  9.349E-01  9.937E-01  8.598E-01 
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APPENDIX  E 

COMPLETE  SET  OP  NUMERICAL  CALCULATIONS  GENERATED 

BY  BATTELLE  MODEL 
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Table  E-1  RESULTS  OF  SAMPLF  CALCULATIONS  USING  HMD  =  5  Atm 
AND  Q  =  70  L/tnin 


Inlet 

Train 

Section 

Local 

Deposition 
Rate, 
mg/mi  n 

Exit 

Concentration, 

ing/in3 

Cumulative 

Sampling 

Efficiency, 

%  (mass) 

Horizontal  pipe 

0.208E-5 

0.197E-2 

98.52 

90®  bend 

0.362E-6 

0.197E-2 

98.26 

Vertical  pipe 

0.272E-8 

0.197E-2 

98.26 

45®  bend 

0.181E-6 

0.196E-2 

98.13 

Inclined  pipe 

0.143E-5 

0.194E-2 

97.11 

Table  E-2  RESULTS  OF  SAMPLF  CALCULATIONS  USING  MMD  =  15 
AND  Q  »  70  L/min 


Inlet 

Train 

Section 

Local 

Deposition 
Rate, 
mg/mi n 

Exit 

Concentration, 

mg/m3 

Cumulative 

Sampling 

Efficiency, 

%  (mass) 

Horizontal  pipe 

0.168E-4 

0.176E-2 

87.99 

90®  bend 

0.136E-5 

0.174E-2 

87.02 

Vertical  pipe 

0.122E-8 

0.174E-2 

87.02 

45®  bend 

0.6S9E-6 

0.173E-2 

86.55 

Inclined  pipe 

0.950E-5 

0.160E-2 

79.76 
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Table  E-3 

RESULTS  OF  SAMPLF 
AND  Q  =  130  L/min 

CALCULATIONS  USING  FWD  =  5  /;ni 

.  Inlet 

Train 

Section 

Local 

Deposition 

Rate, 

mg/min 

Exit 

Concentration, 

mg/m3 

Cumulative 

Sampling 

Efficiency, 

%  mass 

Horizontal  pipe 

0.209E-5 

0.198E-2 

99.19 

90®  bend 

0.212E-5 

0.197E-2 

98.38 

Vertical  pipe 

0.482E-8 

0.197E-2 

98.38 

45®  bend 

0.105E-5 

0.196E-2 

97.97 

Inclined  pipe 

0.143E-5 

0.195E-2 

97.42 

Table  E-4 

RESULTS  OF  SAMPLF  CALCULATIONS  USING  MMD 
AND  Q  =  130  L/min 

*  15  /jm 

Inlet 

Train 

Section 

Local 

Deposition 

Rate, 

mg/min 

Exit 

Concentration, 

mg/m3 

Cumulative 

Sampling 

Efficiency, 

V  (mass) 

Horizontal  pipe 

0,177E-4 

0.186E-2 

93.18 

90®  bend 

0.158E-4 

0.174E-2 

87.11 

Vertical  pipe 

0.220E-8 

0.174E-2 

87.11 

45®  bend 

0.725E-5 

0.169E-2 

84.32 

Inclined  pipe 

0.956E-5 

0a61E-2 

80.65 
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Table  E-5 

RESULTS  OF  SAMPLF  CALCULATIONS  USING  D 
AND  Q  =  70  L/min 

-  3  ftm 

Inlet 

Train 

Section 

Local 

Deposition 

Rate, 

mg/min 

Exit 

Concentration, 

mg/m3 

Cumulative 

Sampling 

Efficiency, 

%  mass 

• 

Horizontal  pipe 

0.548E-6 

0.199E-2 

99.61 

90®  bend 

0.315E-6 

0.199E-2 

99.38 

Vertical  pipe 

0.369E-8 

0.199E-2 

99.38 

45®  bend 

0.157E-6 

0.199E-2 

99.27 

Inclined  pipe 

0.386E-6 

0.198E-2 

98.99 

Table  E-6 

RESULTS  OF  SAMPLF  CALCULATIONS  USING  D 
AND  Q  «  70  L/min 

-  5  fim 

Inlet 

Train 

Section 

Local 

Deposition 

Rate, 

ng/min 

Exit 

Concentration, 

nig/m3 

Cumulative 

Sampling 

Efficiency, 

%  (mass) 

Horizontal  pipe 

0.151E-5 

0.198E-2 

98.92 

90®  bend 

0.342E-6 

0.197E-2 

98.68 

Vertical  pipe 

0.263E-8 

0.197E-2 

98.68 

45®  bend 

0.171E-6 

0.197E-2 

98.55 

Inclined  pipe 

0.106E-5 

0.I96E-2 

97.80 
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Table  E-7  RESULTS  OF  SAMPLF  CALCULATIONS  USING  0  =  7  /tm 
AND  Q  =  70  L/min 


Inlet 

Train 

Section 

Local 

Deposition 

Rate, 

mg/min 

Exit 

Concentration, 

mg/m3 

Cumulative 
Sampl ing 
Efficiency, 

%  mass 

Horizontal  pipe 

0.294E-5 

D.196E-2 

97.90 

90®  bend 

0.385E-6 

0.195E-2 

97.62 

Vertical  pipe 

0.209E-8 

0.195E-2 

97.62 

45®  bend 

0.192E-6 

0.195E-2 

97.48 

Inclined  pipe 

0.204E-5 

0.192E-2 

96.03 

Table  E-8  RESULTS  OF  SAMPLF  CALCULATIONS  USING  0  =  10  /tm 
AND  Q  =  70  L/min 


Inlet 

Train 

Section 

Local 

Deposition 

Rate, 

mg/min 

Exit 

Concentration, 

mg/m3 

Cumulative 
Sampling 
Efficiency, 
k  (mass) 

Horizontal  pipe 

0.594E-5 

0.192E-2 

95.76 

90®  bend 

0.486E>6 

0.191E-2 

95.41 

Vertical  pipe 

0.163E-8 

0.191E-2 

95.41 

45®  bend 

0.242E-6 

0.191E-2 

95.24 

Inclined  pipe 

0.402E>5 

0.185E-2 

92.37 

Table  E-9 

RESULTS  OF  SAMPLF  CALCULATIONS  USING  D  *  15  ;<m 

AND  Q  =  70  L/mln 

i 

Local  Cu-ilative 

Deposition  Exit  Sampling 

Rate,  Concentration,  Efficiency, 

rog/mln  ing/m3  h  (mass) 


Inlet 

Train 

Section 
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Table  E-11 

RESULTS  OF  SAMPLF 
AND  0  =  70  L/min 

CALCULATIONS  USING  0 

=  25 

Local 

Cumulative 

Inlet 

Deposition 

Exit 

Sampl ing 

Train 

Rate, 

Concentration, 

Efficiency, 

Section 

mg/min 

mg/m3 

%  (mass) 

Horizontal  pipe 

0.332E-4 

0.153E-2 

76.29 

90®  bend 

0.265E-5 

0.149E-2 

74.40 

Vertical  pipe 

0.709E-9 

0.149E-2 

74.40 

45®  bend 

0.130E-5 

0.147E-2 

73.47 

Inclined  pipe 

0.179E-4 

0.121E-2 

60.67 

Table  E-12 

RESULTS  OF  SAMPLF  CALCULATIONS  USING  D 

=  3  ^ 

AND  Q  =  130  L/min 

Local 

Cumulative 

Inlet 

Deposition 

Exit 

Sampling 

Train 

Rate, 

Concentration, 

Efficiency, 

Section 

mg/fliin 

mg/m3 

%  mass 

Horizontal  pipe 

0.552E-6 

0.200E-2 

99.79 

90®  bend 

0.148E-5 

0.198E-2 

99.22 

Vertical  pipe 

0.646E-8 

0.198E-2 

99.22 

45®  bend 

0.736E-6 

0.198E-2 

98.93 

Inclined  pipe 

0.388E-6 

0.198E-2 

98.78 
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Table  E-13  results  OF  SAMPLF  CALCULATIONS  USING  D  =  5  /iHi 
AND  Q  =  130  L/min 


Inlet 

Train 

Section 

Local 

Deposition 

Rate, 

mg/min 

Exit 

Concentration, 

ing/in3 

Cumulative 

Sampling 

Efficiency, 

V  (mass) 

Horizontal  pipe 

0.152E-5 

0.199E-2 

99.42 

90®  bend 

0.178E-5 

0.198E-2 

98.73 

Vertical  pipe 

0.466E-8 

0.198E-2 

98.73 

45®  bend 

0.886E-6 

0.197E-2 

98.39 

Inclined  pipe 

0.106E-5 

0.196E-2 

97.98 

Table  E-14 

RESULTS  OF  SAMPLF  CALCULATIONS  USING  D 
AND  Q  =  130  L/min 

*  7  ;tm 

Inlet 

Train 

Section 

Local 

Deposition 

Rate, 

mg/min 

Exit 

Concentration, 

mg/m3 

Cumulative 

Sampling 

Efficiency, 

V  (mass) 

Horizontal  pipe 

0.296E-5 

0.198E-2 

98.86 

90®  bend 

0.231E-5 

0.196E-2 

97.97 

Vertical  pipe 

0.375E-8 

0.196E-2 

97.97 

45®  bend 

0.115E-5 

0.195E-2 

97.53 

Inclined  pipe 

0.204E-5 

0.193E-2 

96.74 
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Table  E-15  RESULTS  OF  SAMPLF  CALCULATIONS  USING  0  =  10  ^ 
AND  Q  =  130  L/min 


Inlet 

Train 

Section 

Local 

Deposition 

Rate, 

mg/min 

Exit 

Concentration, 

ing/m3 

Cumulative 
Sampling 
Efficiency, 
k  mass 

Horizontal  pipe 

0.600E-5 

0.195E-2 

97.69 

90®  bend 

0.388E-5 

0,192E-2 

96.20 

Vertical  pipe 

0.294E-8 

0.192E-2 

96.20 

45®  bend 

0.192E-5 

0.191E-2 

95.46 

Inclined  pipe 

0.406E-5 

0.188E-2 

93.90 

Table  E-16 

RESULTS  OF  SAMPLF  CALCULATIONS  USING  D 
AND  Q  =  130  L/min 

=  15  ;tm 

Inlet 

Train 

Section 

Local 

Deposition 

Rate, 

ng/min 

Exit 

Concentration, 

fflg/m3 

Cumulative 

Sampling 

Efficiency, 

%  (mass) 

Horizontal  pipe 

0.133E-4 

0.190E-2 

94.89 

90®  bend 

0.142E-4 

0.179E-2 

89.43 

Vertical  pipe 

0.213E-8 

0.179E-2 

89.43 

45®  bend 

0.678E-5 

0.174E-2 

86.82 

Inclined  pipe 

0.822E-S 

0.167E-2 

83.66 
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Table  E-17  RESULTS  OF  SAMPLF  CALCULATIONS  USING  D  =  20 
AND  Q  =  130  L/min 


Inlet 

Train 

Section 

Local 

Deposition 
Rate, 
ng/nii  n 

Exit 

Concentration, 

tng/fn3 

Cumulative 

Sampling 

Efficiency, 

%  mass 

Horizontal  pipe 

0.232E-4 

0.182E-2 

91.09 

90®  bend 

0.258E-4 

0.162E-2 

81.17 

Vertical  pipe 

0.162E-8 

0.I62E-2 

81.17 

45®  bend 

0.118E-4 

0.153E-2 

76.62 

Inclined  pipe 

0.127E-4 

0.143E-2 

71.73 

Table  E-18 

RESULTS  OF  SAMPLF  CALCULATIONS  USING  D 
AND  Q  *  130  L/min 

=  25  fun 

Inlet 

Train 

Section 

Local 

Deposition 

Rate, 

mg/min 

Exit 

Concentration, 

fflg/m3 

Cumulative 

Sampling 

Efficiency, 

%  (mass) 

Horizontal  pipe 

0.353E-4 

0.173E-2 

86.44 

90®  bend 

0.357E-4 

0.145E-2 

72.70 

Vertical  pipe 

0.127E>8 

0.145E-2 

72.70 

45®  bend 

0.157E-4 

0.133E-2 

66.68 

Inclined  pipe 

0.170E-4 

0.120E-2 

60.15 
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APPENDIX  F 

COMPLETE  SET  OF  EXPERIMENTAL  DATA  GENERATED  BY 
WIND  TONNEL  TESTS 
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Tabl*  Exp«rlsi«ntal  Valuta  of  Aarosol  Ptnatration. 

Taat  Cast  1.  Wind  spead  ^  3  m/u,  Inlet  Oriented 

Perpendicular  to  Air  Stream  (90*)/  Sampled  Flow  Rate  = 
70  L/nin. 

A£0  Mean  Penetration 

Particle  ±  l  standard 

Size/  Deviation/ 


un 

Penetration  Values/  Percent 

Percent 

3.5 

94.6/  92.4/  92.2 

93.1  ±  1.3 

4.9 

93.5/95.9/  91.7/  89.0 

92.5  t  2.9 

7.2 

76.4/  75.7,  78.6 

76.9  ±  1.5 

7.4 

78.3/  73.3/  73.7/  76.4 

75.4  1  2.4 

10.1 

62.9,  63.1/  60.0 

61.9  t  1.8 

11.5 

44.3/  42.3,  45.9,  46.1 

44.7  i  1.8 

14.9 

21.0/  19.9,  19.9/  21.1 

20.4  ±  0.7 

19.6 

4.1/  4.0,  4.6,  4.2 

4.3  t  0.3 

25.4 

4.9,  0.8,  0.0 

1.9  t  2.6 

xaoie  r-2  Experimental  Values  of  Aerosol  Penetration. 
Test  Case  2.  Wind  Speed  «  3  m/s,  Inlet  Oriented 
Perpendicular  to  Air  Stream  (90*)/  Sampled  Flow  Rate  - 
130  L/min. 


A£D 

Particle 

Size/ 

urn  Penetration  Values/  Percent 


Kean  Penetration 
t  1  standard 
Deviation/ 
Percent 


3.5 

92.5, 

92.8/ 

91.7 

92.4 

t 

0.6 

4.9 

92.4/ 

91.3/ 

91.4 

91.7 

1 

0.6 

7.2 

80.2, 

84.5/ 

84.0 

82.9 

± 

2.4 

7.4 

80.6, 

82.7, 

83.4,  80.2 

81.8 

i 

1.6 

10.1 

64.7, 

65.9/ 

61.4/  61.2 

63.3 

t 

2.4 

11.5 

53,1/ 

58.5/ 

56.4 

56.0 

i 

2.7 

14.9 

24.5/ 

25.8, 

29.8/  26.8 

26.8 

♦ 

2.3 

19.6 

3.8, 

4.3,  4. 

2/  3.9 

4.0 

t 

0.3 

25.4 

0.5, 

0.0,  0. 

7,  0.9 

0.5 

t 

0.4 
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Table  F-3  Experiaental  Values  of  Aerosol  Penetration. 
Test  Case  5.  wind  Speed  »  3  n/S|  Inlet  Oriented 

Parallel  to  the  Air  Strean  (0*),  Saopled  Plow  Rate  » 
70  L/fflin. 

AED  Mean  Penetration 

Particle  t  i  standard 

Size,  Deviation, 

um  Penetration  Values,  Percent  Percent 


3.0 

100.3, 

113. 

6,  97.4,  110.1 

105.3 

i 

7.7 

5.1 

95.3, 

99.5, 

95.9,  94.7 

96.3 

t 

2.2 

7.0 

90.8, 

89.6, 

90.5,  88.8 

89.9 

0.9 

9.8 

83.3, 

87.5, 

80.3,  81.0 

83.0 

t 

3.2 

12.0 

73.7, 

74.1, 

72.4,  74.4 

73.7 

± 

0.9 

15.2 

55.6, 

56.8, 

51.8,  58.8 

55.8 

i 

2.9 

19.4 

29.3, 

34.1, 

35.2,  32.6 

32.8 

t 

2.6 

25.7 

7.1,  7 

.3,  5 

.4,  6.1 

6.5 

♦ 

0.9 

Table  F-4  Experiaental  Values  of  Aerosol  Penetration. 
Test  Case  6.  wind  Speed  »  3  a/s,  Inlet  Oriented 
Parallel  to  the  Air  stream  (0*),  Sampled  Flow  Rate  « 
130  L/ain. 


AED 

Particle 

Size, 


Mean  Penetration 
i  1  Standard 


urn 

Penetration 

Values 

,  Percent 

Percent 

w 

3.0 

94.8,  91.5, 

96.4, 

89.8 

93.1 

t 

3.0 

5.1 

92.9,  89.9, 

91.2, 

87.9 

90.4 

t 

2.1 

7.0 

92.4,  83.5, 

86.3, 

81.7 

86.4 

t 

5.5 

9.8 

69,4,  70.7, 

69.4, 

69.1 

69.6 

± 

0.7 

12.0 

56.8,  59.8, 

55.3, 

58.1 

57.6 

t 

1.9 

15.2 

31.1,  32.7, 

30.7 

31.5 

t 

1.0 

19.4 

11.9,  11.5, 

9.8 

11.1 

± 

1.1 

25.7 

0.8,  1.4,  1 

.2 

1.1 

t 

0.3 

149 


